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CHAPTER  1 


INTRODUCTION 


Integer  programming  problems  in  which  the  variables  are  assigned 
the  values  0  or  1  only  are  called  0-1  programming  problems.   The  implicit 
enumeration  method  is  known  to  be  an  efficient  algorithm  for  solving  0-1 
problems.   ILLIP  (iLLinois  Integer  Programming  Code)  was  developed  by 
T.  Ibaraki,  C.R.  Baugh,  T.K.  Liu  and  S.  Muroga  [k]   to  solve  general  0-1 
problems.   A  programming  manual  for  ILLIP  (which  was  coded  by  T.K.  Liu  and 
C.R.  Baugh)  was  written  by  T.K.  Liu  [6].   The  problem  of  designing  optimal 
networks  was  successfully  solved  by  S.  Muroga,  T.  Ibaraki,  C.R.  Baugh  and 
T.K.  Liu  [l],[i|],[5],[8],[9],[lO],[ll]  by  tailoring  ILLIP  to  the  network  syn- 
thesis problem.   This  approach  allows  optimal  networks  to  be  designed  under 
arbitrary  constraints  within  a  reasonable  computation  time.   This  report  is  a 
programming  manual  for  ILLODIE-AIF  (iLlinois  LOgical  Design  by  Implicit 
Eniimeration  using  the  All- interconnection  Inequality  Formulation)  -  a  FORTRAN 
IV  computer  program  (coded  by  T.K.  Liu  and  the  author)  which  uses  the  basic 
algorithm  of  ILLIP  but  is  designed  specifically  for  solving  network  synthesis 
problems.   Details  of  the  algorithm  used  in  ILLODIE-AIF  can  be  found  in  [l],[U] 

Let  all  of  the  gates  in  a  network  be  arranged  in  a  line  with  the 
output  gate  (one  of  the  output  gates  in  the  case  of  a  multiple  output  net- 
work) on  the  far  right.   If  each  gate  in  the  network. can  have  as  its  inputs 

only  the  outputs  of  the  preceeding  gates  on  its  left  as  well  as  the  connec- 

f 
tlons  of  the  external  variables,  then  the  network  is  called  a  feed-forward 


The  external  variables  are  the  variables  of  the  switching  function  which 
the  network  realizes  at  its  output.   The  term  "variables"  without  the 
adjective  "external"  refers  to  the  unknowns  in  the  inequality  formulation 
of  the  network  synthesis  problem  (see  Chapter  2). 


network.   When  a  function  is  given  to  be  realized,  ILLODIE-AIF  synthesizes 
feed-forward  networks  using  the  all- interconnect ion  inequality  formulation 
which  permits  an  interconnection  between  every  pair  of  gates  (see  [l]  for 
details  of  the  all- interconnect ion  formulation).   Since  this  formulation 
could  easily  derive  networks  with  loops,  a  special  subroutine  to  prohibit 
loops  has  been  incorporated  into  the  program.   The  use  of  the  all-inter- 
connection inequality  formulation  reduces  the  computation  time  by  precluding 
a  large  number  of  partial  solutions  which  are  equivalent  to  each  other  except 
for  permutations  of  the  gates  [l] . 

One  attractive  feature  of  ILLODIE-AIF  is  its  flexibility  which  is 
achieved  through  the  inequality  formulation.   First,  network  synthesis  prob- 
lems with  different  gate  types  can  be  solved  by  using  inequalities  which 
characterize  networks  of  the  desired  types.   ILLODIE-AIF  can  synthesize  net- 
works with  only  one  gate  type  (e.g.,  NOR  networks)  as  well  as  networks  with 
mixed  gate  types  (e.g.,  networks  with  a,  mixture  of  AND  gates  and  OR  gates). 
The  notation  ANP/OR  will  be  used  to  denote  networks  with  a  mixture  of  AND 
gates  and  OR  gates.   Networks  with  a  mixture  of  gate  types  should  not  be 
confused  with  networks  such  as  ECL  networks  in  which  the  output  gate  realizes 
both  the  NOR  and  OR  functions.   Second,  arbitrary  constraints  can  be  specified 
for  the  networks.   For  example,  a  maximum  value  for  the  fan- in  and/or  fan- 
out  of  each  gate  can  be  specified,  or  the  maximum  number  of  levels  in  the 
networks  can  be  restricted.   These  arbitrary  constraints  are  usually  in- 
corporated into  the  program  as  additional  inequalities.   However,  it  is 
sometimes  convenient  to  implement  restrictions  on  the  networks  by  writing 
a  special  subroutine.   For  example,  inequalities  could  be  included  to  prevent 
loops  in  the  networks  synthesized  by  ILLODIE-AIF,  but  it  is  more  practical 


to  use  the  special  subroutine  mentioned  above. 

The  user  of  ILLODIE-AIF  specifies  the  number  of  gates  the  networks 
are  to  contain.   This  value  remains  constant  throughout  the  computation  and 
ILLODIE-AIF  finds  all  the  networks  which  have  a  minimum  number  of  inter- 
connections for  the  specified  number  of  gates  subject  to  any  arbitrary  con- 
straints provided  in  the  additional  inequalities.   The  usual  design  objective 
is  to  find  all  networks  for  a  given  function  which  have  first  the  minimum 
number  of  gates  and  second  the  minimum  number  of  interconnections.   To  achieve 
this  objective  using  ILLODIE-AIF,  the  user  must  eventually  specify  the  min- 
imum value  for  the  number  of  gates.   A  procedure  for  determining  this  min- 
imum value  for  the  number  of  gates  is  discussed  in  section  2.k. 

The  network  synthesis  problem  has  the  following  characteristics: 
the  number  of  inequalities  is  larger  than  the  number  of  variables;  the 
coefficient  matrix  of  the  inequalities  is  very  sparse;  and  most  of  the  non- 
zero coefficients  are  -1  or  1.   For  example,  for  an  8  gate  NOR  network  which 
is  to  realize  a  function  of  three  external  variables,  the  number  of  inequal- 
ities is  1,5^6,  the  number  of  variables  is  592  and  the  number  of  non-zero 
coefficients  is  9^683  (1"^).   In  general,  as  the  number  of  gates  and/or 
external  variables  increases,  the  difference  between  the  number  of  inequal- 
ities and  the  number  of  variables  increases  and  the  percentage  of  non-zero 
coefficients  decreases. 

The  use  of  the  following  programming  gimmicks  helps  to  improve  the 

+ 
speed  of  the  program.    The  non-zero  coefficients  of  the  inequalities  are 

stored  in  memory  both  rowwise  and  columnwise  so  that  each  coefficient  can 


Pseudo-underlining  (discussed  in  [h] )  is  not  used  in  ILLODIE-AIF. 


be  fetched  efficiently.   The  concepts  of  the  maximum  sum  and  the  minimum 
sum  of  each  inequality  (see  Chapter  3)  are  used  to  improve  the  rate  of  con- 
vergence of  the  algorithm.   The  time  required  to  check  the  inequalities  for 
the  various  conditions  which  arise  during  the  computation  is  reduced  throiogh 
the  use  of  chain  checking  (see  Chapter  3)» 

Chapter  2  of  this  report  contains  algebraic  descriptions  of  the 
inequalities  which  characterize  the  networks  for  a  given  switching  function. 
It  also  discusses  the  integer  programming  formulation  of  the  network  synthesis 
problem  and  establishes  the  correspondence  between  the  integer  programming 
formulation  and  the  algebraic  inequalities.   Chapter  3  gives  a  summary  of 
the  functions  of  the  main  subroutines  in  the  program. 

Chapter  k   describes  GEN-IEQ,  the  subroutine  which  generates  the 
inequalities.   Network  synthesis  problems  with  different  gate  types  can  be 
solved  using  ILLODIE-AIF  by  simply  changing  the  inequalities.   Four  different 
versions  of  GEN-IEQ,  have  been  developed  which  allows  the  user  a  degree  of 
flexibility  in  that  NOR,  NOR/NAND,  NOR/AND,  and  AND/OR  network  synthesis 
problems  can  be  solved  by  using  the  appropriate  subroutine.   In  each  case 
additional  inequalities  which  reflect  special  properties  of  the  gate  types 
being  used  are  also  generated.   The  use  of  additional  inequalities  increases 
the  complexity  of  GEN-IEQ,  but  it  also  reduces  the  computation  time  [71- 

Chapter  5  contains  a  discussion  of  AGMT-VAR,  the  subroutine  which 
assigns  a  value  of  0  or  1  to  a  variable  (a  free  variable  in  the  terminology 
of  the  implicit  enumeration  method).   The  method  used  in  AGMT-VAR  is  a 
modified  version  of  the  method  described  by  Davidson  [2],[13]. 

Chapters  6  and  7  provide  the  necessary  information  for  anyone  who 
wants  to  use  ILLODIE-AIF  to  solve  network  synthesis  problems.   Chapter  6 


gives  a  complete  description  of  the  input  data  cards  for  the  program.   Both 
the  format  and  the  interpretation  of  the  printed  and  punched  output  from  the 
program  are  discussed  in  Chapter  7' 

Appendix  A  gives  a  complete  example  and  Appendix  C  gives  a  thorough 
discussion  of  the  additional  inequalities  in  the  current  version  of  the  program. 
Appendix  B  contains  a  summary  of  the  parameters  and  variables  used  in  the 
algebraic  descriptions  of  the  inequalities  and  additional  inequalities  given 
in  Chapter  2,  Chapter  k,   Appendix  A  and  Appendix  C. 

Table  1.1  shovrs  the  number  of  FORTRAN  instructions  in  ILLODIE-AIF 
for  each  of  the  four  versions  of  GEN-IEQ  as  well  as  the  amount  of  core  storage 
required  to  run  the  program. 


VERSION  OF 

NUMBER  OF  INSTRUCTIONS 

BYTES  OF 

CORE  STORAGE 

STORED 

GEN-IEQ 

GEN-IEQ   TOTAL 

PROGRAM 

DATA    TOTAL 

NOR 

327      2135 

67K 

266k    33 3K 

nor/and 

373     2181 

74k 

266K    340K 

nor/nand 

^06      2214 

74k 

266k    34ok 

and/or 

29k               2102 

72K 

266k    338k 

Table  1.1.   The  Number  of  Instructions  and  Core 
Requirement  for  ILLODIE-AIF. 


CHAPTER  2 


PROBLEM  DESCRIPTION 


2.1.   Definitions 

Suppose  a  network  has  R  gates  and  n  external  variables  x, ,  Xp,  . . . , 

X  .   ILLODIE-AIF  was  written  to  design  networks  which  have  o.ompletely  speci- 

+  — » 

fied  output  functions,  and  thus  all  of  the  2  possible  input  vectors  x  = 

(x^ ,  X   ,    ...,   X  )  must  be  considered.   Each  of  the  input  vectors  x  is  numbered 

as  x^'^^  (j-  1,  2,  3,  ...,  r  =  2^)  from  x^"'"^  =  (0,  0,  ...,  O)  through  x^^^  = 
(l,  1^  •-')   l)'   Let  w  represent  the  connection  from  external  variable  x 

to  gate  k  and,  if  the  complements  of  the  external  variables  are  available, 

k  —  k 

let  V  represent  the  connection  from  external  variable  x,  to  gate  k.   If  w. 

(or  V, )  =1,  the  connection  exists.   If  w  (or  v , )  =0,  the  connection  does 

not  exist.   Let  the  interconnection  from  the  i-th  gate  to  the  k-th  gate  be 

denoted  by  a.,  .   The  interconnection  exists  if  a.,  =1  and  does  not  exist  if 
^      ik  ik 

a.   =0.   Let  P.   denote  the  output  value  of  the  i-th  gate  when  the  external 
variables  connected  to  the  network  assume  the  values  of  input  vector  x^"^  . 
Finally,  let  P-^  =  ol       ^^K      (p|^^  =  1  when  the  output  of  gate  i  is  1  for 

input  vector  x  and  there  is  a  connection  from  gate  i  to  gate  k.  3-^  =  0 
when  either  the  output  of  gate  i  is  0  for  input  vector  x  or  gate  i  is  not 
connected  to  gate  k  or  both. ) 


A  summary  of  the  notations  used  in  this  paper  is  given  in  Appendix  B. 

It  is  possible,  though  somewhat  cumbersome,  to  change  ILLODIE-AIF  so  that 
networks  which  have  incompletely  specified  output  functions  can  be  designed. 


The  variables  defined  above  are  sufficient  for  describing  the  in- 
equalities which  characterize  networks  which  consist  of  only  one  gate  type 
(e.g.,  NOR  gate  networks).   However,  to  describe  the  inequalities  for  net- 
works which  consist  of  two  types  of  gates,  two  additional  definitions  are 
necessary.   First,  let  A.,  represent  the  type  of  the  i-th  gate.   For  NOR/AND 
and  NOR/NAKD  networks  X..  =  1  if  the  i-th  gate  is  a  NOR  gate  and  X.    =  0  if 
the  i-th  gate  is  an  AND  or  NAND  gate.   For  AND/OR  networks  X.  =  1  if  gate  i 
is  an  OR  gate  and  X.    =  0  if  gate  i  is  an  AND  gate.   Second,  let  7.P   =  cc 

(l-P.   ).   The  motivation  for  defining  the  7.,   is  given  in  section  2.2. 
^   1   ^  ik 

The  basic  configuration  of  a  feed -forward  network  which  realizes 
some  switching  function  f(x)  is  shown  in  Fig.  2.1.1  in  terms  of  the  w,  v  and 
a  variables  defined  above.   Networks  with  loops  cannot  be  derived  using  this 
configuration  because  there  are  no  a  variables  to  represent  the  connections 
from  the  output  of  each  gate  to  the  gates  on  its  left.  Fig.  2ol.2  shows 
the  all-interconnection  formulation  of  a  network  in  terms  of  the  w,  v  and 
Oi  variables.  In  this  formulation  there  is  an  ex  variable  to  represent 
the  connection  from  each  gate  (including  the  output  gate)  to  all  of  the 
other  gates.  ILLODIE-AIF  uses  the  all-interconnection  formulation  but 
derives  only  feed-forward  networks  by  using  a  subroutine  called  PHTLP  (see 
Chapter  3)  to  prohibit  any  loops  which  could  be  formed.  The  a  variables 
which  represent  the  connections  from  the  output  gate  to  all  of  the  other  gates 
in  the  network  could  be  eliminated  from  the  all -interconnect ion  formulation 
since  PHTLP  will  automatically  set  them  to  zero.  However,  these  o;  variables 
are  included  for  uniformity.  Similarly,  the  P,  P  and  7  variables  correspond- 
ing to  the  output  gate  could  be  eliminated,  but  are  included  for  uniformity. 


8 


Figure  2.1.1.      A  Feed-forward  Network 
(Gates   3  through  R-2  are  not   explicitly  shown. ) 
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Table  2.1.1  shows  the  range  of  values  for  the  subscripts  and  superscripts 
of  the  variables  defined  above  for  both  the  basic  feed-forward  configuration 
and  the  all-interconnection  formulation. 

2.2.   The  Basic  Inequalities 

The  inequalities  discussed  in  this  section  are  called  basic  inequal- 
ities because  they  completely  characterize  the  networks  which  are  to  be  designed. 
A  complete  description  of  the  basic  inequalities  for  a  network  with  R  NOR  gates 
is  given  below.   Also,  some  of  the  basic  inequalities  for  a  NOR/NAND  network 
are  used  to  explain  the  motivation  for  the  7  variables  defined  in  section  2.1. 
Complete  discussions  of  the  basic  inequalities  for  NOR/AM),  NOR/NAND,  and  AND/ 
OR  networks  are  given  in  [l],[5]  and  [12],  respectively. 

Using  the  definitions  given  in  section  2.1,  inequalities  characterizing 
all- interconnection  formulation  networks  with  R  NOR  gates  are  described  as 
follows  (assuming  the  complements  of  the  external  variables  are  not  available): 


.  T   t  t     .  T  "^ik  -       k 
t=l  1=1 


:/k 


(2.1) 


-  Z  w^  xi^^   -     Z  ^^.P   >0  -   U(l-R^^*b 
,  T   t  t     .  ,   ik  —         k 
t=l  1=1 


k  =  1,  2,  3,  •  •  •  >  R 


n 


J  ~  -L^  <— f     j,  •••^  ■^    ~    '—    • 
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Here  U  is  a  positive  number  large  enough  so  that  the  upper  inequality  is  non- 
restrictive  (always  satisfied)  if  P}   =  1  and  the  lower  inequality  is  non- 

K. 

restrictive  if  p}        =   0.   The  first  summation  in  each  inequality  of  (2.1) 
gives  the  number  of  external  variables  which  are  connected  to  gate  k  and  are 
equal  to  1  for  input  vector  x^*^  .   Likewise,  the  second  summation  counts  the 
number  of  gates  which  are  connected  to  gate  k  and  have  outputs  equal  to  1 
for  input  vector  x   .   If  P/;   =  0,  the  first  inequality  requires  that  the 
value  of  one  of  the  summations  must  be  greater  than  0.   (For  a  NOR  gate 
to  have  its  output  equal  to  0,  it  must  have  at  least  one  input  equal  to 


_,  the  second  inequality  requires  that  the  value  of  both  sum- 

K. 


1. )  If  p^^  ^  =  1, 

mations  be  0.   (A  NOR  gate  must  have  all  0  inputs  to  have  a  1  output.) 

Inequalities  (2.1)  characterize  all  of  the  gates  in  the  network. 
Of  course  the  values  of  P  '^   (j  =  1,  2,  3,  . . . ,  r  =  2  )  (i.e.,  the  values  of 
the  output  gate)  are  the  values  of  the  switching  function  f(x)  realized  by 
the  network.   Since  the  values  of  P    are  known,  it  is  actually  not  necessary 
to  include  both  inequalities  for  the  output  gate  for  each  value  of  j.   That 
is,  the  range  of  values  for  k  in  (2.1)  could  be  changed  to  k  =  2,  3j    '-•}    R 
and  the  output  gate  could  be  represented  by  the  inequality 


or 


n  R 

Z  w^  xp^  +  Z  ^[\^  >  1   if  f(x^'^'^)  =  0  {v[^^    =  0) 
t=l  ^  ^     i:.2  ^-^  ^ 

-  Z  wj  x^^^  -  Z  p[^^  >  0   if  f(x^'^'^)  =  1   (Pp^  =  1) 


for  j  =1,  2,  3,  ...,  2^.   Thus,  using  inequalities  (2.1)  for  the  output 
gate  includes  2  unnecessary  inequalities,  but  the  uniformity  of  (2.1)  makes 
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subroutine  GEN-IEQ,  simpler. 

To  complete  the  basic  inequalities  for  an  R  gate  NOR  network, 
the  non-linear  relation  3.!^  =  o;   P.    is  represented  by  the  following 

IK       iK   1 


linear  inequalities:: 


-Pr^  -  a.,    +  ^:P  >  -1 

X      ik    ik  — 

Fr^  +  a.,  -2^)p  >    0 

1      ik    ik  — 


(2.2) 


i  =1,  2,  3^  •  •  •  ^  R 

k  =  1,  2,   3,  ...,  R(i?^k) 
j  =  1,  ?,  3,  ...,  r  =  2"". 

The  following  example  illustrates  the  use  of  the  7.^  •   For  a 

IK 

NOR/NAND  network  one  set  of  inequalities  for  a  NAM)  gate  is  of  the  form 


n  R 

z  w^  (1-x^^^)  +  z  a.^  (i-ppb  >  1  -  u(i-p^^h  -  u  K^ 

t=l  i=l 

-  E  w^  (l-x|J')  -  E  a.^  d-Pp'))  >  0  -  UPj^J)  -  U  X^ 

t=l  1=1 

k  =.  1,  2,  3,  ...,  R 

3     ~   -^f    '-)     iy     '''}    r  =  c: 


(2.3) 


n 


Using  the  7.^  variables  defined  in  section  2.1,  these  inequalities  can  be 
rewritten  as 


n  R 

Z  w!^xp^+  Z  7^i^  >  1  -  U(l-P^^^)  -  U  A.  (2.U) 

t=l  "^  ^     i=l  ^^  ~         ^        ^ 


-A 


lU 


t=l  1=1 

k  =  1,  2,  3,  ...,  R 

i    =  1,    2,    3,    '..,    T   =  2^ 

since  xp^  =  (1-xf^^). 
t         t 

Inequalities  (2.3)  could  be  implemented  directly.   However,  in  the  subroutine 
AGMT-VAR,  which  is  discussed  in  Chapters  3  and  5^  it  is  convenient  to  have 
all  the  inequalities  in  the  form  of  (2.1).   For  inequalities  (2.1)  AGMT-VAR 
searches  for  aw,  v  or  3  variable  which  has  not  been  assigned  a  value  and 
assigns  a  value  to  the  variable  it  selects.  Using  inequalities  (2.U)  rather 
than  inequalities  (2. 3) allows  AGMT-VAR  to  treat  the  7   variables  in  the  same 
manner  as  the  3  variables  of  (2.1). 

2.3.   Additional  Inequalities 

Unlike  the  basic  inequalities  of  section  2.2,  the  additional  in- 
equalities described  below  are  not  essential  for  characterizing  a  network o 
Instead,  additional  inequalities  reflect  special  properties  of  certain  net- 
work configurations  as  well  as  special  properties  of  the  gate  types  being 
used.   These  inequalities  are  usually  effective  in  reducing  the  computation 
time  [7].  It  should  be  noted  that  it  may  be  necessary  to  modify  or  delete 
some  of  the  additional  inequalities  described  in  this  section  if  other 
additional  inequalities  such  as  those  for  fan- in  and /or  fan-out  restrictions 
are  added  (see   the  discussion  at  the  end  of  Appendix  C). 

The  following  discussion  describes  the  additional  inequalities  for 
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a  single  output  NOR  network  with  R  gates  where  gate  1  is  the  output  gate 
[l].   The  network  is  to  realize  a  function  of  n  external  variables  (x  ,  x  , 
...,  X  ),  and  it  is  assumed  that  the  complements  of  the  external  variables 
are  not  available.   In  each  case  the  verbal  description  is  followed  by  an 
algebraic  description  which  is  given  in  terms  of  the  w.  and  a.,  « 

(l)  Each  NOR  gate  except  the  output  gate  must  have  at  least  one 

input  from  the  external  variables  or  at  least  two  inputs  from 
the  other  gates. 


n       R 
2  E  W ,  +  Z  a.,  >  2 

t=l       1=1 


k  =  2,  3, 


,   R 


i^k 


(2)  If  three  NOR  gates  in  a  network  are  connected  as  shown  in 
Fig.  2,3.1>  then  it  can  be  proven  that  the  network  is  not 
optimal.   Inequalities  which  prohibit  triangular  connections 


no  other  output  connections 


Figure  2.3.1.   Triangular  Connection  of  Three  NOR  Gates, 


^  A  complete  discussion  of  the  additional  inequalities  for  NOR,  NOR/AND, 
NOR/NAND  and  AND/OR  networks  is  given  in  Appendix  C. 


16 


are  given  by 


i 

= 

1, 

2, 

.  .  .  , 

R 

J 

= 

1, 

2, 

'  '  '  > 

R 

k 

= 

1, 

2, 

'  '  '  y 

R 

i 

^ 

J 

^k. 

This  property  is  true  even  if  NOR  gate  i  is   replaced  by 


external  variable  x. . 

1 

-"i  -  ^i  -  ^Ok  ^  £  ^Ji 

> 

-2 

i/k 

J    —  ^}   '- )    •  • ' )    R 

k   =  1,    2,    ...,    R 

i   -  1,    2,    ...,   n 

(3)     All  NOR  gates   except  the  output  gate   connect  to  at  least 
one  other  gate. 


R 

Z     Q;  .   >1         k  ^  2,   3,    ...,   R 
3=1     ^ 


(h)     All  NOR  gates  which  are  connected  to  the  output  gat( 
are  not  connected  to  any  other  gate. 
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R 

-U  a. ^  -  z  a   >  -U   i  =  2,  3,  . . . ,  R 
^^      k=2  ^^  ~ 

Mi 


U  is  a  positive  integer  large  enough  so  that  the  inequality 

becomes  non-restrictive  if  a. ^  =  0. 

il 


2.U.   The  Network  Synthesis  Problem  as  an  Integer  Programming  Problem 

Any  0-1  programming  problem  can  be  defined  as: 

minimize  c  z 
subject  to  b  +  A  z  >  0.  (2.5) 

Here  c  =  (c^ ,  c   ,    ...,  c  )  is  an  N-dimensional  row  vector;  z  =  (z  ,  z  ,  ••., 

z,^)  is  an  N-dimensional  column  vector;  b  =  (b, ,  b^,  •  •  •  »  b  )  is  an  m-dimen- 
N  J      \  jL'   2'     '   m 

sional  column  vector;  A  =  [a. .]  is  an  m  x  N  matrix;  and  0  is  the  m-dimensional 
column  vector  whose  entries  are  all  zero.   The  variables  z  for  t  =1,  2,  •••, 
N  can  assume  the  values  0  or  1  only. 

Any  assignment  of  0  or  1  to  all  of  the  variables  z.  is  called  a 
solution,  and  any  solution  that  satisfies  b  +  A  z  >  0  is  called  a  feasible 
solution.   Otherwise  a  solution  is  said  to  be  infeasible.   An  optimal  solution 
is  a  feasible  solution  which  minimizes  c  z.   If  no  feasible  solution  exists, 
the  problem  is  called  infeasible.   Any  assignment  of  0  or  1  to  part  of  the 
variables  is  called  a  partial  solution.   The  variables  which  have  been  set  to 
0  or  1  in  a  partial  solution  are  called  fixed  variables  while  the  remaining 
variables  are  called  free  variables . 

For  the  network  synthesis  problem  each  z  variable  uniquely  represents 
a  w,  V,  a,   V,    p,  X.  or  7  variable  (defined  in  section  2.1).   The  objective 
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function  c  z  represents  the  interconnections  in  the  network.   The  intercon- 
nections include  both  the  connections  from  the  external  variables  to  the  gates 
and  the  connections  among  gates  and  can  be  represented  by 

R   n    .    R   n    .    R   R 

Z  L     wi   +     L       L     vi   +     L       Z  a  (2.6) 

j=l  t=l      0=1  t=l  ^   i=l  k=l   "^ 

Mi 

for  a  network  with  n  external  variables  and  R  gates.   Values  can  be  assigned 

to  the  coefficients  c.  of  the  objective  function  c  z  by  comparing  the  inner 

product  c  z  with  expression  (2.6).  That  is,  c.  =  1  if  z.  represents  a  w,  v 

or  a  variable  and  c.  =  0  if  z.  represents  a  P,  p,  A.  or  7  variable.   Finally, 

the  basic  and  additional  inequalities  described  in  sections  2.2  and  2.3, 

respectively,  are  represented  by  expression  (2.5)-   The  coefficients  of  the 

variables  and  the  constants  in  the  expressions  for  the  basic  and  additional 

inequalities  determine  the  values  of  the  a. .  and  b  in  (2.5). 

ij      k 

ILLODIE-AIF  is  usually  used  to  obtain  networks  which  have  first 
a  minimum  number  of  gates  and  second  a  minimum  number  of  interconnections. 
By  minimizing  c  z  under  the  constraints  b  +  A  z  >  0  which  are  formed  under 
the  assumption  that  the  networks  have  R  gates,  ILLODIE-AIF  obtains  the  net- 
work (or  networks)  with  R  gates  which  have  a  minimum  number  of  interconnec- 
tions assuming  a  feasible  solution  exists.   Thus,  to  achieve  the  desired 
objective,  it  is  necessary  to  find  the  smallest  value  of  R  for  which  a  feasi- 
ble solution  exists.   This  value  of  R  can  be  found  as  follows:   begin  ^-rlth 
an  initial  value  for  R  and  solve  (2.5).   If  there  is  no  feasible  solution, 
set  R  to  R  +  1  and  solve  (2.5)  again.   Repeat  this  procedure  until  a  feasible 
solution  is  obtained.   If  a  feasible  solution  is  found  when  (2.5)  is  solved 
with  the  initial  value  of  R,  then  set  R  to  R-1  and  solve  (2.5)  again.   Repeat 
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this  procedure  until  (2.5)  becomes  infeasible.   Since  the  computation  time  for 
ILLODIE-AIF  increases  exponentially  (roughly  speaking)  as  R  increases,  it  is  in 
the  best  interest  of  the  user  to  specify  an  initial  value  of  R  which  is  too  small 
rather  than  too  large.   That  is,  if  R  is  the  smallest  value  of  R  for  which  a 
feasible  solution  exists,  the  computation  time  for  values  of  R  smaller  than  R 
is  much  shorter  than  the  computation  time  for  R  ,  and  the  computation  time  for 
values  of  R  greater  than  R^  is  much  longer  than  the  computation  time  for  R  . 

The  manner  in  which  the  w,  v,  a,  P,  p,  \  and  7  variables  are  re- 
presented by  the  z  variables  is  indicated  by  the  following  example.   Assume 
the  network  to  be  synthesized  is  to  realize  a  two  variable  function  f(x). 
Assume  also  that  the  complements  of  the  external  variables  are  available  as 
inputs  to  the  network,  that  the  network  can  have  two  gate  types  (e.g.,  NOR/AND), 
and  that  R  =  3  is  being  used  as  an  initial  value.   Table  2.U.1  shows  the  sub- 
scripts i  of  the  variables  z.  which  represent  the  w,  v,  a,  P,  3,  \   and  7 
variables.   For  this  example  N  =  8I.   If  the  assumption  that  the  complements 
of  the  external  variables  are  available  is  removed.  Table  2.^.1  (b)  is  elim- 
inated and  the  subscripts  of  the  z.  are  changed  accordingly  (3-,p  becomes  z  , 
a       becomes  z   ,  etc.;  and  N  =  75).   On  the  other  hand,  if  the  network  is  to 
have  only  one  gate  type  (e.g.,  NOR),  both  (f)  and  (g)  are  eliminated  from 
Table  2.4.1  and  N  =:  ^h.      Finally,  if  the  complements  of  the  external  variables 
are  not  available  and  the  network  is  to  have  only  one  gate  type,  (b),  (f)  and 
(g)  are  eliminated  from  Table  2.U.1;  3   becomes  z  ,  a   becomes  z      ,   etc.; 
and  N  -  kQ. 
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(a) 


(b) 


X"^ 

\ 

1 

2 

1 

1 

2 

2 

3 

J+ 

< 

5 

6 

X"^ 

\ 

1 

2 

1 

7 

8 

2 

9 

10 

^^ 

11 

12 

(c) 


(d) 


1 

1 

2 

3 

h 

^12 

13 

lu 

15 

16 

o(j) 

^13 

17 

18 

19 

20 

^21 

21 

22 

23 

21+ 

o(j) 

^23 

25 

26 

27 

28 

^31 

29 

30 

31 

32 

^32 

33 

3h 

35 

36 

1 

1 

2 

3 

^ 

■)f* 

37 

38 

«2k 

39 

^Hf 

1+0 

^3k 

i+1 

1+2 

■x-x- 

Table  2.1+.1.  The  Subscripts  of  the  z  Variables 
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(e) 


X 

1 

2 

3 

h 

p(j) 

^1 

^3 

kk 

^5 

kG 

2 

hi 

kQ 

h9 

50 

p(j) 

3 

51 

52 

53 

5i+ 

(f) 


X 

1 

2 

3 

^. 

55 

56 

57 

(g) 


\ 

1 

2 

3 

k 

^12 

58 

59 

60 

61 

^13 

62 

63 

6i| 

65 

^21 

66 

67 

68 

69 

^23 

70 

71 

72 

73 

^31 

7^+ 

75 

76 

77 

^32 

78 

79 

80 

81 

Table  2.4.1.      The   Subscripts  of  the  z  Variables    (continued) 
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CHAPTER  3 

PROGRAM  DESCRIPTION 

The  program  for  ILLODIE-AIF  consists  of  twenty  seven  subroutines, 
but  they  can  be  grouped  logically  into  the  following  seven  main  subroutines: 
INIT  (iNITialize),  GEN-IEQ  (GENerate  -  InEQualities),  CHK-IEQ  (CHecK  -  In- 
Egualities),  AGMT-VAR  (AuGMenT  -  VARiable),  BCKTRK  (BaCKTRacK),  CHNGMM  (CHaNGe 
Maximum  sum  and  Minimum  sum),  and  PHTLP  (ProHibiT  LooP).   Fig.  3*1  shows  the 
linkage  between  these  subroutines.   Arrows  which  originate  and  terminate 
at  the  same  subroutine  represent  FORTRAN  IV  subroutine  CALL  statements.   That 
is,  subroutine  CHNGMM  is  called  by  subroutines  AGMT-VAR,  CHK-IEQ,  and  BCKTRK 
and  will  return  to  the  subroutine  from  which  it  was  called  when  its  computa- 
tion is  completed.   Similarly,  PHTLP  is  called  by  CHNGMM.   Arrows  which  orig- 
inate at  one  subroutine  and  terminate  at  another  represent  a  more  general 
kind  of  program  flow.   For  example,  the  program  flow  between  subroutines 
CHK-IEQ,  AGMT-VAR  and  BCKTRK  is  controlled  by  a  subroutine  called  GEOF  (not 
shown  in  Fig.  3.1)'   The  remainder  of  this  chapter  contains  a  brief  descrip- 
tion of  each  of  the  seven  main  subroutines. 

1.   INIT 

This  subroutine  initializes  several  program  parameters  after 
reading  the  problem  parameter  card  (see  Chapter  6).   Then  INIT  does  one  of 
the  following  (depending  upon  the  value  of  one  of  the  parameters): 

(a)   If  the  program  is  to  generate  the  objective  function  and  the 
inequalities,  INIT  transfers  to  GEN-IEQ.   GEN-IEQ  generates 
the  objective  function,  and  GEN-IEQ  and  INIT  work  together  to 
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GEN-IEQ 

1 1 
\  1 

INIT 

-- 

1 

i 

r 

I 

CHK-IEQ, 

^ 

/ 

i 

I 

A(^T-VAR 

v_y 

BCIdRK 

3       CHNGMM     C 

J 

I 

PHTLP 

Figure  3.1.   Linkage  Between  Main  Subroutines 
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generate  the  inequalities  (see  Chapter  h) . 
(b)   If  the  objective  function  and  inequalities  are  provided  on 
cards,  INIT  reads  the  objective  function  first  and  then  the 
inequalities.   The  inequalities  are  read  one  at  a  time, 
expanded  if  necessary  and  properly  stored. 

INIT  also  generates  the  initial  partial  solution  which  consists 
of  the  following  two  parts  (depending  on  what  the  user  provides): 

(a)  The  P  variables  for  the  output  gate  (P|^^  for  j  =1,  2,  3,  ..., 
2  )  which  are  assigned  values  according  to  the  values  of  the 
switching  function  f(x). 

(b)  Any  variables  to  which  the  user  wants  to  assign  specific  values 
(either  0  or  l). 

Part  (a)  is  always  included  in  the  initial  partial  solution  while  part  (b)  is 
optional.   All  variables  in  the  initial  partial  solution  are  underlined. 
If  a  variable  z  in  a  partial  solution  is  underlined,  then  the  partial 
solution  in  which  z,  is  assigned  the  opposite  value  need  not  be  considered. 

2.   GEN-IEQ 

Since  the  number  of  inequalities  in  most  network  synthesis  problems 

i 

is  large,   it  is  impractical  to  punch  the  inequalities  on  cards.   For  this 
reason  subroutines  called  GEN-IEQ,  were  developed  to  generate  the  inequalities. 


The  concept  of  underlining  was  first  discussed  in  [3]« 

Statistics  on  the  number  of  inequalities  for  NOR  and  NOR/ AND  network  synthe- 
sis problems  can  be  found  in  [l].   Similar  statistics  for  NOR/NAND  problems 
are  given  in  [5] . 
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Chapter  k   of  this  report  will  discuss  in  greater  detail  how  the  inequalities 
are  generated.   GEN-IEQ,  is  also  used  to  generate  the  objective  function,  cz. 

3.   CHK-IEQ 

Whenever  a  new  partial  solution  is  formed,  CHK-IEQ  will  check  the 
inequalities  to  determine  whether  any  free  variables  must  be  set  to  0  or  1 
in  order  to  satisfy  some  inequality.   If  any  such  free  variable  is  found, 
it  will  be  set  to  the  required  value  and  underlined.   If  a  fixed  variable  in 
a  partial  solution  is  underlined,  then  during  backtrack  the  opposite  value  of 
this  variable  need  not  be  assigned.   However,  it  is  usually  not  necessary  to 
check  all  of  the  inequalities  each  time  a  new  free  variable  is  fixed.   That 
is,  at  any  point  in  the  sequence  of  partial  solutions,  some  inequalities  may 
be  such  that  they  will  always  be  satisfied  regardless  of  which  free  variables 
are  fixed.   To  make  use  of  this  fact,  a  method  called  chain  checking  is  used 
based  on  the  array  CLSP.   CLSP(l)  consists  of  the  row  number  of  the  next  in- 
equality to  be  checked  after  the  I-th  inequality  is  checked.   If  the  I-th 
inequality  is  not  in  the  chain,  then  CLSP(l)  =  0.   Thus  the  non-zero  entries 
of  CLSP  form  a  chain  of  the  row  numbers  of  the  inequalities  which  must  be 
checked. 

The  remainder  of  this  section  describes  how  CHK-IEQ  checks  an  in- 
equality.  Assume  a  partial  solution,  F,  is  given  byF  =  (d.  ,  d.  ,  ...,  d.  } 

^1   ^2       ^k 

where  d.  is  the  binary  value  (0  or  l)  assigned  to  fixed  variable  z.  for  j  =  j  , 

j  ,  ...,  j  .   Some  of  the  z.  may  be  underlined  which  is  denoted  by  d . .   Let 
^       k  J  —J 

W  be  the  set  of  the  remaining  free  variables.   To  faciliate  the  discussion  of 
CHK-IEQ,,  the  following  three  vectors  are  defined. 
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(a)  Partial  sum:   y(F)  =  (y^(F)  ,  y2(F),  ...,  yj^)) 

where  y.  (f)  =  "b  +   Z   a.  .  d.  for  i  =  1,  2,  . . . ,  m. 

z.€F  ^^   ^ 
J 

The  notation   Z   means  that  the  summation  is  taken  over  all 

z.€F 
J 

fixed  variables. 

(b)  Maximum  sum:   u(f)  =  (u  (F),  u  (F),  ...,  u  (F)) 

where  u. (f)  =b.+  L       a..d.+  L       a.. 

^   z.eF  ^J  J   z.ew  ^J 
J  3 

a.  .>0 


y. (f)  +   Z    a. .  for  i  =  1,  2, 


m. 


Z.€W    ^«^ 

J 

a.  .>0 

The  notation    Z    means  that  the  summation  is  taken  over 

z.€W 

0 
a.  .>0 

all  free  variables  with  positive  coefficients, 
(c)   Minimum  sum:   1(F)  =  (i^(F),  i^iT) ,    ...,    I   (f)) 

where  i. (F)  =  y, (F)  +   Z    a. .  for  i  =  1,  2,  . . . ,  m. 

z.ew   ^ 

a.  .<0 

The  notation   Z    means  that  the  summation  is  taken  over 

z.ew 

J 
a.  .<0 
ij 

all  free  variables  with  negative  coefficients. 

Thus  u. (F)  is  the  maximum  value  that  the  left  hand  side  of  the  i-th  inequali- 
ty can  achieve  for  all  solutions  which  contain  F.   Similarly,  i. (F)  is  the 
minimum  value. 

If  the  i-th  inequality  needs  to  be  checked  (i.e.,  CLSP(l)  ^   O), 


27 

then  CHK-IEQ,  will  do  one  of  the  following  depending  on  the  value  of  u.  (F)  [h] . 

(a)  If  u. (F)  <  0,  then  no  feasible  solution  containing  F  exists. 
The  program  wi.lJL  backtrack  by  transferring  to  BCKTRK. 

(b)  If  u. (F)  -  0,  then  all  free  variables  with  negative  coefficients 
in  the  i-th  inequality  are  set  to  zero  and  underlined;   and  all 
free  variables  with  positive  coefficients  in  the  i-th  inequality 
are  set  to  one  and  underlined. 

(c)  If  u. (f)  >  0;  then,  for  all  free  variables  z.  which  have  co- 
efficients  satisfying  ja..]  >u.(F),  CHK-IEQ  will  either  set 

z.  to  one  (if  a. .  >  O)  or  set  z.  to  zero  (if  a. .  <  O)  and  under- 
line  it. 

In  (b)  and  (c)  above  the  program  also  compares  the  value  of  the  objective 
function  for  the  current  partial  solution,  F,  with  the  current  upper  bound 
on  c  z  .   If  the  bound  is  exceeded,  the  program  transfers  to  BCKTRK. 

The  vectors  u(F)  and  i(F)  must  be  updated  each  time  a  free  variable 
is  fixed.   If  u. (F)  changes  for  some  i,  i  is  added  to  the  end  of  the  chain 
of  inequalities  to  be  checked  because  the  i-th  inequality  must  be  checked. 
All  of  the  inequalities  in  the  chain  are  repeatedly  checked  until  one  of 
the  following  occurs: 

(a)  u. (f)  <  0  for  some  i.   In  this  case  there  can  be  no  feasible 
solutions  which  contain  F  so  the  program  transfers  to  BCKTRK. 

(b)  A  feasible  solution  is  found.   The  feasible  solution  is  print- 
ed only  if  the  value  of  the  objective  function,  c  z,  for  the 
feasible  solution  is  less  than  or  equal  to  ZMAK  where  ZMAX  is 
the  current  upper  bound  on  c  z  .   The  program  then  transfers  to 
BCKTRK. 
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(c)      No  more  free  variables  can  be  fixed  by  CHK-IEQ,  so  the  program 
enters  AGMT-VAR. 

k.      AGMT-VAR 

When  CHK-IEQ,  cannot  fix  any  more  free  variables,  the  program  enters 
AGMT-VAR.   AGMT-VAR  does  one  of  the  following: 

(a)  If  all  of  the  variables  are  fixed,  a  feasible  solution  has  been 
found.   As  before,  the  feasible  solution  will  be  printed  if  it 
gives  c  z  a  value  less  than  or  equal  to  ZMAX  where  ZMAX  is  the 
current  upper  bound  on  cz.  The  program  then  transfers  to  BCKTRK. 

(b)  One  free  variable,  say  z.,  is  fixed  according  to  the  scheme 
discussed  in  Chapter  5«   The  program  transfers  back  to  CHK-IEQ, 
to  check  the  inequalities  for  the  new  partial  solution  which  is 
the  previous  partial  solution  plus  z.. 

(c)  If  no  free  variable  can  be  fixed  by  the  scheme  discussed  in 
Chapter  5^  then  there  are  no  feasible  solutions  which  contain 
the  current  partial  solution.   The  program  transfers  to  BCKTRK. 

5.   BCKTRK 

When  AGMT-VAR  or  CHK-IEQ  either  finds  a  feasible  solution  or  dis- 
covers that  no  feasible  solution  exists  containing  the  current  partial  so- 
lution, the  program  enters  BCKTRK  in  order  to  generate  a  new  partial  solution. 
The  program  returns  to  CHK-IEQ  after  the  new  partial  solution  is  generated  by 
the  backtrack  scheme  described  below.  This  backtrack  scheme  guarantees  that 
all  possible  solutions  will  be  implicitly  enumerated  [3]. 

To  facilitate  the  description  of  the  backtrack  scheme  let  us 
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describe  how  a  partial  solution  is  stored  in  ILLODIE-AIF.   Two  arrays,  X  and 
XS,  are  used  to  represent  a  partial  solution.   X  indicates  what  value  has  been 
assigned  to  each  fixed  variable  and  XS  indicates  whether  or  not  each  fixed 
variable  has  been  underlined.   Initially  X  and  XS  are  empty.   When  the  I-th 
free  variable  is  fixed,  X(l)  and  XS(l)  are  assigned  the  corresponding  sub- 
script with  the  proper  sign.   X(l)  =  J  means  that  z  was  set  to  1  while  X(l)  = 
-J  means  that  z  was  set  to  0.   Similarly,  XS(l)  =  J  means  that  z  was  not 

underlined  while  XS(l)  =  -J  indicates  that  z  was  underlined.   For  example, 

J 

if  the  current  partial  solution  is  given  by  F  =  (3^-^^  2}),  then  X  =  (3^-^^  2) 
and  XS  =  {-3,  h,-2].      If  AGMT-VAR  or  CHK-IEQ  sets  z.   to  0  and  underlines  it, 
then  the  new  partial  solution  is  F  =  (_3^-^^  ^,-6}  whi.ch  is  represented  by 

X  =  (3,-^,  2,-6}    XS  =  {-3,  ^,-2,-6}  (3.1) 

In  BCKTRK  all  underlined  variables  to  the  right  of  the  right  most 
non-underlined  variable  jn  the  sequence  of  numerals  in  F  are  changed  to  free 
variables.   The  right  most  non-underlined  variable  is  then  set  to  its  opposite 
value  and  underlined  thus  forming  a  new  partial  solution.   For  example,  if 
the  current  partial  solution  is  given  by  (3.I),  then  BCKTRK  changes  z  and 
z^  to  free  variables.   The  right  most  non-underlined  variable,  z,  ,  is  then 
set  to  1  and  underlined.   The  new  partial  solution  is  given  by  F  =  {3,  ^] 
which  is  represented  as  X  =  {3,  ^)  and  XS  =  (-3^-^)« 

If,  upon  entering  BCKTRK,  all  the  fixed  variables  are  underlined, 
then  all  possible  solutions  have  been  implicitly  enumerated.   All  optimal 
solutions  to  the  problem  have  been  printed.   If  no  feasible  solution  has  been 
found,  the  given  problem  is  infeasible. 
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6.  CHNGMM 

CHNGMM  updates  u(F)  and  i(F)  each  time  a  free  variable  is  fixed. 
Since  free  variables  can  be  fixed  by  CHK-IEQ,  AGMT-VAR,  and  BCKTRK,  CHNGMM 
can  be  called  by  each  of  these  subroutines  and  will  return  to  the  subroutine 
which  called  it  when  its  computation  is  finished.   If  a  free  variable,  z., 
is  set  to  1,  CHNGMM  will  update  u(F)  and  ^(f)  by  setting 

u.  (F'  )  =   u.  (F)      and  H.  (F' )  =  i.  (F)  -l-a.  .      if  a.  .  >  0, 
1        1  1       1      ij         ij    ' 

u.  (F'  )  =  u.  (F)  +a.  .  and  H.  (F'  )  =  H .  (F)  if  a.  .  <  0, 

1       1      ij       1        1  ij    ' 

u. (F' )  =  u. (F)      and  l. (F' )  =  L  (F)  if  a. .  =  0, 

i  =1,  2,  ' '  • ,   m. 

If  z .  is  set  to  0,  then  CHNGMM  will  set 
J 

u.  (F* )  =  u.  (f)  -  a.  .  and  L  (F'  )  =  St.  (F)  if  a.  .  >  0, 

u.  (F'  )  =  u.  (F)      and  L{Y'  )   =  i.  (F)  -  a.  .  if  a.  .  <  0, 

u.(F')=u.  (F)      and  L{Y')^L(j)  if  a.  .  =  0, 

i  =1,  2,  •  •  • ;  ni. 

CHNGMM  has  two  other  functions.  First,  CHNGMM  assigns  values  to 
the  elements  of  the  array  CLSP.  Second,  CHNGMM  calls  PHTLP  to  ensure  that 
the  network  remains  loop  free. 

7.  PHTLP 

A  gate,  L,  is  called  a  successor  of  a  gate,  J,  (and  gate  J  is  called 
a  predecessor  of  gate  L)  if  and  only  if  there  exists  a  sequence  of  gates 
J  =  J_,  J  ,  J  ,  ...,  J  =  L  such  that  the  output  of  gate  J.  is  connected  to 
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gate  J,  ^  for  i  =0,  ....  s-1  (i.e.,  a^  _  =  a^  ^  =  . . .  =  a^   ^  =  l). 
^+1  Vl    V2         Js-l-^s 

Whenever  the  connection  variable,  ^j}   from  the  output  of  some  gate,  H,  to 
the  input  of  another  gate,  I,  is  set  to  1,  PHTLP  is  called  to  insure  that  the 
network  remains  loop  free.   PHTLP  stores  the  gate  numbers  of  all  the  succes- 
sors of  gate  I  (including  gate  l)  in  the  array  SCR  and  the  gate  numbers  of  all 
the  predecessors  of  gate  H  (including  gate  H)  in  the  array  PRD.   Then  the  a 
variables  representing  the  connections  from  each  gate  in  SCR  to  every  gate 
in  PRD  are  set  to  0. 
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CHAPTER  k 

generatinct  the  inequalities 

^.1.   Generating  the  Basic  Inequalities 

In  ILLODIE-AIF  the  subroutines  GEN-IEQ  and  INIT  are  used  to  generate 
the  inequalities.   How  these  two  subroutines  are  used  to  generate  the  basic 
inequalities  for  a  network  with  R  NOR  gates  and  n  external  variables  is  dis- 
cussed in  this  section.   The  same  method  is  used  to  generate  the  basic  inequal- 
ities for  and/or,  N0R/NAJ)JD  and  NOR/AND  networks. 

The  basic  inequalities  for  a  network  with  R  NOR  gates  and  n  external 
variables  (without  their  complements)  given  in  Chapter  2  are  repeated  here  for 
convenience.   The  R  NOR  gates  are  characterized  by  the  inequalities 


I     wSi^^  +  2  ^[P    >1-  up(^) 


(i+.l) 


t=l  ^  ^     i.l  ^^  -       k 
iA 

-  Z  wl^  xi^^  -   Z  ^^.P    >0  -   U(l-R^^^ 
t=l  ^   ^     i=l  ^^  -  ^ 

k  =  1,  2,  3,  ...,  H 

J  =  1,  2,  3,  ...,  r  =  2^^. 

As  before,  U  is  a  positive  number  large  enough  so  that  the  upper  inequality- 
is  non- restrictive  if  P^"^  "^  =  1  and  the  lower  inequality  is  non- restrictive  if 

Pp^  =  0.   The  relation  ^[^   =  <^-i  ?•   is  represented  by  the  inequalities 

k  ik     ik  i 
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■  p^^^  -  a.,  +  ^^P     >  -1 

1      ik    ik   — 

P?^^  4-  a.,  -  2  ^[p   >    0 
1      ik     ik  - 


(^.2) 


i 

=  1, 

2, 

3, 

. . . , 

R 

k 

=  1, 

2, 

3, 

. . . , 

R(i/^k) 

3 

=  1, 

2, 

3, 

. . . , 

r   =2^ 

n 
For  each  fixed  value  of  k,  there  are  2  inequalities  formed  from 

each  inequality  in  (U.l)  (as  j  goes  from  1  to  2  ).   Similarly,  for  each  com- 

n 
bination  of  values  for  i  and  k,  there  are  2  inequalities  formed  from  each 

inequality  in  {k.2).      Each  such  group  of  2  inequalities  will  be  referred  to 

as  a  block  of  inequalities.   As  an  example,  let  R  =  3  and  n  =  2.   The  first 

inequality  in  (^.l)  yields  the  following  block  of  inequalities  for  k=l: 


(It. 3) 


J  ,  3.  „1  43)  ,  „1  ,(3)  ,  ,(3)  ,   ,(3)  ,  „  p(3)  >  , 
,,k:     .1  4'^  .  .\  X  W  .  4^)  .  4^)  .  U  P(^)  >  1 


Similarly,  the  second  inequality  in  (^.l)  yields  the  following  block  of 
inequalities  for  k  =  1; 
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i   =  l:   -w^  x(l)  -  „^  x(l)  -  ,W   .  3(1)  .  „  ,(1)  ^  .„ 

,  =  .:  -.1   xf )  -  „^  4^)  -  ,(^)  .  4f  -  „  p(^)  >  .„ 

0  =  3:   -w^  xf)  .  „1  ,(3)  .  3(3)  .  3(3)  .  „  p(3)  >  .„ 

d  =  ^:   .„!  xf )  .  .1  x(^)  -  4J)  -  eW  .  „  p(^)  >  .u 

For  i  =  1  and  k  =  2,  the  first  inequality  in  (i+.2)  yields  the  following 
block  of  inequalities: 


J  =  1:  -i^'   -  ^2  ^  ^[2     >  -1 


(h.k) 


J  =  2:   .p|2^  -  a^2  +  ^[f   >  -1 


j  =  3:   -Pp^  -  a^2  +  ^[l^   >  -1 


0=^:    -p1^^  -  a,,  .  el^L> -1 


Finally,  for  i  =  1  and  k  =  2,  the  second  inequality  in  (1|.2)  yields  the 
following  block  of  inequalities: 


3    =   1:   pj^^  +  a^2  -  2PJ2^  >  0 


(i+.5) 


(2)  _   _  0.(2)  ^^ 

(h.6) 


J  =  2:   P^  ^  +  0^2  -  2Pi2^  >  0 


.(3)  ,  ^     0.(3) 
1 


j  =  3:   P;-"'  +  ^2  "  ^^12  -  ° 


(^)  ,  ^     oo(^) 


j  =  k:      P^^^  +  0^2  ■  ^^12  -  ^ 
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Now,  using  Table  2.U.1,  inequalities  (^.3),  (^-^),    (^+-5)  and  (4.6)  can  be 
rewritten  in  terms  of  the  variables  z, .   The  value  nXR+l=7is  used  for 
U.   Inequalities  (4.3)  become 


3=1:      Oz^  +  Oz^  +  Iz^^  +  Iz^^  +  Tz^^  >  1 


j  =  2:   Oz^  +  izg  +  Iz^^  +  Izg^  4-  7^38  >  1 


j  =  3:   Iz^  +  Ozg  +  Iz^^  +  Iz^^  +  Tz^^  >  1 


3   =  h:      Iz^  +  IZg  +  Iz^g  +  Izg^  +  Tz^Q  >  1 


Inequalities  (h.k-)   become 


J  =  1:   -Oz^  -  Oz^  -  lz^5  -  IZ23  -  TZ3^  >  -7 


j  =  2:   -Oz^  -  iz^  -  iz^g  -  iz^^  -  7Z38  >  -7 


j  =3:   -iz^-  OZ2-  iz^^-lz^^-  TZ3^^-T 


3    =  k:      -Iz^  -  IZ2  -  lz^8  -  1^26  -  ^\0  ^  -^ 
Inequalities  (4.5)  become 


j  =  1:   -1Z3^  -  1Z3^  +  Iz^  >-l. 


j   =  2:   -lz3g  -  1Z3^  +  Izg  >  -1 


j  =  3:  -1Z3^  -  1Z3^  +  iz^  >  -1 


j  =  h:      -Iz^Q  -  Iz^^  +  Iz^Q  >  -1 


(4.7) 


(4.8) 


(4.9) 
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Finally,  inequalities  (^.6)  become 


(i+.io) 


j  =  1:   Iz^^  +  Iz^^  -  2z^  >  0 


j  =  2:   Iz  g  +  Iz   -  2Zg  >  0 


j  =  3:  1^33  +  izj^  -  2zg  >  0 


J  =  "*=  I^UO  +  1"31  ■  2"l0  ^  ° 


Observing  the  patterns  of  both  the  coefficients  and  subscripts  of 
the  variables,  z   ,    in  the  blocks  of  basic  inequalities  leads  to  the  definitions 
given  below.  These  definitions  are  important  for  understanding  how  GEN-IEQ 
and  INIT  are  used  to  generate  the  basic  inequalities. 

1.  The  n  variables,  z  ,  which  represent  the  n  variables,  w  ,  in 
a  block  of  inequalities  are  called  type  W  variables.   The  co- 
efficients of  type  W  variables  follow  a  binary  counting  up 
pattern  from  (0,  0,  ...,  O)  to  (1,  1,  ...,  l)  as  j  goes  from 
1  to  2^. 

2.  If  the  complements  of  the  external  variables  are  available  as 
inputs  to  the  network,  then  the  n  variables,  z  ,  which  repre- 
sent  the  n  variables,  v, ,  in  a  block  of  inequalities  are  called 

type  V  variables.  The  coefficients  of  type  V  variables  follow 
the  binary  counting  down  pattern  (1,  1,  ...,  l),  (1,  1,  •••,  1, 
O),  ...,  (0^  0,  ...,  O)  as  j  goes  from  1  to  2  . 

3.  The  r  =  2  variables,  z  ,  which  represent  the  r  variables, 
p-'l   }   or  the  r  variables,  P.   ,  are  called  type  D  variables. 
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Type  D  variables  have  the  following  property:   if  dz  appears 
in  the  first  inequality  of  a  block^  then  dz    appears  in  the 
second  inequality,  etc.,  where  d  is  some  integer  constant. 
Also,  for  networks  which  have  more  than  one  gate  type  (e.g., 
nor/and),  the  2  variables,  z  ,  which  represent  the  2  variables, 
7.,   are  type  D  variables.   The  name  "type  D"  is  suggested  by 

IK 

the  "_diagional"  pattern  followed  by  the  variable  subscripts. 
k.      The  variable,  z  ,  which  represents  the  variable,  oc     ,    in  a 
block  of  inequalities  is  called  a  type  C  variable.   Type  C 
variables  have  the  property  that  dz  appears  in  each  inequality 
of  the  set.   A  non-zero  constant  term  in  a  block  of  inequalities 
is  also  defined  to  be  of  type  C.   In  networks  with  more  than  one 
gate  type,  the  variable,  z  ,  which  represents  the  variable,  k. , 
is  also  a  type  C  variable.   The  name  "type  C"  is  suggested  by 
the  fact  that  the  variable  subscript  remains  _constant  throughout 
the  block  of  inequalities. 

Before  GEN-IEQ,  and  INIT  generate  the  inequalities,  INIT  reads  the 
problem  parameter  card  (see  Chapter  6)  which  contains  values  for  n  and  R. 
Using  these  values  GEN-IEQ  and  INIT  generate  the  basic  inequalities  one  block 
at  a  time  in  the  following  manner.   GEN-IEQ,  provides  the  necessary  information 
about  the  first  inequality  of  a  block  by  specifying  values  to  three  arrays. 
TEMP  is  used  to  store  the  subscripts,  h,  of  the  variables,  z  ;  VAL  is  assigned 
the  values  of  the  coefficients  of  the  variables;  and  TP  is  assigned  the  type 
of  each  variable.   If  the  constant  term,  b,  in  the  first  inequality  of  the 
block  is  non-zero,  this  information  is  stored  in  TEMP(l),  VAL(l)  and  TP(l) 
since  other  subroutines  of  the  program  are  structured  under  the  assumption 
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that  the  constant  term  is  stored  first  in  TEMP,  VAL  and  TP.  The  non-zero 
constant  term  is  represented  by  TEMP(l)  =  0,  VAL(l)  =  b  and  TP(l)  =  C. 
TEMP(l)  is  set  to  0  to  distinguish  the  constant  from  the  other  variables  in 
the  inequality  and  TP(l)  is  specified  as  C  since  the  non-zero  constant  term  is 
defined  to  be  of  type  C.   The  information  about  the  other  variables  in  the  first 
inequality  is  placed  in  the  three  arrays  in  an  arbitrary  order.   If  the  first 
inequality  has  n  type  W  or  V  variables  only  the  variable  with  the  smallest 
subscript  is  specified  by  GEN-IEQ.   The  other  variables  will  be  generated 
automatically  by  INIT. 

When  GEN-IEQ  has  stored  all  of  the  information  about  the  first 
inequality  of  a  block,  the  program  transfers  to  INIT  where  all  the  variables 
except  the  constant  term  are  rearranged  according  to  the  descending  order  of 
the  absolute  value  of  the  variable  coefficients.   This  rearranging  makes  it 
easier  for  CHK-IEQ  to  check  condition  (c)  described  on  page  2?.  INIT  then 
uses  the  information  about  the  first  inequality  of  the  block  provided  by 
GEN-IBQ,  to  form  all  2  inequalities  of  the  block.  The  process  of  forming 
all  2  inequalities  in  INIT  is  based  on  the  types  which  were  assigned,  to 
the  variables  by  GEN-IE)5,.  These  types  define  the  coefficient  and  variable 
subscript  patterns  for  the  entire  block  of  inequalities.  Finally,  INIT 
stores  the  necessary  information  (coefficients,  subscripts,  etc.)  about  each 
inequality  of  the  block  in  the  proper  arrays.  Further  details  about  these 
arrays  are  not  pertinent  to  this  discussion.  The  program  then  transfers  back 
to  GEN-IB3,  which  stores  the  information  about  the  first  inequality  of  the 
next  block  in  the  arrays  TEMP,  VAL  and  TP. 

As  an  example  consider  the  inequalities  in  (^.7).  GM-IIR  stores 
the  information  for  the  first  inequality  in  the  arrays  TEMP,  VAL  and  TP  as 
follows : 
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\ 

.  L 

\ 

1 

2 

3 

h 

5 

TEMP(L) 

0 

1 

15 

23 

37 

val(l) 

-1 

1 

1 

1 

7 

tp(l) 

c 

W 

D 

D 

D 

Notice  that  the  constant  term  is  moved  to  the  left  hand  side  of  the  inequality 
before  the  arrays  are  filled.   INIT  then  rearranges  the  variables  as  follows : 
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L 

\ 

1 

2 

3 

k 

5 

TEMP(L) 

0 

37 

15 

23 

1 

VAL(L) 

-1 

7 

1 

1 

1 

TP(L) 

C 

D 

D 

D 

W 

(^.11) 


Next,  INIT  uses  the  information  in  (4.11)  to  form  all  four  inequalities  of 
(U.7).   The  constant  term,  -1,  is  included  in  a2J.  four  inequalities  since 
TP(l)  =  C.   TP(2)  =  D,  TP(3)  =  D  and  TP(U)  =  D  cause  INIT  to  include  1^^^, 
Iz   and  Iz^o  in  the  first  inequality;  I'^^^Q}   ■^'^-tf.   ^^^  -^^pii  ^^   ^^®  second 
inequality,  etc.  Since  only  the  variables  with  non-zero  coefficients  in  an 
inequality  are  stored,  TP(5)  =  W  causes  INIT  to  include  Iz^  in  the  second 
inequality;  Iz  in  the  third  inequality;  and  both  Iz  and  Iz  in  the  fourth 
inequality.   Thus  the  inequalities  formed  by  INIT  are 


0  <  -1  +  Tz^^  +  Iz^^  +  Iz^^ 

0  <-l  +   7^38  +  lZi6  +  1^21+  ■*"  ^^2 


0  <  -1  +  7z   +  Iz   +  Iz^^  +  Iz^ 

0  <  -1  +  Tz^Q  +  Iz^g  +  IZg^  +  Iz^  +  Iz^ 


ko 


Notice  that  in  each  inequality  the  variables  are  still  arranged  according  to 
the  descending  order  of  the  absolute  value  of  the  variable  coefficients. 

k.2     Generating  the  Additional  Inequalities 

In  general,  additional  inequalities,  like  those  described  in 
Chapter  2,   cannot  be  grouped  into  blocks  (as  defined  in  the  previous  section). 
Thus  the  type  definitions  given  in  section  k.l   are  not  valid  for  additional 
inequalities.   GEN-IEQ  and  INIT  are  used  to  generate  the  additional  inequali- 
ties in  the  same  way  they  are  used  to  generate  the  basic  inequalities  except 
that  the  additional  inequalities  are  generated  one  at  a  time  rather  than  one 
block  at  a  time.   GEN-IEQ  provides  the  information  about  each  additional  in- 
equality by  specifying  the  appropriate  values  to  the  arrays  TEMP  and  VAL  (TP 
is  filled  with  blanks).   As  before,  if  the  constant  term  in  the  inequality  is 
non-zero,  it  must  be  stored  in  TEMP(l)  and  VAL(l).   The  program  then  transfers 
to  INIT  where  all  of  the  variables  (but  not  the  constant  terra)  are  rearranged 
according  to  the  descending  order  of  the  absolute  value  of  the  variable 

coefficients.  However,  upon  discovering  that  the  types  of  the  variables  have 
not  been  specified,  INIT  merely  stores  the  inequality  in  the  proper  arrays. 
The  program  then  transfers  back  to  GEN-IEQ  which  specifies  values  to  TEMP  and 
VAL  for  the  next  additional  inequality. 

4.3  The  Order  of  the  Inequalities 

When  searching  for  a  free  variable  to  fix,  it  is  sufficient  for 
AGMT-VAR  to  consider  only  a  subset  of  the  basic  inequalities  which  character- 
ize the  gates  in  the  network  (e.g.,  inequalities  (U.l)  for  a  network  with 
NOR  gates)  (see  Chapter  5).   To  make  use  of  this  fact,  AGMT-VAR  was  programmed 
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under  the  assumption  that  the  inequalities  are  stored  in  a  certain  order. 
Thus,  GEN-IBQ,  and  IKTT  generate  the  inequalities  in  the  order  described  below. 
Not  only  are  the  basic  inequalities  which  characterize  the  gates  in  the 
network  generated  first,  but  they  are  also  generated  in  a  specific  order 
which  depends  on  the  gate  types  being  used.   The  inequalities  for  gate  1  (the 
output  gate)  come  first  followed  by  the  inequalities  for  gate  2,  etc.   If  gate 
k  is  a  NOR  or  AND  gate,  the  inequalities  for  gate  k  are  generated  in  the 
following  order: 

[2     inequalities  which  are  restrictive  if  P    =0;   J  =  1,  2,  . . . ,  2  } 


k 


{2     inequalities  which  are  restrictive  i 


If  gate  k  is   a  NMD  or  OR  gate,   the  order  is 


if  pf^^  =1;   0  =  1,  2,  ...,  2^'} 


{2  inequalities  which  are  restrictive  if  pA   =1;   J  =  1^  2,  . . . ,  2  } 
[2     inequalities  which  are  restrictive  if  'P-^      =0?   J  =  1^  2,  . . . ,  2  ] 


For  example,  in  a  NGR/^AND  network  synthesis  problem  the  basic  inequalities 
which  characterize  the  gates  are  generated  first  and  in  the  following 
order  (of  course  for  any  gate  k  the  NAND  inequalities  could  come  before  the 
NOR  inequalities): 
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[2  NOR  inequalities  which  are  restrictive  if  P^  =^  0 
(2  NOR  inequalities  which  are  restrictive  if  P^''  =  1 
[2     NMD  inequalities  which  are  restrictive  if  P  "^  "^  =  1 


2  NAND  inequalities  which  are  restrictive 


if  f[^  ^ 


3 

=  1, 

2, 

...,    2^} 

d 

- 1, 

2, 

...,    2^} 

J 

=  1, 

2, 

...,    2^} 

0 

=  1, 

2, 

...,    2") 

n  (  t) 

[2     NOR  inequalities  which  are  restrictive  if  Vlr     =  0 

[2     NOR  inequalities  which  are  restrictive  if  Vr^     -   1 

K 


-,n 


{ 2  NAND  inequalities  which  are  restrictive  if  P. 


(j)^  1 


R 

[2     NAND  inequalities  which  are  restrictive  if  P}>   =  0 

K 


j 

=  1, 

2, 

...,    2"} 

J 

=  1, 

2, 

...,    2^} 

J 

=  1, 

2, 

...,    2^} 

j 

=  1, 

2, 

...,    2"] 

The  remaining  basic  inequalities  and  the  additional  inequalities  are  generated 

in  an  arbitrary  order. 
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CHAPTER  5 
THE  SCHEME  FOR  AGMT-VAR 

The  objective  of  AGMT-VAR  is  to  select  an  inequality  which  is  not 
satisfied  and  then  fix  a  free  variable  in  that  inequality  according  to  the 
scheme  described  in  this  chapter.   AGMT-VAR  does  not  need  to  consider  all 
of  the  inequalities  generated  by  GEN-IEQ,.   All  additional  inequalities  can 
be  ignored  by  AGMT-VAR  since  the  variables  in  these  inequalities  will  be 
assigned  values  by  CHK-IB5,  (specifically,  by  either  condition  (b)  or  (c)  on 
page  27).  Also,  basic  inequalities  such  as  inequalities  (2.2)  for  NOR  gate 
networks  need  not  be  considered  by  AGMT-VAR  since  they  will  be  satisfied  by 
conditions  (b)  and  (c)  of  CHK-IEQ,,   Finally,  basic  inequalities  of  the  form 
(2,1)  which  are  restrictive  when  P^   =  1  for  NOR  and  AND  gates  and  when 
P^^^  =  0  for  NAND  and  OR  gates  will  be  satisfied  by  condition  (b)  of  CHK-IH^ 
and  need  not  be  considered  by  AGMT-VAH.   For  example,  all  of  the  basic  inequal- 
ities generated  by  the  second  expression  in  (2.l)  can  be  ignored  by  AGMT-VAR. 
Thus  AGMT-VAR  only  needs  to  consider  a  subset  of  the  inequalities. 

It  is  important  to  note  the  difference  between  the  subset  of  in- 
equalities which  must  be  considered  by  AGMT-VAR  and  the  inequalities  which 
must  be  checked  by  CHK-IEQ.   AGMT-VAR  must  consider  a  fixed  subset  of  in- 
equalities, while  the  subset  of  inequalities  checked  by  CHK-IEQ,  continually 
changes  during  the  computation.   To  make  the  most  efficient  use  of  these 
facts,  the  basic  inequalities  are  generated  in  a  certain  order  (see  Chapter 
^)  to  simplify  AGMT-VAR,  and  the  method  of  chain  checking  (see  Chapter  3) 
is  used  in  CHK-IEQ. 

For  each  inequality  (call  it  i)  which  it  must  consider,  AGMT-VAR 


takes  the  following  action.  First,  the  value  of  the  P  variable  is  checked. 
If  the  P  variable  is  free,  the  i-th  inequality  is  skipped.   If  the  P  variable 
is  fixed,  AGMT-VAR  checks  ^.(F).   If  ^. (f)  >  0,  the  inequality  is  skipped 
since  it  is  already  satisfied.   If  ^.  (F)  is  negative,  the  X.  variable  is 
checked  if  the  network  has  more  than  one  gate  type.   The  X.  variable  is  set 
to  1  if  it  is  free.   Otherwise,  the  only  free  variables  in  the  inequality  are 
either  w  variables,  v  variables  or  p  variables.   Each  free  w,  v  or  3  variable 
is  assigned  a  type  according  to  the  following  definitions  [l].   The  variable 
type  is: 

k     k 
VC'^  :    if  the  variable  is  w  or  v  (a  Variable  with  Connection  being  *). 

G-^  :    if  the  variable  is  3^^  with  cl^  =   1  and  P^'^  ^  =  *  (a  Gate  h  with 

p(j)  ^  ^). 
h     — 

GC-^  :    if  the  variable  is  a^^  with  P^^^  ^  ^  *   and  0'   =  *  where  gate  h  is 

not  isolated  (a  Gate  h  with  Connection  being  *). 
G^*:  if  the  variable  is  p^^  with  P^'^"  ^  =  *   and  a_  =  *  where  gate  h  is 

not  isolated  (a  Gate  h  with  P/   =  *  and  Connection  being  *). 


NWG 


if  the  variable  is  ^^^^   with  P^"  ^  =  *  and  ol^   =  *  where  gate  h  is 
isolated  (aNeWGate). 

The  above  variable  types  can  be  ordered  from  most  desirable  to 


t 


+ 


During  the  computation  a  variable  can  have  one  of  three  values:   0  (zero), 
1  (one)  or  -^  (a  free  variable). 

A  gate,  h,  is  said  to  be  isolated  if  it  is  not  connected  to  any  other  gates 
in  the  current  partial  solution  (i.e.,  if  all  a   (k^h)  are  0  or  *) . 


h3 
least  desirable  as 

G^,    VC*,    GC-^,  G*C-^  and  NWG  (5-1) 

Why  the  desirability  ordering  in  (5«l)  is  used  is  more  obvious  in  some  cases 
than  others.   For  example,  G-^  is  preferred  to  G-^C*  because  G-^C^  requires  both 
an  additional  interconnection  (cx  v.)  and  the  fixing  of  a  free  variable  (P}      ) 
while  G*  only  requires  that  a  free  variable  (P/^*^  )  be  fixed.   On  the  other 
hand,  it  is  not  intuitively  obvious  why  G-^  is  preferred  to  VC*.   VC*  requires 
an  additional  interconnection  while  G*  does  not.   However,  the  fixing  of  a 
free  variable  required  by  G'^  could  make  some  other  inequality  harder  to  satisfy, 

It  was  empirically  concluded  that  (5.l)  is  the  best  desirability  ordering. 
Finally,  the  type  of  an  inequality  is  defined  to  be  the  same  as  the  best  vari- 
able type  for  that  inequality.   The  best  variable  type  is  chosen  since  it 
tends  to  minimize  the  number  of  interconnections. 

In  AGMT-VAR  the  inequalities  which  must  be  considered  are  scanned 
one  at  a  time  (looking  for  free  w,  v  and  p  variables),  and  each  inequality  is 

given  a  type.   If  at  any  point  the  inequality  being  scanned  is  determined  to 

(1)      t 
be  of  type  NWG,  AGMT-VAR  immediateily  sets  the  type  NWG  variable,  p.^  ,  to  1 

Ik 

and  underlines  it.   The  gate  connected  to  the  network  is  set  to  a  type  dif- 
ferent from  the  gate  it  feeds  (i.e..  A..  ^  K)   if  the  network  has  more  than 
one  gate  type  (e.g.,  NOR/AND).   Also,  if  some  inequality  is  scanned  and  does 
not  have  any  free  variables,  the  current  partial  solution  is  infeasible. 
Otherwise,  the  inequality  with  the  least  desirable  type  is  chosen  and  the 
free  variable  with  the  best  variable  tjrpe  for  that  inequality  is  fixed.   If 


Subsequently,  a   and  P.   are  set  to  1  when  the  inequalities  are  checked 

again  by  CHK-IEQ.   Thus  AGMT-VAR,  in  effect,  connects  the  isolated  gate 
to  the  network. 
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two  inequalities  have  the  least  desirable  type,  the  one  with  the  smallest 
number  of  free  variables  is  chosen.   The  inequality  with  the  least  desirable 
type  is  chosen  for  two  reasons.   First,  if  the  problem  has  a  feasible  solu- 
tion, the  inequality  with  the  least  desirable  type  must  be  satisfied  some- 
time; and  if  it  is  not  satisfied  as  soon  as  possible,  the  number  of  free 
variables  for  that  inequality  may  be  smaller  when  CHK-IEQ  is  called  again 
later  in  the  calculation.   Consequently,  that  inequality  may  be  harder  to 
satisfy.   Second,  if  the  current  partial  solution  may  lead  to  solutions  which 
exceed  the  current  bound  on  the  value  of  the  objective  function,  the  bound 
may  be  exceeded  earlier  since  the  cost  (in  terms  of  the  number  of  gates  and/ 
or  connections)  associated  with  the  inequality  with  the  least  desirable  type 
tends  to  be  greater.   Each  time  AGMT-VAR  fixes  a  free  variable  the  program 
transfers  to  CHK-IEQ,. 


hi 
CHAPTER  6 

INPUT 

The  easiest  and  most  practical  approach  for  generating  the  inequali- 
ties for  ILLODIE-AIF  is  to  use  GEN-IEQ  and  INIT  (see  Chapter  ^4-).   However, 
early  versions  of  the  program  read  the  inequalities  from  cards,  and  this 
capability  is  maintained  in  the  current  version  of  the  program.   Thus,  the 
discussion  of  the  input  data  cards  for  ILLODIE-AIF  given  in  this  chapter 
includes  a  discussion  of  how  to  provide  the  inequalities  on  cards. 

The  following  definitions  facilitate  the  discussion  of  the  input 
data  cards. 

(a)  A  problem  refers  to  using  ILLODIE-AIF  to  synthesize  networks 
which  realize  a  set  of  one  or  more  switching  functions. 

(b)  A  problem- set  refers  to  a  group  of  one  or  more  problems  which 
use  the  same  parameter  card  (to  be  discussed  later  in  this 
chapter)  and  the  same  inequalities. 

(c)  A  multiple-problem- set  refers  to  a  group  of  one  or  more  problem- 
sets.   If  the  inequalities  are  generated  by  the  program,  all 
networks  synthesized  for  a  multiple-problem- set  must  be  of  the 
same  gate  type  (or  types)  since  only  one  version  of  GEN-IEQ, 

can  be  used  at  a  time.   However,  if  the  inequalities  are  read 
from  cards,  the  networks  synthesized  for  each  problem-set  may 
be  of  a  different  gate  type  (or  types )  from  the  networks  for 
the  other  problem-sets. 
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The  order  of  the  cards  in  the  input  deck  for  a  problem  and  a 
multiple -problem- set  are  shown  in  Fig.  6.1  and  Fig.  6.2,  respectively.  Since 
all  arrays  and  nearly  all  variables  in  the  program  are  defined  as  half-word 

integers,  all  numerical  data  with  three  exceptions  must  be  in  the  range 

-15     15 
from  -2    to  2   -1.  All  numbers  are  read  using  a  FORTRAN  I  format  which 

requires  that  they  be  right  justified  in  their  assigned  fields  on  the  data 

cards  to  be  read  correctly. 

1.  THE  TITLE  CARD 

Each  problem-set  must  have  exactly  one  title  card  which  may  contain 
any  suitable  title  for  the  problem- set.  The  contents  of  the  title  card  will 
be  printed  as  a  heading  for  the  output  (see  Chapter  7)-   The  program  does  not 
try  to  obtain  any  information  about  the  problem-set  from  the  title  card. 

2.  THE  PARAMETER  CARD 

This  card  contains  six  problem-set  parameters  each  of  which  occupies 
five  columns. 

(a)  NOVR 

NOVR,  punched  in  columns  1  through  5^  is  the  number,  n, 

of  external  variables  (x  ,  x   ,    . . . ,   x  )  of  the  given  switching 

(output)  function(s).   Theoretically  NOVR  could  be  any  number, 


t 


NOGT,  NOVR  and  NOPT  which  will  be  explained  later  are  defined  as  full-word 

t       -?1     ^1 
integers  and  can  range  from  -2    to  2   -1.   Actually,  all  the  data  in 

network  synthesis  problems  are  well  within  these  ranges. 


These  ranges  are  for  the  IBM  360/75. 
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END  CARD 


INITIAL  PARTIAL  SOLUTION  CARD(S) 


LAST  FUNCTION  TITLE  CARD 


LAST  FUNCTION  CARD 


ZBAR  CARD 


FIRST  FUNCTION  TITI£  CARD 


FIRST  FUNCTION  CARD 


INEQUALITY  CARDS 


OBJECTIVE  FUNCTION  CARD(S) 


PARAMETER  CARD 


TITLE  CARD 


Figure  6.1.      Tlie   Input  Deck  for  One  Problem 
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END  CARD  (LAST  PROBLEM-SET) 


ZBAR  CARD,  FUNCTION  AND  FUNCTION  TITLE 
'  CARDS,  INITIAL  PARTIAL  SOLUTION  CARD(S) 
.  (LAST  PROBLEM,  LAST  PROBLEM-SET) 


ZBAR  CARD,  FUNCTION  AND  FUNCTION  TITLE 
CARDS,  INITIAL  PARTIAL  SOLUTION  CARD(S) 
FIRST  PROBLEM,  LAST  PROBLEM-SET) 


TITLE  CARD,  PARAMETER  CARD,  OBJECTIVE 
FUNCTION  CARD(S),  INEQUALITY  CARDS 

(LAST  PROBLEM- SET) 


END  CARD  (FIRST  PROBLEM-SET) 


ZBAR  CARD,  FUNCTION  AND  FUNCTION  TITLE 
CARDS,  INITIAL  PARTIAL  SOLUTION  CARD(S) 
(LAST  PROBLEM,  FIRST  PROBLEM-SET) 


ZBAR  CARD,  FUNCTION  AND  FUNCTION  TITLE 
CARDS,  INITIAL  PARTIAL  SOLUTION  CARD(S) 
(FIRST  PROBLEM,  FIRST  PROBLEM-SET) 


TITLE  CARD,  PARAMETER  CARD,  OBJECTIVE 
FUNCTION  CARD(S),  INEQUALITY  CARDS 
(FIRST  PROBLEM-SET) 


Figure  6.2.   The  Input  Deck  for  a  Multiple-Problem-Set. 
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but  practical  considerations  (the  amount  of  core  storage  needed, 
the  number  of  inequalities,  the  computation  time,  etc.)  limit 
NOVR.   In  ILLODIE-ArF  NOVR  must  be  less  than  or  equal  to  four 
although  the  program  could  be  modified  to  handle  functions  of 
more  than  four  variables, 
(b)  NOGT 

NOGT  is  the  number  of  gates  in  the  network  and  is  punched 
in  columns  6  through  10.  The  output  subroutines  are  designed, 
to  handle  networks  with  a  maximum  of  ten  gates.  However, 
these  subroutines  could  be  changed  to  allow  for  networks 
with  more  than  ten  gates. 

Once  set,  NOGT  remains  fixed  throughout  the  calculation  and, 
if  the  problem  is  feasible,  all  networks  with  a  minimum  number  of 

interconnections  for  that  value  of  NOGT  are  found. 

What  value  of  NOGT  to  choose  and  how  to  choose  it  depend 
upon  the  objective  of  the  user.   The  usual  objective  for 
using  ILLODIE-AIF  is  to  find  the  network  or  networks  for  a 
given  switching  function  which  have  first  the  minimum  number 
of  gates  and  second  the  minimum  number  of  interconnections.   The 
following  approach  is  reasonable  for  finding  the  minimum  value 
of  NOGT.   Choose  a  reasonable  estimate  of  the  optimal  value  of 
NOGT  (perhaps  by  a  simple  hand  calculation  using  a  Karnaugh  map) 
and  attempt  to  solve  the  problem  using  ILLODIE-AIF.   If  the 
problem  is  infeasible  (i.e.,  no  solution  exists)  increase  the 
value  of  NOGT  by  1  and  run  the  problem  again.   Repeat  this 
process  until  the  first  value  of  NOGT  is  found  for  which  a 
solution  exists. 
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This  value  of  NOGT  is  the  minimum  value.   On  the  other  hand  if 
the  initial  estimate  of  NOVR  results  in  a  feasible  solution, 
decrease  the  value  of  NOVR  by  1  until  the  first  value  of  NOGT 
is  found  for  which  the  problem  is  infeasible.  The  value  of 
NOGT  1  greater  than  this  last  value  is  the  minimum  value 
of  NOGT. 

(c)  NOKGT 

NOKGT  is  the  number  of  gate  types  in  the  network  and  is 
punched  in  columns  11  through  15.   In  ILLODIE-AIF  NOKGT  is 
either  1  or  2.   The  problem  would  have  to  be  reformulated  if 
the  user  wanted  to  synthesize  networks  with  more  than  two 
gate  types  (for  example,  k   variables  which  can  only  assume  the 
values  0  and  1  cannot  be  used  for  networks  with  more  than  two 
gate  types). 

(d)  TYIN 

TYIN,  punched  in  columns  l6  through  20,  indicates  whether 
the  inequalities  will  be  read  from  cards  or  generated  by  the 
program.   TYIN  =  PG  if  the  program  will  be  used  to  generate  the 
inequalities,  and  TYIN  =  CD  if  the  inequalities  will  be  read 
from  cards. 

(e)  CPNCP 

CPNCP  indicates  whether  or  not  the  complements  of  the 
external  variables  (as  a  group)  are  available  and  is  punched 
in  columns  21  through  25.   The  complements  are  available  if 
CPNCP  =  CP  and  are  not  available  if  CPNCP  =  NC. 
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(f)  NOPT 

NOPT  Indicates  the  number  of  output  gates  in  the  network 
(i.e.,  the  number  of  output  functions  the  network  is  to  realize) 
and  is  punched  in  columns  26  through  30,  The  current  version 
of  the  program  is  capable  of  handling  networks  with  any  number 
of  output  functions  subject  to  the  restrictions  on  the  sizes 
of  NOGT,  WOVE,  etc. 

3.   THE  OBJECTIVE  FIMCTION  CARD(S) 

If  TYIN  =  PG,  the  objective  function,  cz,  will  be  generated  by  the 
program.  However,  if  the  inequalities  are  read  from  cards,  the  user  must 
also  specify  an  objective  function  on  cards.   The  objective  function  is 
punched  in  col"umns  1  through  72  of  each  card.   Columns  73  through  80  are  not 
read  by  the  program.   The  first  seventy-two  colimins  are  divided  into  eighteen 
fields  with  four  columns  to  a  field.   The  numbers  are  read  from  the  cards  in 
pairs.  The  first  nimiber  read  is  the  coefficient  of  the  variable  and  the  second 
number  is  the  subscript.   Thus  columns  1-U,  9-12,  ,,,,  and  65-68  contain 
variable  coefficients  and  columns  5-8,  13-l6,  ,,.,  and  69-72  contain  the 
respective  subscripts.   For  example,  if  columns  1  through  8  of  an  objective 
function  card  contain  bbb^+bblO,  then  ^z^^   is  part  of  the  objective  function. 
Only  the  non-zero  coefficients  and  their  respective  subscripts  need  to  be 
specified.  Also,  the  numbers  are  usually  punched  on  the  cards  according  to  the 
ascending  order  of  the  subscripts.   The  end  of  the  objective  function  is  signaled 
by  putting  a  negative  number  in  the  subscript  field  immediately  following  the 
subscript  field  of  the  last  non-zero  coefficient.   If  the  last  non-zero 


The  letter  b  is  used  to  represent  a  blank, 
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coefficient  occurs  in  columns  65  through  68,  then  an  additional  card  with  a 
negative  number  in  columns  5  through  8  is  necessary. 

k.      THE  INEQUALITY  CARDS 

As  in  the  case  of  the  objective  function,  the  inequalities  are 
punched  in  columns  1  through  72  only.   Each  card  is  divided  into  nine  fields 
with  eight  columns  in  each  field.   Each  field  is  used  as  follows:   the  first 
column  contains  the  type  of  the  variable  (here,  the  variable  types  are  those 
defined  in  Chapter  h   rather  than  Chapter  5  and  are  denoted  as  b  (blank),  V,  W, 
D  or  C);  the  next  three  columns  contain  the  value  of  the  variable  coefficent; 
and  the  last  four  columns  contain  the  subscript  of  the  corresponding  variable. 
The  numbers  on  the  inequality  cards  are  read  into  the  array  INPT  which  is 
dimensioned  such  that  each  inequality  is  limited  to  a  maximum  of  ten  cards 
(i.e.,  90  non-zero  coefficients).   If  an  inequality  has  a  non-zero  constant, 
it  must  be  punched  in  the  first  field  of  the  first  card  for  that  inequality. 
The  remaining  fields  must  be  filled  according  to  the  descending  order  of  the 
absolute  value  of  the  variable  coefficients.   The  end  of  an  inequality  is 
indicated  by  punching  a  negative  number  in  the  last  four  columns  of  the  field 
immediately  following  the  field  which  contains  the  last  non-zero  coefficient. 
If  the  last  non-zero  coefficient  occurs  in  the  last  field  of  a  card,  an 
additional  card  with  a  negative  number  in  columns  5  through  8  is  required. 
The  end  of  the  inequalities  is  signaled  by  a  special  card  which  contains 
any  symbol  (except  b,  V,  W,  D  or  C)  in  the  first  column  and  a  negative  number 
in  columns  13  through  I6. 

For  the  basic  inequalities  the  user  specifies  the  information  about 
the  first  inequality  of  each  block  and  INIT  generates  the  entire  block  of 
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inequalities  as  described  in  Chapter  h   (also,  see  Appendix  A  for  an  example). 
Furthermore,  the  order  of  the  inequality  cards  must  be  the  same  as  the  order 
described  in  Section  ^.3.   If  the  program  is  used  to  generate  the  inequali- 
ties, the  objective  function  and  inequality  cards  are  simply  omitted. 

5.  ZBAR  CARD 

ZBAR,  p-unched.  in  columns  1  through  5,  is  an  upper  bo\and  on  the  value 
of  the  objective  function.  If  no  upper  bound  is  known,  ZBAR  =  NOVR  X  NOGT  + 
NOGT  X  (NOGT-l)  is  always  sufficient  since  NOYR  X  NOGT  allows  each  external 
variable  to  be  connected  to  every  gate  and  NOGT  X  (NOGT-l)  allows  each  gate 
to  be  connected  to  every  other  gateo 

6.  THE  FUNCTION  CARD 

The  function  card  contains  the  truth  table  of  a  switching  (output) 
function  and  must  be  punched  (one  column  for  each  value)  beginning  in 
column  6.   Since  the  program  can  solve  problems  with  up  through  four  vari- 
ables, the  number  of  columns  needed  will  vary  (sixteen  columns  for  a  four 
variable  function,  eight  columns  for  a  three  variable  function,  etc.)«   The 
truth  table  values  must  be  punched  in  the  following  order  (in  this  case  for 
a  three  variable  function):   f(0,0,0),  f(0,0,l),  ...,  f(l,l,0),  f(l,l,l). 
For  a  problem  with  more  than  one  output  function,  a  function  card  is  punched 
for  each  function  and  placed  in  the  input  deck  as  shown  in  Fig.  6.1. 

7.  THE  FUNCTION  TITLE  CARD 

Following  each  function  card  there  must  be  exactly  one  function 
title  card  which  may  contain  any  suitable  information  about  the  output 
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function  such  as  its  algebraic  representation.  The  function  title  card  is 
read  and  printed  as  explained  in  Chapter  7«   It  is  convenient  but  not  neces- 
sary to  punch  the  truth  table  of  the  output  function  on  the  function  title 
card  in  exactly  the  same  columns  as  on  the  function  card.   With  this  pre- 
caution the  function  and  function  title  cards  could  accidently  be  inter- 
changed and  the  program  would  still  work. 

8.   THE  INITML  PARTIAL  SOLUTION  CARD(S) 

If  a  partial  solution  to  the  problem  is  known,  the  user  can  specify 
it  on  the  initial  partial  solution  cards.   This  is  done  by  spejif^ing  values 
which  are  read  into  the  arrays  X  and  XS  which  are  defined  on  page  32.   Only 
columns  1  through  72  of  each  card  are  read,  and  each  card  is  divided  into 
twelve  fields  with  six  columns  in  each  field.   The  first  five  columns  of  each 
field  contain  the  subscript  of  a  variable  in  the  initial  partial  solution. 
This  number  is  read  into  the  array  XS  and  must  be  negative  if  the  variable 
is  to  be  underlined.   The  last  column  of  each  field  is  used  to  determine 
whether  the  variable  is  to  be  set  to  zero  or  one.   If  the  number  in  the  last 
column  of  the  I-th  field  is  1,  X(l)  is  set  to  XS(l).   If  the  number  is  0, 
X(l)  is  set  to  -XS(l).  For  example,  if  the  first  variable  in  the  initial 
partial  solution  is  to  be  Xo  and  it  should  be  set  to  1  and  \ander lined,  then 
the  first  six  columns  of  the  first  initial  partial  solution  card  should 
contain  bbb-80.   If  the  next  variable  is  x  ,  and  it  should  be  set  to 
1  and  not  underlined,  then  columns  7  through  12  should  contain  bbblUl. 
Such  a  card  would  set  XS(l)  to  -8,  X(l)  to  8,  XS(2)  to  Ik,   and  X(2)  to  1^. 
Thus  columns  1  -  5,  7  -  11,  ...,  and  67  -  7I  contain  subscripts  (either 
positive  or  negative)  while  columns  6,  12,  ...,  and  72  indicate  to  what 
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value  the  variables  are  to  be  set.   All  columns  to  the  right  of  the  field 
which  corresponds  to  the  last  variable  in  the  initial  partial  solution  must 
be  blank.   If  the  last  variable  is  in  the  last  field  of  a  card,  a  blank 
card  must  be  inserted  next  in  the  input  deck  to  signal  the  end  of  the  initial 
partial  solution.   If  the  user  does  not  wish  to  specify  an  initial  partial 
solution,  a  blank  card  must  be  inserted  in  the  input  deck  where  the  initial 
partial  solution  cards  would  have  been. 

9.   THE  END  CARD 

The  end  of  each  problem- set  is  signaled  by  an  EKD  card  which 
contains  a  negative  number  in  the  first  five  columns.   When  ILLODIE-AIF 
finishes  a  problem,  it  begins  the  next  problem  by  reading  the  ZBAR  card,  the 
function  card(s),  the  function  title  card(s)  and  the  initial  partial  solution 
card(s)  for  the  next  problem.   When  the  program  detects  the  end  of  a  problem- 
set  (i.e.,  reads  and  END  card),  it  begins  the  next  problem-set  by  reading 
the  title  card,  parameter  card,  objective  function  card(s),  inequality 
card(s),  etc.,  for  the  next  problem-set.  The  program  terminates  when  there 
are  no  more  input  cards  to  be  read. 
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CHAPTER  7 

OUTPUT 

In  addition  to  printing  all  feasible  solutions  found  (assuming  the 
problem  is  feasible),  ILLODIE-AIF  also  prints  some  job  statistics  and  other 
information  which  may  be  helpful  in  determining  whether  or  not  the  program 
is  working  properly.   The  purpose  of  this  chapter  is  to  describe  both  the 
printed  and  punched  output  of  the  program.   A  complete  example  of  the  printed 
output  is  given  in  Appendix  A. 

First,  the  title  card  (see  Chapter  6)  is  printed  as  a  heading  for 
the  output.   This  is  followed  by  a  title  for  the  variable  table  (Table  2.i+.l) 
and  then  the  variable  table  itself.   Of  course  w.  is  represented  as  W(l,  l), 
V.  as  V(l,  l),  etc.,  and  each  variable  z  is  represented  by  its  subscript  t. 

Next  the  inequalities  for  the  problem  ai^e  printed  twice.   First, 
they  are  printed  in  exactly  the  same  format  they  would  be  in  if  they  were 
provided  on  cards,  i.e.,  they  are  printed  as  they  appear  in  the  arrays  TP, 
VAL  and  TEMP  after  the  variables  have  been  rearranged  according  to  the  descend- 
ing order  of  the  absolute  value  of  the  coefficients  but  before  INIT  forms 
the  blocks  of  basic  inequalities  (see  Chapter  k) .      The  inequalities  are  also 
printed  after  INIT  forms  the  blocks  of  basic  inequalities.   This  time  they 
are  printed  in  the  form  -b  <  A  z  where  the  symbol  "<"  is  replaced  by  the 
letters  "LE".   In  both  cases  only  the  variables  with  non-zero  coefficients 
are  printed.   In  between  the  two  printings  of  the  inequalities,  the  total 
number  of  variables  N,  the  number  of  inequalities  m,  the  number  of  external 


If  the  inequalities  are  read  from  cards  (TYIN  =  CD),  the  first  printing 
of  the  inequalities  is  omitted. 
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varialDles  n,  the  number  of  non-zero  coefficients  in  the  inequalities,  and 
the  objective  function  are  printed. 

To  separate  the  printing  of  the  inequalities  from  the  printing  of 
the  solutions,  the  function  title  is  printed  ten  times  at  the  bottom  of  two 
consecutive  pages.   If  the  problem  has  more  than  one  output  function,  each 
function  title  is  printed  ten  times  with  the  first  function  title  printed 
last. 

Next  the  value  of  each  partial  sum  y.  (F^^)  (i  =1,  2,  ...,  m)  (see 
Chapter  3)  is  printed  where  F_  is  the  initial  partial  solution.   F„  is  then 
printed  in  the  form  of  the  arrays  X  and  XS  followed  by  the  time  required 
for  initialization  (i.e.,  the  time  required  by  INIT  and  C^EN-IEQ).   The  rest 
of  the  output  depends  upon  the  outcome  of  the  problem. 

ILLODIE-AIF  prints  the  time,  the  current  partial  solution  (X  and 
XS)  and  a  line  of  statistics  every  1000  iterations  where  each  entry  into  the 
subroutine  CHK-IEQ,  is  counted  as  one  iteration.   The  time  value  printed  is 
the  amount  of  execution  time  elapsed  since  the  initialization  time  value  was 
printed.   The  line  of  statistics  (hereafter  referred  to  as  LINESTAT)  includes 
OF  (the  value  of  the  objective  function),  NO.  ITER  (the  number  of  iterations), 
EKTDX  UNFS  (the  number  of  times  the  program  transfers  from  CHK-IEQ  to  BCKTRK 
because  the  current  partial  solution  is  infeasible),  EXTDX  FSBL  (the  number 
of  times  CHK-IEQ,  finds  a  feasible  solution),  CGPSNS  (the  number  of  times 
CHNGMM  is  called  which  is  equal  to  the  total  number  of  variables  fixed  by 


University  of  Illinois  timing  subroutines  are  used  to  provide  the  time 
values  printed  by  the  program. 

Values  called  GETYS,  CYCLES  and  OFBND  are  also  printed  even  though  they 
are  no  longer  meaningful. 
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the  program)  and  FIXJ  (the  number  of  times  the  subroutine  AGMT-VAR  is  entered). 
All  of  these  statistics  are  cumulative  for  each  computer  run  except  OF. 
Also,  the  program  punches  the  current  partial  solution  every  2000  iterations 
in  exactly  the  same  format  as  the  initial  partial  solution  cards  described 
in  Chapter  6.   This  allows  these  cards  to  be  used  as  initial  partial  solution 
cards  ^f  the  program  is  interrupted  before  the  problem  is  finished  and  the 
user  wants  to  resume  the  computation.   After  each  partial  solution  is  punched, 
the  program  also  punches  the  corresponding  value  of  the  objective  function. 
This  card  serves  two  purposes.   First,  it  makes  it  easier  to  separate  the 
different  partial  solutions  which  may  be  punched.   Second,  it  can  be  used 
as  a  ZBAR  card  if  the  problem  is  to  be  continued.   In  this  case  it  must  be 
remembered  to  move  the  ZBAR  card  in  front  of  the  partial  solution  cards  to 
maintain  the  order  of  the  input  cards  described  in  Chapter  6.   The  order  of 
the  cards  punched  by  the  program  is  shown  in  Fig.  T'l-   If  a  problem  is  not 
finished  after  40,000  iterations,  the  program  goes  on  to  the  next  problem. 
The  frequency  (in  terms  of  the  number  of  iterations)  with  which  partial 
solutions  are  printed  and./ or  punched  can  be  changed,  by  changing  the 
appropriate  FORTRAN  DO  loop  parameters  within  the  program. 

If  the  problem  is  feasible,  ILLODIE-AIF  will  print  all  feasible  solutions 
which  give  c  z  a  value  less  than  or  equal  to  the  current  value  of  OF.   The 
value  of  OF  is  updated  each  time  a  feasible  solution  which  lowers  the  value  of 
c  z  is  found.  For  each  such  feasible  solution  the  program  prints  the  values  of 
the  arrays  L,  W,  A  and  V  (if  the  complements  of  the  external  variables  are 
available).  For  example,  let  the  problem  to  be  solved  be  specified  as  a  two 
variable  (NOVR  =2),  two  gate  (NOGT  =2),  NOR  problem  with  no  complements  of 
the  external  variables  available  (CPNCP  =  NC).  Also,  let  the  output  function 
to  be  realized  be  x  sy   x  (Olll).   Then,  one  feasible  solution  as  it  would  be 
printed  by  the  program  is  given  as  follows : 
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ZBAE   CARD     AT  4000  ITEEIATIONS 


PARTIAL  SOLUTION  AT  4000  ITERATIONS 


ZBAR  CARD  AT  2000  ITERATIONS 


PARTIAL  SOLUTION  AT  2000  ITERATIOl^ 


Figure  7.1.   The  Order  of  Partial  Solutions  Punched  by  the  Program. 
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K=12  J=12 
L(l)  0  W(1,K)  0  0  A(1,J)  0  0 
L(2)   0    W(2,K)     1   1    A(2,J)     1  0 

Since  the  problem  has  only  one  gate  type,  the  values  for  the  array  L  can  be 
ignored.   The  values  of  the  W  array  show  that  both  x^  and  x  are  connected 
to  gate  2  while  the  values  of  the  A  array  indicate  that  gate  2  is  connected 
to  gate  1  (the  output  gate). 

LINESTAT  is  printed  immediately  before  each  feasible  solution  is 
printed.   In  addition,  after  the  last  feasible  solution  is  printed,  the 
value  of  each  variable  for  the  last  feasible  solution,  the  time  and  LINESTAT 
are  printed.   This  time  value  is  also  the  amount  of  execution  time  elapsed 
since  the  initialization  time  value  was  printed.   The  final  printing  of 
LIKESTAT  is  useful  for  estimating  how  much  of  the  computation  time  was  spent 
in  proving  the  optimality  of  the  feasible  solutions  found. 

If  the  problem  is  infeasible,  then  the  value  of  each  variable  is 
printed  as  0.   LINESTAT  and  the  amount  of  time  spent  trying  to  find  a  feasible 
solution  are  also  printed. 


63 

APPENDIX  A 

AN  EXAMPLE 

To  provide  an  example  of  the  input  to  and  the  output  from  the  pro- 
gram, let  us  obtain  all  optimal  AND/OR  networks  for  the  switching  function 
A  B  N/  A  B  using  ILLODIE-AIF.   This  example  has  the  following  parameters: 
NOVR  =  n  =  2  (two  external  variables),  NOKGT  =  2  (two  types  of  gates  are 
available),  CPNCP  =  CP  (the  complements  of  the  external  variables  are  avail- 
able), and  NOPT  =  m  =  1  (one  output  function).   Let  us  assume  that  the  net- 
works have  three  gates  (i.e.,  NOGT  =  n  =  3)*   In  addition,  the  objective 
function  and  the  inequalities  are  to  be  read  from  cards  (TYIN  =  CD).   The 
arbitrary  constraint  that  the  output  gate  should  be  an  OR  gate  is  added  to 
illustrate  how  to  provide  an  initial  partial  solution. 

Using  the  variables  w,  ,  v,  ,  a.^  ,   P.?  ,  P^   ,  \    and  7.,   which 

t'  t   ik'   ik  '  k  '  k      xk 

were  defined  in  Chapter  2,  the  basic  inequalities  for  this  example  are  given 
by: 

2  2^ 

Z  vj  xp-^  +  E  ^Jt   x^^^  +  Z  7??^  >  1  -  U  p/^^  -  U  \ 
t=l  ^     ^  t=l  "^  ^     i=l  ^^  -        ^       ^ 

2   V  f.\  2   ,,   ,  .     3 


-  Z  v^  xp^  -  Z  wj;  xp^  -  Z  7??^  >  -U(l-p/^"^)  -  U  A 

i=l    '  - 
iA 


,  .   t  ^t     ,  %  ^t  ^t     .%  ^ik  ^  -"^"-"k   ^  -  ^  ^^ 
t=l  t=l  1=1 


(Al) 


2  2'^ 

Z     w^  x'J)  +  E  v^  x^3)  +  Z  p(j)  >  1  -  U(l-p(^>)  -  U(l-X^) 

t=±  t=l  X=l 

,l^\\  tfl  ^  ^     iSl'i^  >-UP^^)-U(l-^) 

i^ 

k  =  1,  2,  3 

j  =  1,  2,  3,  h. 


6k 


When  A.  =  1^  the  first  two  inequalities  of  (Al)  become  non- restrictive  and  the 
last  two  inequalities  represent  an  OR  gate.   When  >v  =  0,  the  last  two 
inequalities  become  non-restrictive  and  the  first  two  inequalities  represent 
an  AND  gate. 

The  two  non-linear  equations  3.^^  =  a   P;'^  ^  and  7.'^^  =  o:  (l-P.^'^)  = 

lie       XK   1  XK       XK      X 

a   -  P.|!  ai'e  represented  by  the  following  linear  inequalities: 

^)P  -   2a.,  +  ?r^   +  27:,  >  0 
xk      xk    X      xk  — 

-3^[l^  +  3a^^   -  Pp^  -  k7\_l^  >  -1  (A2) 

-  z  s:^-^  +  u  P.  ^  >  0 
t=i  "^     ^  - 
t^i 

i  =  1,  2,  3 

k  =  1,  2,  3 

i  ^^ 

j  =  1,  2,  3,  k. 


The  following  additional  inequalities   are  also   included; 
1.      Each  gate  has  at  least  two  inputs. 
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Z     w^  +     Z     v^  +     E     a       >  2  k  ^  1,    2,    3.  (A3) 

Ul     ^        i=l     ^        i=l     ^^ 

i;^k 


2.   Each  gate,  except  the  output  gate,  has  at  least  one  output 
connection  to  another  gate. 

3 

Z  O^  >  1     k  =  2,  3.  (aU) 

i=l  '^^ 
i^k 
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3.   The  triangular  condition  (see  page  15)  holds  even  if  gate  j  con- 
nects to  other  gates  in  addition  to  gate  k. 

-^ij  -^ik-^Jk>-2  (A5) 

i  -  1,  2,  3 

J  =  1,  2,  3 

k  =  1,  2,  3 

i  ^  J  ^  k. 

k.      If  two  gates,  i  and  k,  are  of  the  same  type  (i.e.,  both  are 

AND  gates  or  OR  gates)  and  gate  i  is  connected  to  gate  k,  then 
gate  i  must  have  at  least  one  more  output  connection,  assuming 
that  no  fan- in  restriction  is  placed  on  the  network. 


3 

"it  '  '"i  '  '"k   ^ik  — 


Z  q;_.^  +  A..  +  >v  -  a_. ,_  >  0 


t=l 

'^  (A6) 


3 

it   "i   '\       ""ik  — 


Z    a,^  -  K   -  K    -  a,,    >  -2 


t=l 

t^i 
t;^k 


i  =  1,  2,  3 
k  =  1,  2,  3 

i  -^  k. 
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5.   Each  optimal  network  must  have  at  least  one  AND  gate  and  one 
OR  gate. 


3 


L     \.    >  1 
i=l  ^- 


-  Z  \.  >  -2 
i=l  "- 


(AT) 


The  input  cards  for  this  example  problem  (shown  in  Table  Al)  are  to 
be  prepared  in  accordance  with  the  discussion  in  Chapter  6.   The  discussion 
on  page  l8  and  Table  2.4.1  are  used  to  provide  the  information  on  the  objec- 
tive function  cards.   The  inequality  cards  are  punched  using  the  description 
of  the  inequalities  given  above,  the  discussion  in  Chapter  ^4,  and  Table  2.U.I. 
The  initial  partial  solution  requires  the  output  gate  to  be  an  OR  gate  as  is 
desired. 

The  output  for  this  problem  follows. 


H  CVJ 

cn 

o  o 

§ 

ON  CVl 
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APPENDIX  B 


SUMMARY  OF  PARAMETERS  AND  VARIABLES 


The  parameters  and  variables  used  in  the  algebraic  descriptions 
of  the  inequalities  and  additional  inequalities  given  in  Chapter  2,   Chapter  h, 
Appendix  A  and  Appendix  C  are  listed  below. 


1.  n  is  the  number  of  external  variables  (x. )  in  the  network. 

2.  R  is  the  total  number  of  gates  in  the  network. 

3.  m  is  the  number  of  output  functions  (gates)  in  the  network. 

k.     w.  represents  the  connection  of  external  variable  x  to  gate  k. 
X)  t 

5.  When  the  complements  of  the  external  variables  are  available, 

k  — 

V,  represents  the  connection  from  external  variable  x  to  gate  k. 

6.  a.,  represents  the  connection  from  gate  i  to  gate  k, 

7.  P.   denotes  the  output  value  of  the  i-th  gate  when  the  external 
variables  connected  to  the  network  assume  the  values  of  input 


vector  X 


(j) 


8.  ^[p   =  a.,    v[^\ 

ik     ik  1 

9.  7.,   =  OL.,   (1-P.  0  (used  only  for  mixed  gate  networks). 

ik     ik    1 

10.   \.  indicates  the  type  of  gate  i  in  mixed  gate  networks.   The 
A.  definitions  are  summarized  in  Table  Bl. 


Gate  types  in  the  network 

A..  =  0          \.  =  1 
1              1 

nor/nand 

nor/and 

and/or 

NAND            NOR 
AND            NOR 

AND             OR 

Table  Bl.   The  \  representation  in  mixed  gate  networks. 
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APPENDIX  C 

ADDITIONAL  INEQUALITIES 

Section  2.3  contains  a  discussion  of  the  additional  inequalities 
for  NOR  networks  in  which  the  complements  of  the  external  variables  are  not 
available.   Also,  the  additional  inequalities  for  AKTO/OR  networks  which  have 
2  external  variables  and  3  gates  are  discussed  in  APPENDIX  A.   The  purpose 
of  this  appendix  is  to  give  a  complete  general  description  of  the  additional 
inequalities  for  NOR,  NOR/NAND,  NOR/AND  and  AND/OR  networks  as  they  are 
generated  by  the  current  versions  of  GEN-IEQ,.   (Also,  differences  between  the 
current  versions  and  earlier  versions  of  the  additional  inequalities  will  be 
mentioned  whenever  appropriate. )  Other  discussions  of  these  additional  in- 
equalities are  given  in  [l],  [5] ,  [12] ,  [1^+] . 

For  each  additional  inequality  given  below,  a  verbal  description 
is  followed  by  an  algebraic  description  in  terms  of  some  of  the  parameters 
and  variables  defined  in  Appendix  B.   Some  of  these  additional  inequalities 
may  have  to  be  modified  or  deleted  if  other  inequalities  such  as  those  for 
fan- in  and/or  fan-out  restrictions  are  added  to  the  program. 

I.   Input  Inequalities 

A.   NOR  Networks 

When  the  complements  of  the  external  variables  are  not  avail- 
able, each  gate  except  the  output  gates  has  at  least  one  input  from 
the  external  variables  or  at  least  two  inputs  from  the  other  gates. 
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n       R 
2  Z  w^  +  Z  a.,  >  2 
t=l  *   i=l  ^^- 

k  =  m+1,  . . . ,  R 

When  the  complements  of  the  external  variables  are  available, 
each  gate  except  the  output  gates  has  at  least  two  inputs. 

n   ^    n       R 
Z    w  +  z  V .  +  z  a.,  >  2 
t.l  *   t=l  ^   1=1  ^^- 

k  =  m+1,  . . , ,    R 
B.   nor/and  Networks 

Each  AND  gate  has  at  least  two  inputs  and  each  NOR  gate  has 
at  least  one  input. 


n        R 

Z  w ,  -f  Z  a.,  +  \.  >  2 
t=l  ^    i=l  ^^    ^- 


or 


k  =  1,  ..._,  R 


"^   k    ^   k    ^ 

Z  w ,  +  Z  V ,  +  z  a.,  +  \.    >  2 

t=l  ^    t=l  ^    iA  ^^     ^" 


k  =  1,  •  •  •  ,  R 


when  the  complements  are  available. 


C.   NOR/NAND  Networks 


Each  gate  has  at  least  one  input. 


n 


k 


R 


z  w  +  z  a   >  1 
t=l     1=1 


iA 


k  =  1,  . .. ,  R 


or 


n   ,     n   ,     R 
Z  w  +  Z  V  +  Z  a   >  1 
t=l  ^   t=l  "^   i=l  "-^  ~ 


i^ 


when  the  complements  are  available. 


k  =  1,  ...,  R 


D.   AKD/OR  Networks 


Each  gate  has  at  least  two  inputs, 


^   k    '^   k    ^ 

Z  w  +  Z  V  +  Z  OL       >2 

t=l  ^   t=l  ^   i=l  ^^ 


iA 


k  =  1,  ...,  R 


II.   General  Triangular  Condition  -  NOR,  NOR/NAMD  and  NOR/AND  Networks, 
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no  other  output  connections 


Figure  CI.   General  Triangular  Condition 
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A.   If  three  gates  are  connected  as  shown  in  Fig.  CI,  then  it  can 
be  proven  that  the  network  is  not  optimal. 

R 

-a.  .  -  a.,  -  a.,  +  L    a.  „  >  -2 
ij    ik    jk   ^^^  j^  - 

i/i 


i 

=  h,    ..., 

R 

J 

=  h,    ..., 

R 

k 

=  1,    • • • , 

R 

i 

;^d  /k 

where  h  =  1  if  m  >  1  and  h  =  2  if  m  =  1.   That  is,  gate  i  or 
gate  j  can  be  gate  1  (always  an  output  gate)  if  the  network  is 
a  multiple -output  network.   (It  should  be  noted  that  in  earlier 
versions  of  the  program  h  =  1  was  always  used  for  NOR  networks 
while  h  =  2  was  always  used  for  NOR/NAITO  and  NOR/AND  networks.) 
B.   Even  if  gate  i  in  Fig.  CI  is  replaced  by  an  external  variable, 
the  network  is  not  optimal.   When  the  complements  are  not 
available: 


k         ^ 
-^.    -  w^  -  Q.,  +  Z  a.  ,,  >  -2 
t    t    jk   ^^^  ji  - 


j 

= 

h, 

. . . , 

R 

k 

= 

1, 

. . . , 

R 

J 

^ 

k 

t 

= 

1, 

0      .      .      , 

n 
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Also,  when  the  complements  x  are  available: 


-vi  -  v_^  -  a.,  +  E  a.n>  -2. 


jk 


i=l 

m 


ji 


j 

= 

h, 

'  •  * } 

R 

k 

= 

1, 

. . . , 

R 

J 

^ 

k 

t 

= 

1, 

...  J 

n 

Again  h  =  1  if  m  >  1  and  h  =  2  if  m  =  1.   (Also,  as  above,  in 
earlier  versions  of  the  program,  h  ^  1  was  always  used  for  NOR 
networks  and  h  =  2  was  always  used  for  NOR/NAND  and  NOR/AND 
networks.,   In  addition,  the  inequalities  for  the  case  when  the 
complements  are  available  were  not  included  previously  for  NOR/ 
MND  and  NOR/AMD  networks  even  when  the  complements  were  avail- 
able. ) 


III.  Extended  Triangular  Inequalities 


other  output  connections  permitted 


Figure  C2.   Extended  Triangular  Condition 
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A.   nor/and  Networks 

If  three  gates  are  connected  as  in  Fig.  C2  and  at  least  one 
of  gates  j  and  k  is  an  AND  gate,  then  the  general  triangular  con- 
dition holds  even  if  gate  j  has  other  output  connections. 


-a.  .  -  a.,  -  a.,  +  \.  >  -2  ">  i  =  h,  ...,  r 

^^     ^^    ^^  ^   ~  [  j  =  h,  ...,  R 

-a.  .  -  a.,  -  a.,  +  A.  >  -2  (  ^  j  '^'I'y  ^ 

ij    ik    jk  k  -  J  1  /  J  /  k 


This  condition  is  still  true  if  gate  i  is  replaced  by  an  external 
variable.   When  the  complements  are  not  available: 

'wi      -  vr       -  a   +  X.  >  -2    ~^     j  =  h,  ...,  R 
^      *      ^^  ^  "        I     k  =  1,  ...,  R 

t     t     jk    k—      J  t=l,  ...,n 


When  the  complements  are  available: 

-t  -\  -"ok ^^5 ^-2  ^     i:i::::;R 

t     t     ik    k—      J  t=l,  ...,n 

Again  h  =  1  if  m  >  1  and  h  =  2  if  m  =  1.   (Only  h  =  2  was  used  in 
earlier  versions.   Also,  the  inequalities  for  the  case  when  the 
complements  are  available  were  omitted  previously  even  when  the 
complements  were  available. ) 

B.   NOR/NAND  Networks 

If  gate  j  is  a  NOR  gate  and  gate  k  is  a  NAND  gate  or  vice  versa, 
then  the  general  triangular  condition  holds  even  if  gate  j  has  other 
output  connections. 
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-a.  .  -  a.,  -  a.,  -  \.  +  \    >  -3 
ij    ik    jk    J    k  - 


■a.  .  -  a.,  -  oc.,    +  k.   -  \    >  -3 
xj    ik    jk    J    K  - 


i  =  h,  . . . ,  R 
j  =  h,  . . .,  R 
k  =  1,  ...,    R 


This  condition  is  still  valid  even  if  gate  i  is  replaced  by  an 
external  variable.   When  the  complements  are  not  available: 


t     t 


3  k 

t     t 


a 


■jk 


^j  +  \  >  -3 


a.,  +  X.  -  \    >  -3 
ok    J    k  -  ^ 


j  =  h,  • • • ^  R 

k  =1,  ...,  R 
J  /^  k 

t  =1,  . . . ,  n 


When  the  complements  are  available: 


■^  -  ^t  -  °jk  -  ^0  ^  \  >  -3 


■^t  -  ^t  -  «Jk  +  ^j  -  \  >  -3 


j  =  h^  • •  •  ^  R 

k  =  1,  ...,  R 

J  ^  k 

t  =1,  . . . ,  n 


As  before,  h  =  1  if  m  >  1  and  h  =  2  if  ra  =  1.   (Only  h  =  2  was 
used  in  earlier  versions  of  ILLODIE-AIF.   Also,  the  inequalities 
for  the  case  when  the  complements  are  available  were  omitted  pre- 
viously even  when  the  complements  were  available. ) 

C.   AKD/OR  Networks 

If  three  gates  are  connected  as  in  Fig.  C2,  then  the  network  is 
not  optimal. 


-a.  .  -  a.,  -  a.^  >  -2    i  =  h,  . . . ,  R 

13      IK     JK  — 

j  =  h,  ... ,    R 

k  =  1,  ...,  R 

i  ^  3    ^^ 


Qk 


Again  h  =  1  if  m  >  1  and  h  =  2  if  m  =  1.   (Only  h  =  1  was  used 
in  earlier  versions  of  the  program.  ) 

IV.  Output  Condition  -  NOR,  NOR/AIJD,  NOR/NAKD  and  AND/OR  Networks 

All  gates  except  the  output  gates  connect  to  at  least  one  gate. 

R 

k=l 
Mi 

V.  Connections  to  Outputs  Condition  -  NOR,  NOR/ AND  and  NOR/NAND  Networks 

Gates  which  are  connected  to  all  of  the  output  gates  are  not  con- 
nected to  any  other  gates. 

m         R 
-U  Z  a.  .  -   Z   a   >  -(raXU) 
j=l  ^^   k^+1  ^^ 
k;^i 

i  =  m+1,  .  .  . ,  R 

Here  U  =  n+R+2. 

VI.  Other  Conditions 

A.   NOR/NAND  Networks 

If  gate  i  has  only  one  input,  it  can  be  a  NOR  gate  (i,e., 
\  can  be  set  to  1  and  underlined). 

n   .    R 

Z  w^  +  Z  a   +  \  >  2       i  =  1,  . . . ,  R 
t=l  ^   k=l  ^^    ^  ~ 
k/i 

or 
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n   .    n   .    R 

Z  w  +  Z  V  +  Z 

t=l  ^   t=l  "^   k-1 


OL  .    +  \.    >  2 

Kl      1  — 


X  ~  X  f    %  .  ,  y    R 


k^i 


when  the  complements  are  available.   (Early  versions  of  the  program 
did  not  generate  any  inequalities  for  this  condition  when  the  com- 
plements were  available.   A  later  version  of  the  program  generated 
only  the  second  set  of  inequalities  even  when  the  complements  were 
not  available.   This  had  the  effect  of  generating  the  following 
inequalities  when  the  complements  were  not  available: 


n   .    R 
2Z  wj+  Z  a.+A..  >2 

t=.l  ^    k==l^^     ^- 


i  =  1,  ...,  R 


Mi 


These  inequalities  are  not  incorrect,  but  they  are  less  restrictive 
than  desired. ) 

B.   and/or  Networks 

1.   If  gate  i  connects  only  to  gate  j,   then  gates  i  and  j 
must  not  be  of  the  same  type. 


R 

X.  +  \.  -  a.  .  +    z  a,,  >  0 

kii 

R 

■\.   -  X.  -  a.  .  +    z  a.,  >  -2 
1   J    ij   ^^1  ik- 


k;^i 


>l 


J 


i  =  m+l,  . . . ,  R 
j  =1,  . .  • ,  R 

i  ^  J 
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The  first  inequality  keeps  both  gates  from  being  AND  gates 
while  the  second  inequality  keeps  both  gates  from  being 
OR  gates. 
2.   Each  network  has  at  least  one  AND  gate  and  at  least  one 
OR  gate. 


Z  X..  >  1 
i.l  ^- 


R 

-  Z  X.  >  -(R-1) 
i=l  ^- 


This  concludes  the  discussion  of  the  additional  inequalities  in- 
cluded in  the  current  versions  of  GEN-IE^.   However,  other  additional  inequal- 
ities can  be  added.   For  example,  consider  using  ILLODIE-AIF  to  design  networks 
under  fan-in  and  fan-out  restrictions.   The  following  inequalities  will 
incorporate  such  restrictions: 


n    .    R 
-  Z  w^  -   Z  a.  .  >  -FI  j  =  1,  .  .  . ,  R 

t=l  ^    i=l  ^J  - 


or 


i;^j 


n   .    n   .    R 
■  Z  w^  -  Z  v^  -  Z  a.  .  >  -FT        j  =  1,  . . . ,  R 
t=l  ^       t=l  ^   i=l  ^^ 


i^J 


when  the  complements  of  the  external  variables  are  available. 


8? 


and 


R 

■  Z  a.  .   >  -FO  i=  m+1,  .  . . ,  R 

3=1     "-'     - 


R 
-La...      >     -FOO  i  =1,  . . . ,  m 


^    k 

-  Z  w    >  -FOX  t  =1,  .  . . ,  n 

k=l  ^ 


^   k 
-  Z  V    >  -FOX  t  -  1,  .  . . ,  n 

k=l  ^ 


when  the  complements  are  available.   In  the  above  equations  FT,  FO,  FOO  and 
FOX  denote  the  maxim-am  fan-in  for  gates,  the  maximum  fan-out  for  non-output 
gates,  the  maximum  fan -out  for  output  gates  and.  the  maximum  fan -out  for 
external  variables,  respectively. 

The  user  must  be  careful  when  adding  other  inequalities  such  as 
those  for  fan-in  and  fan-out  restrictions  into  the  program  since  the  new 
inequalities  may  make  some  of  the  current  ones  invalid.   For  example,  the 
input  inequalities  given  in  I  for  NOR,  NOR/AM)  and  MD/OR  networks  do  not  hold 
in  general  when  the  fan -in  restriction  is  imposed. 
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CLLT( 

XS   ( 

"^UP  ( 

NCT2 

NCT6 

NFIXJ 


4000) 
2000) 
2000) 
2000) 


40 


50 


Ui=Ui+l 

8EG  =  CL 

Ei^JD  =  CL 

IE{  (END 

IF(  KS  . 
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S(JAD)  = 

JX  =  JX+ 

DO  I 

SUP 

IF  IC 

PS(S 

(PS(S 

( N  S  (  S 

(CLSP 

CLSPCNCT 

NCT3=SUti 

CLSP(NCT 


IF 
IF 
IF 


PT 
LT 

-ilE 
LT. 
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CONTINUE 
RETURN  1 

X(L)=0 CHANGE  POS  SUM 
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IF(  CLCF( I)  )   240,300,260 
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GU    TO    318 

PS( SU8)=PS( SU8)-CLCFIK) 
CONTINUE 
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IF  (COL(NBEG)  .EJ.  0)  NBEG=NBEG+1 
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(20) ,    CNKGT (20) 
T,    CNKGT,     lEND, 


4000) 
2000) 
2000) 
2000) 


/I. 2, 3, 4, 
2, 3. 4, 5 

CEND) 


NIPT 

U4 

N 

NCX 

COL  (20000) 

RWLT(  4000) 


4000) 
2000) 
2000) 

200C) 


CLVT( 
OF   ( 
S    ( 
TP   ( 
NCT3 
NGETYS 
RWLS(  4000) 
5,6     / 
/ 


T=l,NGGT 

.EQ.  0)  GO  TO  25 
=1, lEND 
EO.  ISGT( AA) )  GO  TO  110 


NOPT 

JAO 
LTS 

RWCF(200G0) 

YS   (  4000) 

CLSP(  4000) 

XHAT(  2000) 

SJ-^K  2000) 

NCT4 
NNEWX 


25    CONTINUE 

SRT=DINP*(GT-l) 

IF     ($L(ONE)     .LE.    N»     SRT=SRT»2 

DO    100    K=l,NIPT 

SUB=$P1GT,K) 

IF     (S(SUB)     .Eg.    0)    GO    TO    100 

I I=SRT+K 

IF     (NS( I  I »     .GE.    01    GO    TO    21 

BEG=RWPT  (II  ) 

EN0=8EG+RWLT( I I)-l 

GO  TO  29 
27     IF     (NSdI+DINP)     .GE.    01    GO    TO    100 

IF     (iL(ONE)     .GT.     N    I     GO    TO    100 

8EG=RWPT( II+OINP) 

END=BEG+KWLT(  II«-DINP)-l 

29  CONTINUE 
COONT=0 
ISCONT=0 
SMCT=5000 
ISB=5000 
WB=N+l 
ODR=100 

3i>    DO    90    L=BEG,END 

SUB  =  COHL) 

IF     ( SUB    .LE.     0)     GO    TO    90 

IF  ISISUB)  .NE.  0)  GO  TO  90 

IF  (SUB2(SUB)  .NE.  01  GO  TO  30 

IF  (TP(SUB)  .NE.  LLT)  GO  TO  36 

KS  =  l 

JAD=SUR 

XS( JX)=JAO 

CALL  CGPSNS  ( Z X , &200, &220  ) 

RETURN  I 
16  CONTINUE 

COUNT=CUUNT+I 

IF  (ODR  .LE.  VC$)  GO  TO  90 

Wrt=SUB 

OOK=VC$ 

GO  TO  90 

30  IF  ( lEND  .EG.  0)  GO  TO  ^3 
on  ^0  AA  =  l, I  END 

IF  (SUBKSUBJ  .NE.  ISGT(AA))  GO  TO  ^0 

ISCUNT=ISC0NT+1 

IF  ( ISR  .LT.  SUB)  GO  TO  90 

IS8=SUB 

GO  TO  90 
40  CONTINUE 
45  CONTINUE 

C0UNT=C0UNT+1 

A2=$AI SUBl (SUb) ,  SU82(SU3)» 

IF  (  S(A2))  50,  50,  60 
50  PP=$P(SUBl(SUB)  ,K) 

IF  {S{PP>  .NE.  0)  GO  TO  55 

IF  (ODR  .LE.  G$C$)  GO  TO  90 

Wb=SUB 

ODR=G$C$ 

GO  TO  90 
55  IF  (ODR  .LE.  GC$)  GO  TO  90 

WB=SUB 

OOR=GC$ 

GO  TO  90 


60  IF     (ODR    .LE.    G$l    GO    TO    90 
UOk=Gt 

WB=SUB 
90    CONTINUE 

6t)     IF     (WB    .L6.    N)     GO    TO    95 
70    IF     (ISB    .EQ.     bOOO)    RETURN 

JAD=ISB 

GO  TO  1^0 
95  CONTINUE 

IF  (ISCONT  .GE.  l>  COUNT=COUNT+l 

IP  (TYPE-OOR)  152,151,100 
151  IF  (SMCT  .LE.  COUNTJ  GO  TO  100 

GO  70  153 
15?  TyPE=ODK 
153  S^CT=CCUNT 

JAn=WH 
100  CONTINUE 
110  CONTINUE 

IF  (JAD  .EQ.  0)  GO  TO  201 

KS  =  1 

XS( JX)=JAD 

CALL  CGPSNSCZX,  &200,  £220) 

RETURN  1 

120  KS=1 

X<^  (JX)=-JAO 
GTKT=JAD 

CALL  CGPSNS  (ZX,  f.?00,  &220I 
IF  ( iL(ONE)  .GT.  N)  RETURN  1 
JAl)  =  $L(SUB1(  JAD)  J 
IF  (S(JAD)  .NE.  0)  GO  TO  210 
GTFD=$L(SUB2(GTKT )) 
IF  (S(GTFO))  51,51,61 
51  KS=1 

GO  TO  121 

61  KS=-1 

121  XS(JX)=JAO 

CALL    CGPSNS     ( ZX ,     L2Q0fL2ZO) 
210    RETURN     I 

200  RETURN 

201  IF     ($L(ONF)     .GT.     N    )     RETUi-^N    2 
DO    205    I=1,N0GT 

AB=$L{ I > 

IF     (S( AB)     .NE.    0)    GO    TO    205 

JAO=AB 

KS  =  -1 

XS( JX»=JAD 

CALL    CGPSNS     (ZX,     G200,     &220) 

RETURN  1 
205  CONTINUE 
220  RETURN  2 

END 


INTEGER   FUNCTION  FUNCT*2(XX,  Y  ) 
IMPLICIT  INTEGER*?  (A-T,V-Z,t) 
INTEGER*^  XX, Y,  NV,NG,  NRG,  NIPT 
INTEGER*^  NOGT,  NOVR,  NIPT,  NOPT 
COMMON     NOVR        ,    NOGT 

I      ,     02  ,    U3 

COMMON     ZBAR         ,     M 

I       ,     KS  ,    JX 


f 

NIPT 

» 

NOP 

t 

U4 

f 

N 

f 

JAD 

t 

NCX 

f 

LTS 

I         , 

JTS 

,          CLCF(20000), 

COL  (20000), 

RWCF(20000) 

3      . 

ROW  (20000)i 

,          RWPT( 

4000) 

f          RWLK  4000), 

YS   ( 

4000) 

*      f 

PS   C  4000), 

NS   ( 

4000) , 

r    CLVTi  4000), 

CLSPC 

4000) 

5      , 

CLPT{  20001 

CLLT( 

2000), 

OF   {  2000), 

XHAT( 

2000) 

6      , 

X    (  2000) 

,          XS   ( 

2000) 

r    S    (  2000), 

SUBK 

2000) 

7      , 

SUB2(  2000), 

»    SUP  ( 

2000), 

»    TP   (  2000) 

COMMON 

CLLS(  2000), 

NCT2 

NCT3         , 

NCT4 

I 

NCT5         , 

,          NCT6 

,          NGETYS       , 

NNEMX 

2 

NUNDLX      , 

NFIXJ 

,          RWLS(  4000) 

DATA   CP 

/•CP« 

/ 

NV=XX 

NG=Y 

NIPT=2**NV 

NBG=NIPT*(NG-l) 

WB=1 

V8=WB+NV*NG 

IF  (NCT5  .NE.  CP)  VB=WB 

6ETA=V8+NV*NG 

ALPHA=BETA+NBG*NG 

PB=ALPHA+(NG-1 )*NG 

LAMDA=P6+NIPT*NG 

GAMA=LAMDA+NG 

FUNCT=NflG 

RETURN 

ENTRY     i^i I ,K) 

$W  =WB+NV*(  I-l)«-K-l 

RETURN 

ENTRY    $V(I,K) 

iV       =V8+NV*(  I-l  )+K-l 

RETURN 

ENTRY     i3( I , J,K) 

JJ  =  J 

IF     (J.GT.I)     JJ=JJ-1 

$H  =BtTA*NBG*(  I -I ) ♦N I PT* ( J  J- 1 ) +K- I 

RETURN 

ENTRY    $ACI , J) 

JJ  =  J 

IF        (J     .GT.     I)     JJ=JJ-l 

$A  =ALPHA+ (NG-1 )♦ ( I-l ) 

RETURN 

ENTkY     SP(I,K) 

tP  =PB+NI PT«( I-l) +K-1 

RETURN 

FNTRY     $L ( I ) 

$L  =LAM0A+I-1 

RETURN 

ENTKY     $G(I,J,K) 

JJ  =  J 

IF     (J.     GT.     I)     JJ=JJ-l 

$G  =GAMA+NbG*i I-l )+NIP" 

RETURN 

END 


♦  J  J-1 


'T«(  JJ-i  )*K-l 


I 


SUBROUTINE  GEOF  (UB) 
IMPLICIT  INTEGER+2  (A-T,V-Z,S) 
INTEGER*4  NUGT,  NHVR,  NIPT,  NOPT 

INTtGER*4  TIME  ,ITIMEZ  ,  CT  1 ,  CT2  ,  CT3  ,  CT4  ,  CT5 , CT6 , C T 7 ,C T d 
INTEGER*4   KT I M EZ , JT IMEZ , I T I  ME 

COMMON     NOVR        ,    NOGT        ,    NIPT        ,    NOPT 
»     U2  ,    U3  ,    U4 


CUMKON 


CGMMflN 


1  iO 


30U 


10^0 


15 


25 


35 


30 

40 

1010 


50 
1030 


DIMFNS 
RrAf  1 
IF  (16 
CALL    P 

MND=JX 

on  5   I 

li-  (X( 
JAD=X( 
CALL     P 

CONTIN 
FORMAT 
CT2  = 
CT3  = 
CT4  = 
CT5  = 

cr6  = 

CT7  = 

CT3  =0 
NFS  =  0 
U=KT  IM 
'J=U-5  5 
IJ=U/10 
PRINT 
FORMAT 
ITI,VE  = 
CALL  S 
UTIN:E  = 
DO  8  0 
DO  7  5 
DO  70 
CT2  = 
hi 'P  VAT 
CALL  E 
CGMTIN 
GOTL;  5 
CONTIN 
CT4  = 
GOTO  6 
CONTIN 
CT5  = 
GOTT  5 
CONTIN 
CT7  = 
FORMAT 
CALL  F 
GO  TO 
CONTIN 
FORNAT 
CALL  N 
CTl  =  i|2 


lU 
30 
AR 
Rt 
-I 
=  1 
I  ) 
I) 
HB 
UE 
( 
0 
0 
0 
0 
0 
0 


ZtiAR 

KS 

JTS 

ROW  (20000) 

PS   (  -^000) 

CLPT(  2000) 

X    (  2000) 

SUB2I  2000) 

CLLS(  2000) 

NCT5 

NUNULX 

N  DUMMY(2500) 

,ZrtAR 
.LT.  0 )  GO  TO 

SET  (ZX,  &300) 


JX 

CLCF( 20000) 

RWPT(  4000) 


NS  ( 
CLLT( 
XS  ( 
SUP  ( 

NCT2 
NCT6 
NFIXJ 


4000) 
2000  ) 
2000) 
2000) 


N 

NCX 

COL  (20000) 

RWLT(  4000) 

CLVT(  4000) 

OF   (  2000) 

S    (  2000) 

TP   (  2000) 

NCT3 

NGtTYS 

RWLS(  4000) 


JAO 
LTS 

RWCF(20000) 

YS   (  4000) 

CLSP(  4000) 

XHAT(  2J00) 

SUBK  2000) 

NCT4 
NNFWX 


80  0 


,NND 
.LE.  0)  GO  TO  5 

LP(ZX,t5,C5) 

115) 


£Z 
0 

0. 
30 
(  • 
36 
TI 
KT 
KK 
II 
LK 
CT 
( 
XT 
UE 
0 

UE 
CT 
0 

ue 

CT 
0 

UE 
CT 

( 
IX 
60 
UE 

( 


(•CANC  ) 


o,u 

OTI 
000 
MFZ 
I  ^E 


2  +  1 
//• 

DX 


ME  IN  SECONOS  I  S  •  fi-H,Z     ) 

0 

(  ITIME) 

ZI  0) 

I  ,20 

1,  2 

1000 

ENTER  EXTDX' ) 
(ZX,&25,&30,&35) 


4+1 


5  +  1 


7+1 

//•ENTER  FIXJ') 

J(ZX,&70,&15) 


//•ENTER  NFWX*  ) 
X  (ZX,   f.60) 


CT3    =    U4 

CT6=U3 

PRINT    310,    NCX    ,CT1  ,    CT2,CT3,    CT4,    CT5,CT6  ,    CT7     ,CT8 

OF'.llX,     'GETYSSex,     'NO.     ITER*     ,7X,     •    CYCLES       •     ,5X 
UNFS«,5X,     'EXTDX    FSBL«,6X,     •       CGPSNS    •,9X,     'FIXJ*     ,10X 
2,'0FBND«,    /     19,8115     ) 
ST=1 
PTS=0 
N8  =  0 

CALL  PRNTCXHAT,  ST,  N8,  PTS) 
60  CONTINUE 
1020  FORMAT  (//'ENTER  UNOLX') 

CALL  UNDLXCZX,665,&90,650) 
65  CONTINUE 
70  CONTINUE 

U=KTIMEZ(0)-UTIME 
U=U/100. 
PRINT  300, U 
ST=0 
NR  =  0 
PTS  =  1 

CALL  PRNT  (XHAT,  ST,  NB ,  PTS) 
CTUU2 
CT3  =  UA 
CT6=U3 

PRINT  310,  NCX  ,CTl     ,  CT2,CT3,  CT4,  CT5,CT6    ,  CT7  ,CT8 
7b  CONTINUE 

BEG  =  NI  PT*NOPT-H 
ENO  =  jx 

00  UO  IN=BEG,ENO 

IF  (  XS(IN)  .EQ.  X(IN))  GO  TO  100 
DUMMY  (IN)  =  0 
GO  TO   110 
100  DUMMY  (IN)  =  1 
110  CONTINUE 

nUMMY(JX)  =  0 

PUNCH  120,  ( XS( IN),DUMMY( IN) ,  IN=dEG,tNDJ 
120  FORMAT  (  12(15,11)) 

PUNCH  350,ZBAR 
350  FORMAT  (215) 
80  CONTINUE 
90  CONTINUE 

U=KTIMEZ(0)-UTIME 
U=U/100. 
PRINT  300, U 
CT1=U2 
CT3  =  U-V 
CT6=U3 
ANCX=ZBAR-1 

PRINT  310,ANCX  ,CT1    ,  CT2,CT3,  CT't,  CT5,CT6    ,  CT7  ,CT9 
360  FORMAT  ( // /5X , • P SEU= • , I  1 5 , 5X , • ZERO= • , I  1 5 ) 

GO  TO  1 
800  RETURN 
END 

SUBROUTINE  GETYS 

COMPUTE  NEW  YS(I)  BY  APPLYING  COLUMN  JAO 

IMPLICIT  INTEGER*2  (A-T,V-Z,$) 

INTEGER*^  NOGT,  NOVR,  NIPT,  NOPT 

COMMON     NOVR         ,    NOGT         ,    NIPT         ,    NOPT 


1  t 

COMMON 
L 
2 
3 
4 
!> 
6 
7 


U2 

ZBAR 
KS 
JTS 
ROU     ( 
PS   I 
CLPT( 
X    ( 
SUB2( 
CLLS( 
NCT3 
NUNDL 
GETYS 


COMMON 
1 

Z  f 

42  FGRMATC 

U2=U2*l 

NCF=CLPT 

NCFFND  = 

1F(CLLT( 

UO  10  1= 

J  =  ROW( I  ) 

YS( J)=YS( J)+CLC 
10  CONTINUE 
43    FCRMATC 
^0  RETURN 

END 


(JAD+l) 

NCF  *■ 

IJAD+l ) . 

=NCF,NCF 


GETYS 


20000) 
4000) 
2000) 
2000) 
2000) 
2000) 


ENTER 

CLLT( JAD 

EJ.O)     GO 
E!>iD 

F( I)*KS 

EXIT     •     , 


U3 
M 

JX 

CLCF(2U000) 
RWPTl  4000) 


NS  ( 
CLLT( 
XS  ( 
SUP  ( 

NCT2 
NCT6 
NFIXJ 


+  1)  -I 

TO  20 


l>i) 


4000) 
2000) 
2000) 
2000) 


U4 

N 

NCX 

COL  (20000) 

RWLT(  4000) 

CLVT(  4000) 

OF   (  2000) 

S    (  2000) 

TP   (  2000) 

NCT3 

NGETYS 

RWLS(  4000) 


JAD 
LTS 

Rv<CF(20000) 

YS   (  4000) 

CLSP(  4000) 

XHAT(  2000) 

SU81(  2000) 

NCT4 

NNEWX 


SUBROUTINE  GTINEO 

GENFPATt  INcQUALlT  lES  FOR  NOR  NETWORKS. 

IMPLICIT  INTEGEK*2  (A-T,V--Z,i) 

CALL  Q1201 

CALL  Q12Q2 

RETURN 

END 


cc 


SUBR 
IMPL 
INTfc 
COMM 


OUT  I 
ICIT 

GER* 
ON 


NE   GTPSNS 

INTEGER*2  (A-T,V-Z,$) 
4  NIUGT,  NOVK,  NIPT,  NOPT 


C.IMMHN 


COMi^ON 


CALC 

NCT2 

NCTi 

DO 

CLSP 

PS(  I 

NS(I 

NBEG 

NrNO 

IF  ( 


ULAT 

=  1 

=  1 

50 

(  I  ) 

)  =  YS 

)=YS 
=  R 
=  N 

COLI 


NCVR 

U2 

ZBAR 

KS 

JTS 

ROW  (20000 


PS   < 

CLPT( 
X    ( 

SUB2( 
CLLS( 
NCT5 
NUNOLX 


4000) 
2000) 
2  000) 
2000) 
2000) 


POSITIVE  SUM 


NCiGT 

U3 

M 

JX 

CLCF(20000) 

PWPTC     ^000) 


NS  ( 
CLLT( 
XS  ( 
SUP  ( 
NCT2 
NCT6 
NFI  XJ 


4000) 
2000) 
^000) 
2000) 


NIPT 

U4 

N 

NCX 

COL  (20000) 

KWLT(  4000) 

CLVT(  4000) 

OF   (  2000) 

S    (  200U) 

TO   (  2000) 

NCT3 

NGETYS 

RWL?(  4000) 


NEGATIVE  SUH 


1  =  1, M 
=    0 
(I  ) 
(I  ) 

UPT(  I  ) 

BEG  +  kWLTI  I  )  -  I 

NBEG)  .NE.  0)  GO  TO 


NOPT 

JAD 

LTS 

KWCF (20000) 

YS   (  4000) 

CLSPl  4000) 

XHAT(  2000) 

SJriK  2000) 

NCT4 
NNErtX 


10 


10 


20 

30 
40 


50 


NBEG=NBEG+1 

IF(NBEG    .GT.    NENDI       GO    TJ    50 

on    40    K=NBEG,NENO 

IF     (S(COL(K)).NE.OI    GO    TO    40 

IF     (KWCF(K))     30,40,     20 

PS( I l=PS( I )+RWCF(K) 

GO    TO    40 

NS(  I)=NS(I  )«-RrtCF(K» 

CONTINUE 

IF     (NS (I  )     .GE.    OJ     GO    TO    50 

CLSPtNCT3)=I 

NC  T  3=  f 

CLSP(NCT3)=I 

CONTINUE 

RETURN 

ENO 


C 
C 
C 
C 
C 
C 
C 
C 
C 


FUNCTIONS 


SUBROUTINE     INTL 

NOVR     :     NO.     CIF    VARIABLES    OF    THE     BOOLEAN 

NOGT     :     NO.     OF       GATES     IN    THE    NETWORK 

NIPT    :     NO.    OF     INPUT    VECTORS 

NOPT  :  NO. OF  OUTPUT  FUNCTIONS 

BEGINNING  POINTER  FOR  CHAIN  CHECKING 
ENO  POINTER  FOR  CHAIM  CHECKING 
NO.  OF  TYPES  OF  GATES  IN  THE  NETWORK 

NCT5=CP  IMPLIES  COMPLEMENTS  OF  VARIABLES  ARE  AVAFLAdLE 
CLLS,KrtLS,NGFTYS,NNEWX,NUNDLX,NFIXJ,ARF  NOT  USED  IN  PROGRAM 
INTEGER*2  (A-T,V-Z,$) 
NOGT,  NOVR,  NIPT,  NOPT 


NCT2 

NCT3 

NCT4 

NCT5 

CLVT 

IMPLIC  I  T 

INTEGER*4 

COMMON 

t 
COMMON 

t 

» 


COMMON 


OIMENSION 


NOVR 
U2 

ZBAR 
KS 
JTS 

ROW  (20000) 
PS  (  40001 
CLPT(  2000) 
X  (2000), 
SUB2(  2000) 
CLLS(  2000) 
NCT5 
NUNDLX 
INPT(270) 


NOGT 
>J3 
H 

JX 

CLCF( 20000) 
KWPT(  4000) 
NS  (  <,000) 
CLLT<  2000) 
XS  (  2000), 
SUP  (  2000) 
NCT2 
NCT6 
NFIXJ 


NIPT         , 
U4 

N 

NCX  , 
COL  (20000), 
RWLT(  4000), 
CLVT(  430J), 
OF  (  2000), 
S  (  2000), 
TP  (  2000) 
NCT3  , 
NGETYS  , 
R><LS(  4000) 


FWUI VALENCE  ( TYI N , JX ) , ( I NPT ( 1 ) , XS ( I ) ) , I NOKGT , NCT4 ) 


/•O 

IB(  12)  , 
1=1,20) 

1=1,20) 


IC( 12) 


DATA    D,W,C,8L,V 

DATA    CARD/«CD'/ 

OIMENSION       UA(20) 

READ       10,       (UA(I) 
10    FORMAT(20A4) 

PRINT    20,        (UA( I ) 
20    FORMAT (1HI,20A4) 

READ    3  0,N0VAR,N0GT,NOKGT,TYIN,NCT3,N0PT 
30  FORMAT     ( 31 5 , 3X , A2  ,3X , A2 , I  5  ) 

NOVR=NOVAR 

NIPT=2**N0VR 

ZZXXYY=FUNCT( NOVR, NOGT) 
C      INITIAL  SET 

DO  50  1=1,2000 

CLLTd  )  =  0 

0F(  I)  =  0 


NJPT 

J  AD 
LTS 

RWCF(200J0) 
YS  (  4000) 
CLSP(  4000) 
XHAT(  2000) 
SUBK  2000) 

NCT4 

NNEWX 


50 


55 


60 
HO 


71 
70 


72 
73 


no 


9b 
iOJO 


100 


10^ 


104 


410 


420 

430 


10« 


CONTINUE 

DO    55    1=1,4000 

RWPT( I )=0 

RWLTd  )=0 

CONTINUE 

1)0    60     1  =  1,20000 

RWCK  I  )=0 

COL( I )=0 

CONTINUE 

J       =1 

TWON    =    2**NUVAR 

RWPT{l)=l 

K    =     I 

TOTl.VR=0 

CALL    TBPT(NOVR,NO 

RETURN 

ENTRY    INTLC 

READ    OBJECTIVE    FU 

READ    70,     (IB(I ) ,1 

FORMAT     (1814) 

DO    73     1=1,9 

IF     (ICdl)     90,72, 

0F(  IC(  I  )  )=IB( I) 

CONTINUE 

GO    TO    71 

READ    MATRIX    AND    G 

READ    EXPANDABLE    D 

R=l 

T  =  ?7 

REAL 

FORMAT ( 

RR=P+2 

00    100    I =    RR  ,T,3 

IF (  INPT(  I  )  .LT.     0 

CONTINUE 

R=R+27 

T=T+27 

GOTO    9  5 

trJTRY    INTLP(T) 

IF{     INPTd)     .£Q. 

GOTO    108 

NON    -    EXPANDABLE 

MAXVR=-100 

DO  420  1=1, T, 3 
IF(  INPT(I+2)) 
RrtCFlKl  =  INPT 
COL     (K)     =     INPT 

CLLTICOL (K)+l )=CL 

IF  (COLIK)  .GT.  M 
RWLT(J)  =  RWLT 
K       =    K-H 

J    =     J     +1 

PvgPT(  J)     =    K 

IF     (MAXVR    .GT.     TO 

IF     ( TYIN.EU.CARDI 

RETURN 

DO       320    JJ=1,TW0N 

MAXVR=-100 

on       300       I=l,T,3 

IF(     INPT( 1+2) .LT. 

IF     (INPT(l)     .EQ. 


1000,     (INPT( 
9(  Al  ,13,1 


GT) 


NCT  ION 

C  I  I  )  ,  I  =  1  ,  9  J 


72 


ENERATE    t<WCF,CbL 
ATA    CARDS 


I)     ,     1=    R,T» 
4)  ) 


)    GOTO    10  2 


BL     )    GOTO     104 

ROWS 


430,410,410 
(     I  +  l     ) 
(     1+2     ) 

LT(COL(K)+l )+l 
AXVR)    MAXVK=COL(K) 
(J)     +1 


TLVR»    TOTLVR=MAXVR 
G(3    TO    90 


0    )    GOTO    110 
V)     GO    TO    1210 


IF(  INPT(I)  .EO.  W     IGOTO  210 
IF{  INPTCII  .EQ.  0     )GOTO  UO 
IF(  INPTdJ  .NE.  C  I  GOTO  400 
C      COLUMN  ENTRIES 

RWCF(K)=  INPT(I+1) 
COL  (K)=  INPT(H-2) 
CLLT(C0L(KK-1  )=CLLT(C0L(K)4-l)*l 
IF  (COL(K)  .GT.  MAXVR)  MAXVR=COL(K) 
RWLT(J)=  RWLT(J|+1 
K=  K  +  I 
GOTO  300 
:       DIAGONAL  ENTRIES 
no  RrtCF(K)  =  INPT(  I*  1  ) 
COL  (K)  =  INPT(I*2) 
CLLT(C0L(K)«-I)=CLLT(C0L(K)+l)4-l 
IF  {COL(K)  .GT.  MAXVR)  MAXVR==COL(  K I 
RWLT( J )  =  RWLT( J)+l 
INPT( 1+2)  =  INPT(  H-2»+I 
K  =  K+1 
GOTO  3  00 
:       V  SECTION 
1210  L=  JJ-1 

CHCK  =  2**(N0VAR-1) 
DO  1230  KK=I,NOVAR 
L  =  L-  CHCK 

IF  (  L  .GE.  0  )  G0T01230 
RWCF(K)  =  INPT(I+1) 
COL  (K|  =  INPT(  1*2)  ♦  KK-I 
CLLT(C0L(K)*1)=CLLT(C0L(K)«-1  )4-l 
IF  (COL(K)  .GT.  MAXVR)  MAXVft=COL(KI 
K=  K  +  l 

RWLT( J)  =  RWLT( J)+l 
1220  L  =  L+CHCK 
1230  CHCK  =  CHCK/2 
GO  TO  300 
W  SECTION 
210  L=  JJ-1 

CHCK  =  2**(N0VAH-1) 
DO  230  KK=l,NOVAR 
L  =  L-  CHCK 

IF  (  L  .LT.  0  )  GOTO  2Z0 
RWCF  (K  )  =  INPTi  K-1) 
COL  (K)  =  INPT(  1+2)  +  KK-1 
CLLT(COL(K)«-l  J=CLLT(C0L{K)+1  H-1 
IF  (COUK)  .GT.  MAXVR)  MAXVR=COL(K) 
K=  K*-l 

RwLT(J)  =  RwLT( J)*l 
GOTO  230 
220  L  =  L+CHCK 
230  CHCK  =  CHCK/2 
300  CONTINUE 
310  J  =J«-1 

RWPT(J)  =  K 

IF     (MAXVR. GT.TOTLVR)    TOTLVR=MAXVK 
320    CONTINUE 

IF     ( TYIN.EQ.CARD)    GO    TO    90 
RETURN 
400    M    =    J-1 
N=TOTLVR 
42     PRINT    40,N,M     ,NOVAR     ,K 
40    FCRMATdH     a7HN0    OF    VARIABLES    = ,  15 .  5X,  20HNC    OF     INEQUALITIES    =,Ib, 


I    5X, 

RETURN 
END 


•NO  OF  INPUT  VARIABLES  •,  I!>f5X,'N0  OF  NGN  ZERO  COEF  =«,I6) 


10 


30 


40 


50 


SUBROUTINE 
IMPLICIT  I 

INTEGfcR*4 
COMMON 

I      t 

COMMON 
1 
2 
3 

5 
6 

7 
COMMON 

I 

2 

DIMENSION 
DO  b    AA=1, 
ISGT{AA)=0 
CNKGT( AA)= 
CONTINUE 
DO  10   1=1 
CNKGT( I )=I 
END1=NCPT 
EN02=0 
BEG=N0PT+1 
DO  50   I=B 
DO  30   J=l 
IF   (J  .EQ 
SU8  =  $A( I  ,J 
IF   (SI  SUB 
CONTINUE 
ENu2=END2+ 
ISGT(END2) 
GO  TO  50 
END1=ENDI+ 
CNKGT( ENDl 
CONTINUE 
RETURN 
END 


LSTISGdSGTt  CNKGT,  ENQ2,  ENDIJ 


NTEGER*2  (A-T,V 

NOGT,  NOVRt  nip 

NOVK 

U2 

ZflAR 

KS 

JTS 

ROV^  (20000) 


4000) 
2000) 
2000) 
2000) 
2000) 


PS   ( 

CLPT( 

X    { 

SU32( 

CLLS( 

NCT5 

NUNOLX 

ISGT(20) 
20 


,NOPT 


-Z,$) 

T,  NOPT 

NOGT 

113 

M 

JX 

CLCF(20000) 

RWPT(  4000) 


NS  ( 
CLLT( 
XS  ( 
SUP  i 
NCT2 
NCT6 
NFIXJ 
CNKGT(20) 


EG, NOGT 

,NOGT 

.  I)  GO  TO  30 

) 

)  .GT.  0)  GO  TO  40 

1 
=  1 

1 
)  =  I 


4000) 

2000) 
2000) 
2000) 


NIPT 
U4 

N 

NCX 

COL  (20000) 

RWLTI  4000) 

CLVT(  4000) 

OF   (  2000) 

S    (  2000) 

TP   (  2000) 

NCTJ 

NGETYS 

RWLS(  4000) 


NOOT 

JAD 
LTS 

RWCF(20000) 

YS   (  4000) 

CLSPt  4000) 

XHAT(  2000) 

SUBU  2000) 

NCT4 

NNEv^X 


MAIN  PROGRAM 

IMPLICIT  [NTEGER*2  (A-T,V-Z,$) 
INTfcGFR*4  NOGT,  NUVR ,  NI^T,  NOPT 
TIME  ,  ITIMEZ 

KTIMEZ ,JTIMEZ,ITIME 
NUVR         ,    NOGT 
U2  ,    U3 

ZDAR         ,    M 
KS  ,    JX 

JTS  ,  CLCFI20000) 
ROW  (20000),  RWPT(  4000) 
PS  (  4000),  NS  (  4000) 
CLPT(  2000),  CLLT(  2000) 
X    (  2000),    XS   (  2000) 


INTEGER*4 
INTFGER*4 
COMMON 
, 
COMMON 


NIPT  , 

U4 

N 

NCX  , 

COL  (20000), 

RWLT(  4000), 

CLVT(  400U), 

OF   (  2000), 

S    (  2000), 


NOPT 

JAD 

LTS 

RWCF( 20000) 

YS   (  4000) 

CLSP(  4000) 

XHATI  2000) 

SUBU  2000) 


20 


10 
15 


7      ,     SUB2(  2000), 

SUP  (  20001, 

TP   (  2000) 

COMMON     CLLS{  2000), 

NCT2         , 

NCT3         , 

NCT4 

I      ,     NCT5         , 

NCT6         , 

NGETYS      , 

NNEWX 

2      ♦     NUNDLX      , 

NFIXJ        , 

RMLS(  4000) 

EQUIVALENCE  (TYIN,JXI 

DATA  CARD/»CD'/ 

CONTINUE 

ITIME=360000 

CALL  STIMEZ(ITIME) 

CALL  INTL 

IF  (TYIN.EQ.CARL;»  GO  TO 

10 

CALL  GTINEQ 

GO  TO  15 

CALL  INTLC 

CALL  TAbPT 

CALL  CLFVAR 

CALL  RWTUCL 

CALL  GEOF  (UB) 

GO  TO  20 

END 


SUBROUTINE        MAX    Z     ( NUM, 
IMPLICIT     INTEGER*2     (A-T, 


lNTbGER*4 
COMMON 

COMMON 
, 
f 

f 
f 
* 


NJGT,     NOVR,     NI 

NOVR 

U2 

ZBAR 

KS 

JTS 

ROW    (20000) 


COMMON 


PS  ( 
CLPT( 
X  ( 

SUB2( 
CLLS( 
NCT5 
NUNOLX 


4000) 
2000) 
2000) 
2000) 
2  000) 


♦,*) 
V-Z,S) 

PT,     NOPT 
NOGT 
U3 
M 

JX 

CLCF(20000) 
RWPT(     4000) 


IF     (NCX-Z^AR+1) 

15  CONTINUE 
RETURN 

16  DO  17  K=l,  N 
IF  (OF(K)  .EO. 
IF  (  S(K)  .NE. 
XS( JX)=-K 
NUM=-1 

JAD  =  K 
KS  =  -1 
CALL  CGPSNS 

17  CONTINUE 

19  RETURN  1 

20  RETURN  2 
END 


19,  l(j 


NS  ( 
CLLT( 
XS  ( 
SUP  ( 
NCT2 
NCT6 
NFIXJ 
15 


0) 
0 


GO  TO 
)  GO  TO 


17 
17 


(ZX,tI5,£20) 


4000) 
2000) 
2000) 
2000) 


NIPT 

U4 

N 

NCX 

COL  (20000) 

RWLTC  400J ) 


CLVT( 

OF   ( 

S    ( 

TP   ( 

NCTi 

NGtTYS 

RWLSI  4000) 


^000) 
2000) 
2000  ) 
2000) 


NOPT 

JAD 

LTS 

RrtCF( 20000) 

YS   (  4000) 

CLSP(  4000) 

XHAT(  2000) 

SJBK  2000) 

NCT4 

NNFWX 


SUBROUTINE   NEWX  (ZX,   *) 
COMPARE  ZrtAR  AND  OBJECTIVE  FUNCTION 
IMPLICIT  INTEGER*2  (A-T,v-Z,$) 
INTEGER+4  NOGT,  NOVR,  NIPT,  NOPT 
COMMON     NOVR         ,    NOGT 
I      f     U2  ,    U3 


TO  FIND  A  NEW  SOLJTION 


NIPT 
U4 


NOPT 


COMMON 


ZBAR 

KS 

JTS 

ROW  (200001 


12 


12 


PS  t  4000) 
CLPT{  2000) 
X  (  2000) 
SUB2(  2000) 
COMMON  CLLS(  2000) 
NCT!j 
MUNDLX 
A2  f-ORMAT(«  MEV^X   ENTER') 

if(ncx.(;e.zbar  )  go  to 

ZBAR  =  .MCX  +  l 
DO  10  K=1,N 
XHAT(K)=0 

10  CONTINUE 

DO   II   L=1,N 
IF(X(L).LE.O)   GO  TO  I  I 
J=X(L) 
XHATC J)=l 

11  CONTINUE 
LTS=1 

Rt  TllrJN 

CUNTINUE 
I    Fi.)RMAT(«  NEWX  EXIT   •  , 
LTS=1 
RETURN  1 
ENO 


M 
JX 

CLCF(200aO) 
RWPT(  4000) 


NS  ( 
CLLT( 
XS   ( 

SUP  ( 

NCT2 

NCT6 

NFIXJ 


13) 


4000) 
2000) 
2000) 
2000) 


N 

NCX 

COL  (20000) 

RWLT(  4000) 

CLVT(  4000) 

OF   (  2000) 

S    (  2000) 

TP   <  2000) 

NCTi 

NGETYS 

RWLS(  4000) 


JAD 
LTS 

RWCF(20000) 

YS   (  4000) 

CLSP(  4000) 

XHAT(  2000) 

SJBK  2000) 

NCT4 
NNEWX 


SU 
IM 

IN 
C  J 


BKOU 
PLIC 

TfcGE 
MMON 


Ct.'MMON 


COMMON 


IF 

TH 
TO 
01 
OA 
IF 
IF 
FI 
RE 
EN 
SC 
IF 
1  = 

on 


» 

OOF 
E  GA 

GAT 

MEN'", 

TA  W 

(T 

(S{ 
NP 

Gl  =  l 
01  =  1 
R(  1) 
(8 
SCR( 

10 


TINE 
IT  IN 

R*4  N 
N 
U 
Z 
K 
J 
R 
P 
C 
X 

s 

C 

N 

N 

INTE 

TES  P 

E  J 

ION 

f  B,A, 

P(  JAD 

JAD) 

THE 


PHB 
TEGE 
OGT, 
OVP, 
2 

BAR 
S 

TS 

uW  ( 
S  ( 
LPT( 
( 
UB2( 
LLS( 
CT^ 
UNDL 
kCON 
RECE 
AND 
SCk( 
PfL, 
)  .N 
.Lfc. 
GATE 


LP  (7X,  * 

R*2  (A-T, 

NUVP,  NI 


20000) 
4000) 
2000) 
^000) 
2000) 
2000) 


NECTKJN  B 
DING  TO  G 
THESE  INT 
20),  PRn( 
G.  V,H 
E.  A)  -^ET 
0)  RETUR 
S  SUCCESS 


,  *) 

V-Zt  i) 

NOGT 
0  3 
M 

JX 

CLCF(2 
RWPT( 
NS   ( 
CLLT( 
XS   ( 
SUP  ( 
NCT2 
NCT6 
NFI  XJ 
ETWFEN 
ATE  I  t 
FRCONNE 
20) 

URN 

N 

IVF  TO 


000  0) 
4000) 
4000) 
^000) 

2000) 


GATE  I 
AN  NOT 
CTIONS 


NIPT 

U4 

N 

NCX 

COL  (20000) 

RWLT(  4000) 


NOPT 

JAO 
LTS 

^■wCF(  2000  0) 
YS   (  4000) 


CLVT(  400U),  CLSP{  4000) 

OF   (  2000),  XHAT(  2000) 

S    (  2000),  SUBK  2000) 
TP   (  2000) 

NCT3         ,  NCT4 

NGETYS       ,  NNEwX 
KWLS(  4000) 
Ar40  GATt  J  IS  SET  TJ  ONE 

FEED  THOSE  GATE  SUCCESSIVE 
ARE  SET  TO  ZERO 


/I, 2, 3, 4, 5, 6,  7,8 


GATE  J  AND  PUT  THEIR  GATE  NO. IN  SC^ 


=SUrt2( JAD) 

EGl  .GT.  LNCl) 

BEGl) 

J =1, NOGT 


GO  TO  2  0 


GO    TO    10 
.    0)    GO    TO    10 
GO    TO    10 


IF     (I     .EO.    J) 
SUB  =  $A( I  ,J» 

IF  (S(SUB)  .LE 
DO  8  JK=l,ENDl 
IF     (SCR( JK).EU.J) 

8  CONTINUE 

eN01=ENDl+l 
SCR(END1)=SUB2<SUBI 
10    CONTINUE 

BEG1=BEGI+1 
GO    TO       5 
C  FIND    THE    GATES    PRECEDING 

20    rt£G2=l 
END2=1 

PR0(1)=SUB1( JADI 
25    IF     (BEG2    .GT.    END2I 
J=PRD( BEG2) 
DO    iO       I=ltNOGT 
IF     (  I     .EQ.     J)     GO    TO 
SUB=$A( I ,J) 
IF     (S(SUB)     .LE.    0) 
DO    Zti    JK=1,END2 
IF     IPRD(JK).EQ.I )     GO    TO    30 
28  CONTINUE 

END2=END2+l 
''f<D(ENn2)  =  SUBUSUB) 
30    CONTINUE 

BEG2=fiEG2*l 
GO  TO  2b 
^0  CONTINUE 

DO  60  I=l,ENDl 
DO  50  J=1,END2 
SUB  =  $A(SCK(  F  )  ,PHD( J) ) 


TO  GATE  I  AND  PUT  THEIR  GATE  NC.  IN 


PRO 


GO  TO  40 


30 


GO  TO  30 


50 

oO 


100 


110 


IF  (S(SU6)  , 

JAD=SUB 

KS=-l 

XS{JX)=-JAD 

CALL  DTPSNS 

CONTINUE 

CONTINUE 

RETURN 

INFFASIBLE 

RETURN  1 

FEASIBLE 

RETURN  2 

END 


LT.  OJ  GO  TU  5  0 


(ZX,  tlOO,  tllOJ 


SUBROUTINE  PRESET  CZX,  ♦) 
IMPLICIT  INTEGER*2  (A-T,V-Z,$) 
INTEGER*^  NOGT,  NUVR,  NIPT,  NOPT 


COMMON 
1 

COMMON 
1 
2 
3 
4 
5 
6 
7 


NOVR 

U2 

ZSAR 

KS 

JTS 

ROW  (20000) 

PS   (  4000) 

CLPT(  2000) 

X    (  2000) 

SUB2(  2000) 


NOGT 

U3 

M 

JX 

CLCF(20000) 

RWPTC  4000) 


NS  ( 
CLLT( 
XS  < 
SUP  ( 


4000) 
2000) 
2000) 
2000) 


NIPT 

U4 

N 

NCX 

COL  (20000) 

RWLT(  4000) 


CLVTI 
OF  ( 
S  ( 
TP   ( 


4000) 
2000) 
2000) 
2000) 


NJPT 

JAD 
LTS 

RWCF(20000) 
YS  (  4000) 
CLSP(  4000) 
XHAT(  2000) 
SUBK  2000) 


COMMON 


CLLS( 

NCT3 
NUNDLX 
1  =  1, M 


2000) 


NCT2 
NCT6 
NFIXJ 


on  10 

YS{  I  )=0 

10  CONTINUE 

DO  13  1=1, N 
S( I )=0 
X( I »=0 

xs(  n  =  o 

XHAT  ( n  =  0 
13  CONTINUE 

JX=1 

NCX=0 

U2  =  0 

U3=0 

U4=  0. 

JAD  =  0 

KS  =  1 

CALL  GETYS 

CALL  SINTL  (ZX, 

CALL  GTPSNS 

LTS  =  0 

RETURN 
20  RETUKN  1 

END 


L20\ 


NCT3        ,    NCT4 
NGETYS       ,    NNEWX 
RWLS(  4000) 


SU 
IM 
IN 
CO 


hPUU 
PLIC 
TtGt 
MM  ON 


COMMON 


COMMON 


1000 


1009 


10 


CI 
DA 
IF 
IF 
IF 
I^ 
PR 
Fu 
PR 
FfJ 
DO 
SU 
CN 
IF 
BU 
GO 
BU 


MENS 
TA 

UE 
WE 
WE 

INT 
R^'AT 
INT 
RMAT 
90 
t^=$L 
T=l 

(OU 
F(CN 

TO 
FICiM 


tin:  PRNT 
IT  INTPGER 
R*4  NOGT, 

NOVR 

U2 

ZbAR 

KS 

jrs 

ROvJ  (2 
PS   ( 
CLPT( 
X     ( 
SU8  2( 
CLLS( 
NCT5 
NUNDLX 
ION  OUT (20 
CP         / 
WANT  TO  PR 
WANT  TO  PR 
WANT  TO  PR 
T  .NE.  1) 
lUOO,  (K, 
( //20X, 'K 
1009,  (J, 
(•♦• ,  41X 
1=1, NOGT 
(  I  ) 


(  OUT, 
*2  (A-T, 

NOVR,  NI 


0000) 
4000) 
2000) 
2000) 
2000) 
2000) 


00) ,BUF( 
•CP* 

INT  LAMO 
INT  BETA 
INT  PART 
GO  TO  10 
K=l ,NOVR 
=•,413) 

J=1,NUG 
,  '0=', 


ST,  NR,  PTS) 

V-Z,t) 

PT,  NOPT 

NOGT 

0  3 

y 

JX 

CLCFI 20000) 

RWPT(  4000) 


4000) 
2000) 
2000) 
2000) 


NS   ( 

CLLT( 

XS   { 

SUP  ( 

NCT2 

NCTt> 

NPIXJ 

/ 

A, ALPHA, W, 
,  GAMA,P, 


THEN 
THEM 


NIPT 

U4 

N 

NCX 

COL  (20000) 

KWLT(  4000) 


4000  ) 
2000) 
2000) 
2000) 


CLVT( 
OF  ( 
S  I 
TP   ( 

NCT3 
NGETYS 
PWLS(  4000) 


ST=1 

^18=1 


I  AL 

0 

) 


SOL.  THEN  PTS=l 


T  ) 
1213/) 


T( SOB) 
T)=l 
20 
T)=0 


.LE.  0)  GU    TO  IJ 


NOPT 

JAD 
LTS 

DrtCF(2a0u0) 

YS   (  4000) 

CLSP(  4000) 

XHAT(  2000) 

SUBK  2000) 

NCT4 

NNEWX 


20 


30 
1010 


t>0 
60 

1020 
90 


102!j 


9^ 
95 

I02d 

98 

100 


120 


1030 

130 
1040 

140 


0)  GO  TO  30 


I,  BUFdl 


I* 
13, 


(BUF(AA»,  AA=2,CNT) 


0)  GO  TO  50 


I  , (8UF( AA),AA=l,CNT) 
,36X,'A(  '  ,  12,  •,  J)  •, 1213) 


CONTINUE 

DO  40  K=l,NOVR 

SUti=SW(  I  ,K) 

CNT=CNT<-1 

IF  (OUT(SUBI  .LE 

BUF(CNT)=1 

GO  TO  40 

BUF(CNT)=0 

CONTINUE 

PKINT  1010, 

FORMAT  (5X,  'H   •  ,  I2,«  »• 

CNT=0 

00  60  J=I,NOGT 

CNT  =  CNT-H 

IF  (  I  .EO.  J)  GO  TO  50 

SUB=$A( I ,J) 

IF  (0UT(SU8)  .LE. 

HUF(CNT»=1 

GO  TO  60 

BUF(CNT)=0 

CONTINUE 

PRINT  1020, 

FORMAT  (•+• 

CONTINUE 

IF  (NCT5  .NE.  CP)GU  TO  100 

PRINT  1025,  (SUB,SUB=l,NdVR» 

FORMAT  ( //10X,'K=»,30I3) 

DO  98  I=1,N0GT 

DO  95  J=1,N0VR 

SU'3=$V(I  ,JI 

IF  (OUT(SUB) .LE.O)  GO  TO  94 

bUF( J|=l 

GO  TO  95 

BUF( J)=0 

CONTINUE 

PRINT  1028, I , (6UF( AA),AA=l,NUVP ) 

FORMAT  (5X, 'V( •, 12, • ,K)' , 3013) 

CONTINUE 

IF  (NP  .NE.  1) 

nDR  =  0 

NI PT=2*«N0VR 

DO  190   I=l,NOGT 

DO  180   J=1,N0GT 

CNT  =  0 

IF  (I  .EU.  J)  GO 

0DR=0DR*-1 

IF  (ODR  .GT.  3)  UDR^l 

DO  120  K=l,  NIPT 

SUB  =  $B( I  ,J,K) 

CNT=CNT+I 

6UF(CNT)=aUT(SUB) 

CONTINUE 

IF  (OCR  .EQ.  3)  GO 

IF  (ODk  ,EQ.  2)  GO 

PRINT  1030,  I,  J, 

FORMAT  (5X,«B(« 

GO  TO  180 

PRINT  1040,  I, J,  (Bl)F(  SUB)  ,S'JB=1,CNT) 

FORMAT  (•  +  • ,48X,«B(  •,  12, •,  •,  12,  •  ): '  ,1612) 


GO  TO  400 


TO  180 


TO  140 
TO  130 
(BUF(S'J8),SU6  =  1,CNT) 
•,1612) 


GO  TO 
PRINT 


180 
1050, 


I,  J,  (BUF(SUB) ,SUB=1,CNT) 


1050 
IHO 
190 


220 


L060 

230 
1070 

240 
1080 
280 
2^0 
293 


32U 


los^o 

330 
1100 

340 

1110 

3  90 

400 

1120 

1130 
500 


.LT.  2)  GO  TO  293 


PURMAT  C  +  '.giX.'BC 

CONTINUE 

CONTINUE 

IF  (NCT4 

ODK  =  0 

on  290  I=1,N0GT 

00  280  J=1,N0GT 

IF  ( I  .EQ.  J )  GO  TO  280 

CNT  =  0 

U0R=0()K+1 

IF  (ODR  .GT.  i)    0DP=1 

DO  220  K=1,NIPT 

SJ3=$G( I tJfK) 

IF  (SUB  .GT.  N)  GO  TO  293 

CNT=CNT+1 

aUF(CNTJ=nUT( SUB» 

CONTINUE 

IF  (OUR  .EU.  3)  GU  TO   240 

IF  (ODR  .EQ.  2)  GO  TO   230 

PRINT  1060,  I , J, (BUF(AA) , AA=l,CNT) 


12,  •  ):  •,  1612) 


FORMAT  (5X, 
GO  TO  280 
PRINT  1070, 


G( 


12,'  )  :• ,1612) 


I,     J,(8UF(AA),     AA=1,CNT) 


FOKr^AT     (  •♦■•  ,48X,«G(  •  ,12,  •  , 
GO    TO    280 

I,     J,     (BUFCAAJ 
,91X, 'Gi • 


•  ,  I2,»  ):•  ,  1612) 


PRINT     1080, 


FORMAT     ( 

CONT INUE 

CONTINUE 

CONTINUE 

CDK=0 

DO     390     I=1,N0GT 

00R=0DR+1 

IF     (ODR    .GT.     3)     0DR=1 

CNT=0 

DO    320    K=1,NIPT 

SUB=$P( I ,K) 

CNT=CNT+1 

3UF(CNT)=0UT(SUB) 

CONTINUE 

IF        (ODR    .EQ.     3)       GO    TO 

IF        (UOR     .EQ,     2)        GO    TO 

PKINT     1090, I ,     ( BJF( AA) , 

FORMAT     (8X,«P(  •  ,  12, •  ):• , 1612) 

GO    TO    390 

PRINT     1100,1,     (BUF(AA),    AA=1,CNT) 

FORMAT     (•*-»,51X,«P(«,I2,M:»,loI2) 

GO    TO    390 

PRINT     1110,1,     {BUF(AA),     AA=1,CNT) 

FORMAT     (•+',94X,'P(',I2,'):',16I2) 

CONTINUE 

IF     (     PTS    .NF.     1)    G(]    TO    50J 

PRINT     1120,     (X(AA),     AA=1,N) 

FORMAT     (3X,     •     X«/(5X,     2016)) 

PRINT     1130, (XS( AA) ,    AA=1,N) 

FORMAT     (3X,     'XS'/(5X,     20Io)) 

RETURN 

END 


AA=1,CNT) 

•  ,1612) 


340 
330 
AA=1 ,CNT) 


SUBKOUTINE    J12Q1 


GENERATE  BASIC  IN 
IMPLICIT  INTEGER* 
INTEGER*^  NQGT,  N 


EJUALIT 
2  (A-T, 
OVR,  NI 


NOR  NETWORKS. 


COMMON 
I      t 

COMMON 
1 
2 
3 
4 
5 
6 
7 

COMMON 
I 
2 

01  MENS  ION 


NOVR 
U2 

ZBAR 
KS 
JTS 

ROW  (20 
PS  C  4 
CLPT(  2 
X  (20 
SUB2(  2 
CLLS(  2 
NCT5 
NUNDLX 
TEMPdO 


EQUIVALENCE  (TEMP 
102) ), (CMT3,X(203) 

DATA  E/'E'/ 
DATA  B1,82.G1,G2 
DATA   W,V,D,C,B 
DATA     CP 
DATA   NR,ND,AO,GT 


000) 
000) 
000) 
00), 
000) 
000) 


0) ♦VAL( 
(1  )  ,X(1 
)  .( SQ,X 


/•CP' 

.ST, IN, 
AD« , 'GT 


NOPT 

JAO 
LTS 

RWCF(20000) 
YS  (  4000) 
CLSP(  40001 
XHAT(  2000) 
SUBK  2000) 


1011 


DATA 
PRINT 
FORMAT 
DO  70 


CL  ,SI 
1011 
(  •  1»  ) 

1=1, NQGT 


lES  FOR 
V-Z,$) 
PT,  NOPT 

NOGT 

U3 

M 

JX 

CLCF(20000) 

RWPT(  4000) 

NS   (  4000) 

CLLT(  20001 
XS   (  2000), 

SUP  (  200J) 

NCT2 

NCT6 

NFIXJ 
100), INPT(270) 

)  )  ,(VAL (1) ,X(101) ) ,(CMTl,X(20l) ),(CMT2 
(204) ) , (INPT(l) ,XS( 1 )) 


G2' 


GT,OU,FI ,FO,BLK 

•,'ST«,'IN«,«GT«,'UU»,'FI«,«Fn«,« 
'  / 


NIPT 
U4 

N 

NCX 

COL  (20000) 
RWLT(  4000) 
CLVT(  4000) 
OF  (  2000) 
S  (  2000), 
TP  (  2000) 
NCT3 
NGETYS 
RWLS(  4000) 


NCT4 

NNEWX 


X{2 


73 
71 


72 
70 


DO  71  J=1,N0VR 
IF  (NCT5  .NE.  CP)  GO  TO  73 
OF  (  $  V  (  I  ,  J »  )  =  1 
CONTINUE 
OF($W(  I, J)  )  =  l 
DO  72  JJ=1,N0GT 
IF  ( JJ.EO.I  )  GO  TO  72 
OF ($A(  I, JJ))=1 
CONTINUE 
CONTINUE 
U=N0GT*NOVR+2 
0NE=1 
CMT2=RLK 
CMT3=GT 
CCCCCCCBASIC  INEOULITIESCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO   100  1=1, NOGT 
S0  =  1 

on   10  AA=1,100 
TEMP( AA)=0 
VAL(AA)=0 
10  TP(AA)=BLK 
N=l 

TEMP(N)=0 
VAL(N)=-1 
TP(N)=C 
N=N+1 

TEMP(N)=$W( I, ONE) 
VAL(N) =1 
TP(N)=W 

IF  (NCT5  .NE.  CP  )  GO  TO  15 
N=N+1 
TEMP(N)=$V( I ,ONEi 


VAL(N)=1 

TP(N)=V 
15    CONTINUE 

00    ^0         J=l,NOGT 

IF     (J     .FQ.     I )    GO    TO    20 

JJ=I 

IF     (     I     .GT.     J)        JJ=JJ-l 

N=N+1 

TFMP(N )^$B( J, I, ONE) 

VAL(N) =1 

TP(N)=0 
?0    CONTINUE 

N=N+1 

TFMP(N)=$P( I ,ONE) 

VAL(N)=U 

TP(N)=D 

CALL     SPUNCH 

DO    30    AA=1,N 
30    VAL( AA)=-VAL(AA) 

VAL( 1)=U 

CALL    5;PUNCH 
LOG    CONTINUE 

DO       200    I=l,NOGT 

DO     190       J=1,N0GT 

IF     ( I     .EQ.     J)     GO    TO    190 

C^T3=R1 

sg=i«io+j 

N=4 

DO     110    AA=1,100 

T£MP< AA)=0 

VAL( AA)=0 
110    TP{AA)=bLK 

TEMPd  )  =  0 

VAL( 1)=1 

TP  (1)  =  C 

TEMP(2 )=$P( I ,ONE) 

VAL(2) =-1 

TP(2)=0 

TFMP(3)=$A(  [  , J) 

VAL( 3) =-l 

TP(3)=C 

Tt^MPI^fMSdl  I  ,  J,ONE) 

VAL(^)=1 

TP(4)=D 

CALL     SPUNCH 

on     120    AA=1,100 

TEMPI AA)=0 

VAL(AA)=0 
120    TP(AA)=t^LK 

N=3 

CMT3=B2 

TF(MP{  1  )  =  tP(  I  ,ONt) 

VAL{ 1)=1 

TP(1)=D 

TEMP(2)=$A( I , J» 

VAL     (2)=1 

TP(2)=C 

TEMP(3)=$a( I , J,ON£) 

VAL(3) =-2 

TP( 3)=D 

CALL     SPUNCH 


190  CONTINUE 
200  CONTINUE 

RETURN 

END 


SUBROU 

C      GENERA 

C    ■  THE  CA 

C      ADDITI 

C      THESE 

C      MAKE  A 

C      THESIS 

IMPLIC 

INTEGE 

COMMON 

1 

COMMON 
1 
2 
3 

5 

6 

7 
COMMON 

1 

2 

1)1  MEN*: 
EQUIVA 

102)), ( 
DATA  E 
DATA 
DATA  R 
DATA 
DATA 

I 
DATA 
DATA    L 
U=Nt,GT 
N0P1=N 
0NE=l 
TW0  =  2 
THREE= 
CCCC       EACH    G 
CCCC       EXTERN 
CMTl=N 
CMT2=B 
CMT3=I 
DU    100 
SQ=I 
DO    10 
TEMP{A 
VAKAA 
10    TP{AA) 
N=l 

TEMP(N 
VAL(N) 
DO    20 
N=N+1 
TEMPiN 


LITIES  FOR  NOR  NET 
UN  THE  TUP  IMPLEM 
OSE  EFFECTIVENESS 
THE  PROGRAM  ON  MAY 
GTINEU  UNIFORM  ( SE 


/•CP« 
2 

GT.ST, IN.GT 
,• AO' , 'GT' , 
/•CL* ,'SI» 


Z,t) 

,  NOP- 

OGT 

3 

X 

LCF(2 
WPT( 
S   ( 
LLT< 

(  2 
UP  ( 
CT2 
CT6 
FIXJ 
0)  ,  IN 
,  (  VAL 
04)  )  , 


0000) 
4000) 
4000) 
2000) 
000)  , 
2000) 


PT(270) 
(1)  ,X(1 
(  iNPTd 


NIPT 
U4 
N 

NCX 
COL 
RWLT 
CLVT 
OF 
S    ( 
TP 

NCT3 
NGET 
RWLS 

01)  ) , ( 
)  ,XS(  1 


WORKS. 

ENT       CHANGES     IN    THE 

HAVE    NOT    BEEN    TESTED. 

1,  1974  IN  ORDER  TO 
E  HOHULIN  MASTER 


NOPT 


(20000) 
i  4000) 
(  4000  ) 
(  2000) 
2000) , 
(  2000) 

YS 

(  4000) 


CMT1,X(201  )  ),(CMT2,X(2 
)) 


TINE  U12Q2 

TE  ADDITIONAL  INEQUA 

RDS  MARKED  WITH  BLUE 

ONAL  INEQUALITIES  WH 

CHANGES  WERE  PUT  IN 

LL  FOUR  VERSIONS  OF 

-  APPENDIX  B). 

IT  INTEGER*2  (A-T,V- 

R*4  NOGT»  NOVR,  NIPT 

NUVR         ,    N 

U2  ,    U 

ZBAR         ,    M 

KS  ♦    J 

JTS         ,    C 

ROW  (20000),    R 

PS   (  4000),    N 

CLPT(  2000),    C 

X    (2000),    XS 

SU62(  2000),    S 

CLLS(  2000),    N 

NCT5         ,    N 

NUNOLX      ,    N 

ION  TEMH(IOO) ,VAL( 10 

LENCE  (TEMP(l) ,X( 1 ) ) 

CMT3,X(203)  )  ,( S3,X(2 

/•  E'/ 

CP 
l,82,Gl,G 
M,V,D,C,B 
NR ,N0, AD, 
•Nk' , 'NO' 
CL  ,SI 
V/ 'LV'/ 
♦NOVR+2 
OPT+1 


ATE  EXCEPT  THE  OUTPUT  GATES  HAS  AT  LEAST  ONE  INPUT 
AL  VARIABLES  OR  AT  LEAST  TxO  INPUTS  FROM  THE  OTHER 
R 

LK 
N 
I=NOP1,NOGT 

AA=1,100 

A)=0 

)=0 

=  BLK 

)=0 
=-2 

K=1,N0VR 

)=$W(I ,K) 


JAD 
LTS 

RWCF(20000) 
YS  (  4000) 
CLSP(  4000) 
XHAT(  2000) 
SUtlK  2000) 

NCT4 
NNEWX 


.•G2 


,nu,FI,FU,aLK 
•ST«,  •  IN*  ,  'GT' , 'OU 
/ 


FROM  THE 
GATES  CC; 


VAHN)=2 

IF  (NCT5  .NE.  CP)  GO  TO  20 

VAL (NJ=l 

N=N-H 

TEMP(N)=$\/(  I  ,K) 

VAL(N)=1 
^0  CONTINUE 

DC)  3U  J  =  1,N0GT 

IH  ( I.EO.JIGO  TO  30 

N=N*-1 

TE^1P(N)=$A(  J, I  ) 

VAL(N)=l 
30  CONTINUE 

CALL  SPUNCH 
100  CONTINUE 
CCCCCCTRAINGULAR  INEQUAL I T I ESCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCC 
C      IJ  =  2  FOR  SINGLE-OUTPUT  NETWORKS;   IJ  =  1  FOR  MULT  I  PL t -OUTPUT 
C      NETWORKS.   PREVIOUSLY  ONLY  IJ= I  WAS  USEO. 

IJ  =  2 

IF  (  NOPT  .GT.  I  )  IJ  =  I 

CMT3=ST 

SQ  =  0 

no  200  I=IJ,NOGT 

00  190  J=IJ,NOGT 

IF  (J  .EQ.  1  )  GO  TO  I'^O 

DO  180  K=l,NOGT 

IF  (  K.EQ.I  .OR.  K.EO.J)  GO  TO  180 

DO  110  AA=l,100 

TEMPI AA)=0 
110  VAL(AA)=0 

Sy=SQ+l 

IF  (SQ  .GE.  100  J  sg=i 

TEMPd  J=0 

VAL(1)=2 

TEMP(2  )  =  $A(  I, J) 

VAL(2) =-1 

TFMP(3)=$A( I ,K) 

VAL(3)=-1 

TEMP(^»=$At J,K) 

VAL(4I=-1 

N=4 

DO  120  L=1,N0GT 

IF  (L  .EQ.  KJ  GO  TO  120 

IF  (  L  .Eg.  J)  GO  TO  120 

IF  (  L  .EQ.  I  )  GO  TO  120 

N=N-H 

TEMP(N)=$A(  J,L) 

VAL(N)=1 
120  CfJNTINUE 

CALL  SPUNCH 
180  CONTINUE 
190  CONTINUE 
200  CONTINUE 

00  300  J=IJ,NOGT 

00  290  K=L,NUGT 

DO  280  L=l,N0VR 

IF  (J  ,EQ.  K)  Gf)  TO  290 

SQ=S0+1 

IF  (SQ  .GF.  lOOJ  SO=l 

DO  210  AA=1, 100 

TEMP(AA)=0 


210  VAL<AA)=0 
TEMPd  )=0 

VAL( l)=2 

TEMP(2)=$W( J,L) 

VAL(2)=-1 

TEMP(3)  =  $v^(K,LI 

VAL(3l=-l 

TEMP(4I=$A(J,K| 

VAL(^J=-1 

N=4 

DO  220  TT=1,N0GT 

IF  (TT.EO.K)  GO  TO  220 

IF  (  TT  .EQ,  J)  GO  TO  220 

TEMP(N)=$A< J,TT| 

VAL(N) =1 
220  CONTINUE 

CALL  SPUNCH 

IF  (NCT5  .NE.CP  )  GO  TO  280 

TFMP(2)=SV( J,L) 

TEMP(3)=$V(K,L) 

CALL  SPUNCH 
280  CONTINUE 
290  CONTINUE 
300  CONTINUE 

CCCC   ALL  GATES  EXCEPT  THE  OUTPUT  GATES  CONNECT  TO  AT  LEAST  ONE  GAT^  CCC 
CMT3=0U 

sy=o 

no  400  I=N0P1,N0GT 

SO=SQ+I 

DO  310  AA=l,100 

TEMP(AA)=0 
310  VAL( AA)=0 

N=l 

TEMPd  )  =  0 

VAL(1)=-1 

IF(  I.EQ.U  VAL(  1)=0 

DO  320  K=l,NOGT 

IF  (I.tO.K)  GO  TO  320 

N=N*l 

TEMP(N)=$A( I,K) 

VAL(N)=l 
320  CONTINUE 

CALL  <;  PUNCH 
400  CONTINUE 
CCCC   GATES  WHICH  APE  CONNECTED  TO  ALL  OF  THE  OUTPUT  GATES  ARE  NOT 
CCCC   CONNECTED  TO  ANY  OTHER  GATES  CCC 

SQ^O 

CMT3=0U 

DO  600  I=NOPl,NOGT 

sg=so+i 

DO  510  AA=1, 100 
TEMP(AA»=0 
510  VAL(AA»=0 
TEMP(1)=0 
VAL(l) =NOPT*U 
DO  515  J=2,N0P1 
JMl=J-l 

TEMP(J)=$A(I,JM1) 

VAL(  J)=-IJ 
N=J 


515    CONTINUf 

DO  520  K=N0P1,N0GT 

IF  (K.FQ.I )  GO  TG  520 

N=N+l 

TFMP(N)=SA( I ,K) 

VAL(N»  =-1 
520    CONTINUei 

CALL  SPUNCH 
600  CONTINUE 

CALL  CLEAR 

INPTd  )  =  E 

IMPT(6)=-99 

N9  =  6 

CALL  INTLP(lMB) 

RETURN 

END 


10 
1000 


SURROU 
IMPLIC 

INTEGE 
COMMON 

1 

COMMON 

1 

2 

3 

4 

5 

6 

7 
COMMON 

I 

2 

IF(  RW 
K=RWPT 
J=  RWL 
IF(COL 
PR  INT 
FORMAT 


1 


20 


1010 


1 


30 


GLJT(i    3 
K    =    K  + 
IZ    =    - 
PRINT 
FORMAT 

RETURN 

ENO 


TINE  ROWP 
IT  INTEGE 
R*4  NOGT, 
NOVR 
U2 

ZBAR 
KS 
JTS 
ROW  ( 
PS   ( 
CLPT( 
X    ( 
SU82( 
CLLS( 
NCT5 
NUNDL 
LT(I)  .EO 
(I  ) 

T(  n+K-1 
(RWPTd  )  ) 
1000,     I, 
(     15,     6X, 
(16X,        9{ 
0 
1 

RWCF(K-II 

1010,  I, 

(  I5,I7,« 

(16X,   9( 


T{  I) 

R*2  (A-T,V-Z,$) 

NOVR,  NIPT,  NOPT 
NOGT 
U3 
M 

JX 

CLCF(20000) 
RWPT(  4000) 


20000J 
4000) 
2000) 
2000) 
2000) 
2  000) 


0)  GOTO 


N*;  ( 
CLLTI 
XS  { 
SUP  i 
NCT2 
NCT6 
NFIXJ 
30 


4000) 
«^000) 
2000) 
2000) 


10,20,10 


NIPT 

U4 

N 

NCX 

COL  (20000) 

RWLTI  4000) 

CLVT(  4000) 

OF   (  2000) 

S    (  2000) 

TP   (  2000) 

NCT3 

NGETYS 

KWLSi  4000) 


)  10,20,10 

(RrtCF{  L)  ,  COL(L)  ,  L  =  K,J) 

•0  LE  ',  9(  •♦•  ,  I  i,  •  X( 
•  +  ',  13,  •  X(«  ,  14,  MM  ) 


,14, 


)/ 


IZ,  (RWCF(L), 
LE  •  ,   9(  •♦■•, 

X(»  ,  14 


COL(L) , 
13,  • 


L  = 
X( 


)•  ) 


K,  J) 

•,  i4,')' 
)  ) 


)/ 


NOPT 

JAD 
LTS 

RWCF(20000) 

YS   (  4000) 

CLSP(  4000) 

XHAT(  2000) 

SU81(  2000) 

NCT4 
NNFWX 


SUBROUTINE   RVPSNS 

IMPLICIT  INTEGER*2  (A-T,V-Z,$) 

INTEGER*4  NOGT,  NOVR,  NIPT,  NOPT 


COMMON 

NOVR 

,    NOGT 

( 

,          NIPT 

» 

NOPT 

1 

U2 

f          'J3 

1 

»    U4 

COMMON 

ZBAR 

M 

t 

r     N 

f 

JAO 

1 

KS 

,          JX 

1 

NCX 

» 

LTS 

2      , 

JTS 

CLCFI20000), 

,          COL  (20000), 

RWCF(20000) 

3      , 

ROW 

120000), 

KWPT( 

4000) , 

RWLT( 

4000), 

YS   (  -tOOO) 

4       , 

OS 

I  4000), 

NS   ( 

4000) , 

,          CLVT( 

4000  )  , 

CLSP(  4000) 

5  , 

6 

7 

COMMON 
1 
2 
BEG    =    CL 
ENO=CLLT 
IF     ((END 
IF    <KS.G 
C  ERASE    X( 

NCX=NCX- 
DO    30     1= 
SUB=ROW( 
IF     (CLCF 
10    PS(SUB)= 
GO    TO    30 
20    NS(SUB)= 
IF     INSIS 
IF     (CLSP 
CLSPINCT 
NCT3=SU8 
CLSP(NCT 
30    CONTINUE 
RETURN 
C  ERASE 

5  0    DO    80     1= 
SUB=ROW( 
IF     (CLCF 
6U    NS(SUB)= 
IF     (NS(S 
IF     (CLSP 
CLSPINCT 
NCT3=SUB 
CLSPINCT 
GO    TU    80 
70    PS( SUB  )  = 
80    CONTINUE 
RETURN 
END 


2000) 
2000  » 
2000) 
2000) 


KS=-l 


CLPT( 
X    I 

SUB2( 

CLLS( 

NCT5 

NUNOLX 
PT(JAD«-l) 
(JAD+D-l+BEG 
-BEG)  .LT.  0) 
T.O)  GO  TO  50 
L)=l 
OFIJAD) 
BEG, END 
I) 

(I))  10,30,20 
PS(SUB)-CLCF( I) 


NS(SUB)-CLCF( I) 
UB).GE.  0)  GO  TO 
(SUB)  .NE.  0)  GO 
3)=SUB 


3)=SU8 

X(L  )  =  0 
BEG, END 
I) 

(I))  60,80,70 
NS(SUB)«-CLCF(  I) 
UB)  .GE.  0)GO  TO 
(SUB)  .NE.  0)  GO 
3)=SUB 

3)=SUB 

PS( SU8)+CLCF( I) 


CLLT( 
XS   ( 
SUP  ( 
NCT2 
NCT6 
NFIXJ 


2000) 
2000) 
2000) 


RETURN 


30 
TO 


30 


KS=  1 


SO 
TO 


80 


OF   (  2000), 

XHAK 

2000) 

S    I  2000), 

SUBK 

2000) 

TP   (  2000) 

NCT3         , 

NCT4 

NGETYS      , 

NNEMX 

RWLS(  4000) 

SUBFOUT 
IMPLICI 

INTEGER 
COMMON 

I 
COMMON 

1 

2 

3 

5 

6 

7 


COMMON 

1 

2 
CLPT(l) 
S( 1)=1 
NN=N+1 


INE   RWTUCL 
T  INTEGER*2  (A 
*4  NOGT,  NOVR, 
NOVR 
U2 

Z8AR 
KS 
JTS 

ROW  (20000) 
PS  (  4000) 
CLPT(  2000) 
X  (  2000) 
SUB2(  2000) 
CLLS(  2000) 
NCT5 
NUNDLX 
=  1 


T,V-Z,$) 

NIPT,  NOPT 
NOGT 
U3 
M 

JX 

CLCF(20000) 
RWPT(  4000) 
NS  (  4000) 
CLLT(  20001 
XS  (  ZOOO) 
SUP  (  2000) 
NCT2 
NCT6 
NFIXJ 


NIPT 

U4 

N 

NCX 

COL  (20000) 

RWLT(  4000) 

CLVT(  4000) 

OF   (  2000) 

S    (  2J00) 

TP   (  2000) 

NCT3 

NGETYS 

RWLS(  4000) 


NOPT 

JAD 
LTS 

RWCF{20000) 

YS   I  4000) 

CLSP(  4000) 

XHAT(  2000) 

SUBK  2000) 

NCT4 

NNEWX 


EU.  0)  GO  TO  50 


no  10   L=2,NN 

CLPKL  )  =  CLPT(L-l)+CLLT(L-l  I 

S(L)=CLPT(L) 
10  CONTINUE 

00  50   K=l,M 

IF  (KrtLT(K) 

REG=RWPT(K) 

EN0=8tG+RWLT(K)-I 

00  40   J=8EG,END 

SUB=COL( J)+l 

CLCF(S(SUB) J=PWCF(J) 

RO»(/(S(  SUB)  )=K 
40  S(  SUB»=S(SUB)-H 
50  CONTINUE 

RETURN 

END 


SJdROU 
IMPLIC 
INTEGE 
COMMON 
1 

COMMON 
I 
2 
3 
4 
5 
6 
7 

COMMON 
I 
2 
X(  I)=XS 
X(  I  )  =  -XS 
S ( I )=  1 
S( I )=  0 
S( I )=-l 
01  MENS 
DATA 
DO  200 
DO   20 
UA{  I  I  , 
205  CONTIN 
200  CONTIN 
K=0 

DO  290 
READ  2 
210  FORMAT 
AB=5-A 
READ 
220  FORMAT 
^0  260 
K=K+l 
XS(K)= 
X(K)=l 
2oO  CONTIN 
290  CONTIN 
K=K-»-l 
00   30 


TINE  S 
IT  INT 
R*4  NO 
NO 
U2 
Zfl 
KS 
JT 
RU 
PS 
CL 
X 

SU 
CL 
NC 
NU 
IF 
IF 
X(  I 
X(  I 
X(  I 
UA 


INTL  (ZX,  *J 
EGER*2  (A-T, 
GT,  NOVR,  NI 
VR 


AR 

S 

W  (20000) 
(    4000) 

PT(  2000) 
2000) 
2000) 
2000) 


I  ) 
I  ) 
IF 
IF 
IF 
ION 
BLK  / 
I  1=1, 
5  JJ  = 
JJ)=BL 
UE 
UE 


( 
82{ 
LS( 
T5 

NULX 
X( I )=1 
X(  I  )=0 

)  =  1 

)=UNASSIGNED 

)=0 

(4,20),  IT(2 

•   •    / 

4 

1,20 

K 


V-Z,$) 

PT,  NOPT 
NOGT 
1)3 
M 

JX 

CLCF(20000) 
RWPTI  4000) 
NS  (  4000) 
t:LLT(  2000) 
XS  (  2000) 
SUP  (  2000) 
NCT2 
NCT6 
NFIXJ 


NIPT 

U4 

N 

NCX 

COL  (20000) 

R«vLT(  4000) 

CLVT(  4J00 ) 

OP   (  2000) 

S    (  2000) 

TP   (  2000) 

NCT3 

NGFTYS 

RWLS(  4000) 


NOPT 

JAD 
LTS 

RWCF(20000) 

YS   (  4000) 

CLSP(  4000) 

XHAT(  2000) 

SJBK  2000) 

NCT4 
NNEWX 


0) 


AA=1,N0PT 
10,  L,  (  IT{  I  ),  1  =  1, NIPT) 

(  13,  2X,  161  I) 
A 
220,  (UA{AB,I  ),  1  =  1,20) 

(20A4) 

1=1, NIPT 

-$P(AA,I ) 
-IT(  I  ) 

UE 
UE 

0  EE=1,2 


310 


340 
330 
320 
300 
10 

20 


30 

40 
50 

60 
70 
80 
90 
92 

100 


llO 
120 

86 


GU  TO  60 


PRINT  310 

FORMAT  i*iW///////////////////) 

DO  320  11=1,4 

on   330   JJ=1,I0 

PRINT  340,  (UA( II.KI ),KI=1,20) 

FORMAT  130X,  20A4I 

CONTINUE 

CONTINUE 

CONTINUE 

J=K+11 

READ  20,  (XS( n ,  X( I ),  I=K,J) 

FORMAT  (12(15,11)1 

IF  (XS(K)  .GT.  N)  RETURN  1 

IF  i  XS( J).EQ.O»  GO  TO  30 

K=J  +  1 

GO  TO  10 

DO  90  I=l,N 

XX  =  XS(  1) 

IF  (XS(I ))  40,80,70 

XX  =-XS( I) 

S(XX)  =  -1 

IF(X(I I.EQ.l) 

X( I )=-XS( I) 

S(XX)  =  -S(XXI 

GO  TO  90 

X(  I )  =  XS(  I) 

GU  TO  90 

S(  XS(  I  )  J=l 

GO  TO  50 

JX=I 

GO  TO  92 

CONTINUE 

JX  =  N*-1 

CONTINUE 

KS=1 

DO  110   I  =  1,N 

IF  (X( I n  110, 

JAD  =  X( I) 

NCX  =  NCX«-OFIJAO) 

CALL  GETYS 

CONTINUE 

CONTINUE 

PRINT  86, 

FORMAT  (• 

B=JX-1 

PRINT  82 


120,  100 


(YSdl,  1  =  1, M) 
YS«/(  2615)) 


iX(I)  ,  1=1, B) 
82  FORMAT   </•  X'/(  26l5)) 

PRINT  85  ,  (XS( I) ,  1=1,8) 
85  FORMAT   (•  XS«/(  2615)) 

RETURN 
END 


SUBROUTINE  SPUNCH 

IMPLICIT  INTEGER*2  JA-T,V-Z,$) 

INTEGER*4  NOGT,  NOVR,  NIPT,  NOPT 


COMMON 

NOVR 

f          NOGT 

,          NIPT 

, 

NOPT 

1 

U2 

U3           , 

,          U4 

COMMON 

ZbAR         , 

'     M             , 

f          N 

t 

JAD 

1 

KS           , 

f          JX           , 

r     NCX 

, 

LTS 

2 

JTS          , 

r    CLCF(20000), 

COL 

(20000) 

f 

RrtCF(20000) 

8 
10 


lb 

zo 


lb 


101 


3 

5 

6 

7 
COMMON 

1 

2  , 

D  I  MENS 
EQUIVA 

102  )  )  ,  ( 
ORDER 
DO  10 
IF  (TE 
DO  8 
/\\/ALK  = 
AVALJ= 
IF  (AV 
IP  (  AV 
ir  <AV 
ZZ =TEM 
TE:MP(  J 
Tf^MPCK 
ZZ=VAL 
VAL(J) 
VALIK) 
ZZ=TP( 
TP(J)= 
TP(K»= 
CONTIN 
CUNTIN 
CALCUL 
DO  Id 
NC=I 
IF  (  I 
CONTIN 
Tr-MP  (N 
VAL(N+ 
NA=N+1 
DO  26 
11=1-1 
INPT( 3 
INPT(3 
INPT(3 
COM  IN 
NB=NA* 
CALL  I 
DO  30 
G=  (K- 
H=K*q 
PR  INT 
FORMAT 
CONTIN 
RETURN 
END 


ION 
LENC 
CMT3 
THE 

K=l 
MP(K 

J=K 
VAL( 
VAL( 
ALK 
ALJ 
ALK 
P(  J) 
)  =  TE 
)=ZZ 
(J) 
=  VAL 
=  11 
J) 

TP{K 
ZZ 
UE 
UE 
ATE 

1=1 


RU 

PS 

CL 

X 

SU 

CL 

NC 

NU 

TE 

E 

tX 

CO 

fN 

) 

fN 

K) 

J) 

.L 

.L 

.G 


W  (20000)t 


( 

PT( 
( 

B2( 
LS( 
T5 

NULX 


^000), 
2000), 
2000) , 
2000), 
20001, 


RWPTI 
MS  ( 
CLLTC 
XS  ( 
SUP  ( 
NCT2 
NCT6 
NFIXJ 


4000) 
4000) 
2000  ) 
2000) 
2000) 


RWLTI 
CLVT( 
DF  ( 
S    ( 

TP   { 

NCT3 
NGETYS 
RWLS(  4U00) 


4000)  , 
4000) , 
2000) , 
2000  ), 
2000) 


YS  ( 
CLSP( 
XHAT( 
SUdK 

NCT4 
NNEWX 


4000) 
4000) 
2000) 
2000) 


MP(100),VAL(100),INPT(270) 

( TEMP(l)  ,Xi  1)  )  ,(VAL(1)  ,X(101) ) , (CMT1,X(2J1 ) ), (CMT2,X(2 

(203)  ),(SQ,X(204)), ( INPT( I) ,XS( 1) ) 

EFFICIENT 

.Eg.    0)    GO    TO    10 


T.  0) AVALK=-AVALK 
T.  0)AVAL J=-AVALJ 
E.    AVALJ)    GO    TO    8 


MP{K) 


(K) 


HOW    MANY    CARDS    NEEDED 
,10 


«9     .  GT  .    N )     GO    T  n    2  0 

UE 

+  1  )=-99 

1)=    0 

1=1, NA 

*I  I«-l)  =TP(  I  ) 
*I I+2)=VAL(  I) 
*I  I  +  3)=TFMP(I  ) 

UE 

3 

NTLPINtJ) 

K=1,NC 
1)*9  4-1 

101,     (TP(L),VAL(L),TEMP(L),L=G,H) , CMT 1 , C MT2 , C MT ^ , S y 
(i5X,9(Ai,I3, 14) ,3A2,I2) 

Ut 


SUBROUTINc    TABPT 

IMPLICIT  INTEGER*2  (A-T,V-Z,S) 
INTbGER*4  NUGT,  NOVR,  NIPT,  NOPT 
COMMON     NOVK         ,    NOGT 


NIPT 


NOPT 


L      f 

COMMON 
I 
2 
3 

5 
6 
7 

COMMON 
1 
2 

PRINT 


U2 

ZBAR 

KS 

JTS 

ROW  (20000) 


U3 

M 

JX 

CLCF(20000) 

RWPT(  4000) 


PS   (  4000),    NS   (  4000) 
CLPT(  2000),    CLLT(  2000) 
X    (  2000),    XS   (  2000) 
SUd2(  2000),    SUP  (  2000) 
CLLS(  2000),    NCT2 
NCT5         ,    NCT6 
NUNDLX      ,    NFIXJ 
1000, (0F(  I)  ,1,  1  =  1, N) 

•OBJECTIVE  FUNCTION'/ 


U4 

N 

NCX 

COL  (20000) 

RWLT(  4000) 


CLVT( 

OF   ( 

S    ( 

TP   ( 

NCT3 

NGETYS 

RWLS(  4000) 


4000) 
2000) 
2000) 
2000) 


1000  FORMAT ( 

1  10(13,  •  X(»,I4,  •)+•)/  (IIX,  9(13, •  X(* 

PRINT  1010 
lOlO  FORMAT  (  IHO) 
DO  100  1=  1,M 
100    CALL  ROWPT( I ) 
RETURN 
END 


14 


JAO 
LTS 

RWCF(20000) 


YS   ( 

CLSP( 
XHAT( 
SUBU 

NCT4 

NNEWX 


4000) 
4000) 
2000) 
2000) 


)♦• )  ) ) 


IM 
IN 

cn 
1 

CO 

1 

2 

i 

4 

5 

6 

7 
CO 

1 

2 
DA 
DI 
NI 
Nb 
WB 
VB 
IF 
BE 
IF 
AL 
PB 
LA 
GA 
PR 
10  FO 

1 
PR 
J= 
NB 
NE 
DO 
J= 


SU 
PL  IC 
TEGE 

MMON 

MMON 


MMON 


TA 

MENS 
PT  =  2 
G=NI 
=  1 
=  WB-»- 

INC 
TA=V 

(NC 
PHA  = 
=  ALP 
MDA  = 
MA  =  L 
INT 
RMAT 

INT 
0 

EG=W 
ND=W 

20 
J  +  1 


8R0UTINE 
IT  INTEGE 
R*4  NUGT, 
NUVR 
U2 

ZbAR 
KS 
JTS 
ROW  ( 
PS   ( 
CLPT( 
X    ( 
SU32{ 
CLLS( 
NCT5 
NUNDL 


T6PT 
R*2  (A- 

NOVR, 


20030) 
4000) 
2000) 
2000) 
2000) 
2000) 


CP 
ION 


X 

/' 

in  12) 


CP 


T,V-Z,S) 

NIPT,     NOPT     flL 
NOGT 
U3 
M 

JX 

CLCF(20000) 
KWPT(  40JJ) 
NS   (  4000) 
CLLT(  2000) 
XS   (  2000) 
SUP  (  2000) 
NCT2 
NCT6 
NFIXJ 
/ 


**NOVK 
PT*(NOGT- 


1) 


NIPT 

U4 

N 

NCX 

COL  (20000) 

R^LT {     4000) 


CLVT( 

OF   ( 

S    ( 

TP   ( 

NCT3 

NGETYS 

RWLS(  4000) 


4000) 
2000) 
2000  ) 
2000) 


NOPT 

JAO 
LTS 

RWCF(20000) 

YS       (  4000) 

CLSP(  4000) 

XHAT(  2000) 

SUBK  2000) 

NCT4 

NNEWX 


NOVR*NOGT 
T5    .NE.    CP)    VB=WB 
B+NOVR*NOGT 
T5    .NE.    CP)    VB=Wb 
8ETA+NBG*N0GT 
HA-KNOGT-l  )*NOGT 
PB«-NIPT«NOGT 
AMOA+NOGT 
10    , NOVR, NUGT 
(  5X,'THE    VARIABLE    TABLE    FOR • , I  5 , 2X , • VAR I ABLES  •  ,  I  5  , 

•       GATES    ALL    INTERCONNECTION    CASE*     //) 
1010,     (SUB,     SUB=1,N0VR) 

B 

3+N0VR*N0GT  -1 

I=NBEG,NENO,NOVR 


K=I+NOVk-l 

20  PRINT     1000,     J,( SUBf SUB=I ,K) 

1000  FORMAT     ( /9X  ,  •  W  (  '  ,  12  ,  •  ,  Y  )  =     .••,^17       ) 
IF     (NCT5    .NE.    CP)    GO    TO    IZ 

PRINT     1010, ( SUB, SUB=1,NUVR) 

J=0 

N8EG=V6 

NEND=BETA-1 

DO  21  I=NBEG,NEN0,NOVR 

J=J+1 

K= I+NOVR-1 

21  PRINT  1001,  J, (SU8,SUB=I,K) 

1001  FORMAT  ( /9X, 'VC  •, I2,« ,Y)=  :«,4I7) 

22  CONTINUE 

PRINT  1010  , ISUB,SUB=1,NIPT> 
1010  F0RMAT(*0«//15X,*Y=  :«,16I7) 

DO   40  I=l,NOGT 

DO   30  J=1,N0GT 

IF   ( I  .EQ.  J)  GO  TO  30 

JJ  =  J 

IF   (J  .GT.  I)  JJ=J-1 

flEG=8ETA+NBG*(  I - 1 ) +N I PT* ( J  J- 1) 

ENO=BEG*NIPT-l 

PRINT  1020,  I,  J,  (SUB,  SUB=BEG,END» 
1020  FORMAT  ( /6X , ' B (  • ,  12  ,  •  , ' ♦  I  2 , ' , Y  )  = 
30  CONTINUE 
40  CONTINUE 

PRINT  1010  ,  (SU8,SUB=l,N0GT  ) 

dEG=ALPHA-NOGT+l 

ENO=ALPHA-l 

DO       50        I=1,NUGT 

BEG=BEG+NGGT-1 

END=EN0+NGGT-1 

DO       4b     IZ=l,NOoT 

IF     { IZ.EU.     I)     GO    TO    43 

IF     (  I Z     .  GT  .     n     GO    TO    44 

ZZ(  IZ)=dEG-H-IZ 

GO    TO    45 

44  ZZ ( IZ)=3EG-2+IZ 
GO    TO    45 

43    ZZ( IZ)=99999999 

45  CONTINUE 

PRINT     1033,     I,     UZ{SUB> ,SUB=1,N0GT) 
1J3U    FORMAT     ( /9X ,  • A(  •  ,  I  2  ,  •  ,  Y ) =     :«,16i7) 
50    CONTINUE 

PRINT     1010, ( SUB,SUB=1,NIPT J 

BEG=P8-NIPT 

END=PB-1 

DO       60       1=1,     NOGT 

BEG=BEG+NIPT 

£NO=END+NIPT 

PRINT     1040,     I,     ( SUB,SUB=BEG,ENDI 
1040    FORMAT     ( /9X , • P (  •  ,  12 ,  •  , Y  )  =     :',16I7) 
60    CONTINUE 

IF     (NCT4    .LT.     2J     RETURN 

PRINT     1010,     (SUB,SUB=1,NUGT) 

BEG=LAMDA 

EN0=8EG*-N0GT-I 
2050    FORMAT     (/12X,»L(X)=     :*,     1217) 

PRINT    2050,     (SUB,SUB=BEG,END) 

PRINT     1010,     (SUB, S'J8  =  1, NIPT) 


DO  90  I=l,NOGT 
DO  80  J=l,NOGT 
IF   (I  .EQ.  J) 
JJ  =  J 
IH   (J  .GT.  [) 


GO  TO  80 


JJ=J-l 


BEG=GAMA*N8G*(  I-l  )*-N  I  PT*(  J  J- I) 

END=BEG+NIPT-l 

PRINT  1060,  I,  J,  ( SUB,SJB=BEG,END) 
1060  FORMAT  (/6X,  • Gl • , I  2 , • , • ,  I  2, • , Y )=  : 

80  CONTINUE 
90  CONTINUE 

RETURN 

END 


1617) 


SUBROUTINE  UNOLX  (2X, *,♦,*) 
UNDERLINE  THE  RIGHTEST  NOT  UNDERLINED 
IMPLICIT  INTEGER*2  (A-T,V-Z,i) 
INTEGER*^  NOGT,  NGVR,  NIPT,  NOPT 


XS  AND  DROP  THOSE  AT  ITS   RISHT 


NOVR 

U2 

Z6AR 

KS 

JTS 

RUW  (200001 
4000) 
2000) 
2000) 
2000) 
2000) 


NOGT 

U3 

M 

JX 

CLCF(20000) 

RWPTI  4000) 


NS   ( 

CLLT( 

XS   ( 

<^UP  ( 

NCT2 

NCT6 

NFIXJ 


COMMON 
1 
COMMON 

1  , 

2  , 

3  , 

4  ,     PS   ( 

5  ,     CLPTI 

6  ,     X    < 

7  ,     SUB2( 
COMMON     CLLS( 

1  ,     NCT5 

2  ,     NUNDLX 
42  FORMAT  (•  UNDLX   ENTER*  ) 

1  IFIJX.6U.I)  GO  TO  12 

SUB=-XS(JX-I ) 

IF  (SUB  .LT.  0)  SUB=-SUB 

BEG=CLPT(SUB+l) 

EN0=CLLT(SUB+1 )+BEG-l 

DO  3  AA=BEG,END 

IF  (PSCROWCAA))  .GE.  0)  GO  TO  3 

IF  (CLSPIROWCAA) )  .N£.  0)  GO  TO 

CLSP{NCT3)=R0W( AA) 

NCT3=R0W(AA) 

CLSP(NCT3)=R0W(AA) 
3  CONTINUE 
5  L  =  JX-l 

IF(XS{L).GT.O)  GO  TO  15 

S(-XSI  L)  )  =  0 

XS(L)=0 

IF( X(L).LE.U) 

J^D=X(L) 

KS  =  -l 

GO  TO  20 

10  JAD  =  -X(L  ) 
KS  =  1 

20  X(L)=0 

CALL  RVPSNS 

JX=JX-1 

IF  (JX.GT.l) 

11  CONTINUE 

12  PRINT  13,  (I, 


4000) 
2000) 
2000) 
2000) 


GO  TO  10 


GO  TO  5 


NIPT 

U4 

N 

NCX 

COL  (20000) 

RWLT(  4000) 


CLVT( 

OF   ( 

S    ( 

TP   ( 

NCT3 

NGETYS 

RWLS(  4000) 


4000) 
2000) 
2000) 
2000) 


NOPT 

JAO 
LTS 
RWCF(20000) 


YS  ( 
CLSP( 
XHAT( 
SUBU 

NCT<, 
NNEWX 


4000) 
4000) 
2000) 
2000) 


XHAT(I),  1=1, N» 


0)  GO  TO  16 


13  FORMAT  (  11(2X,«X( 
PRINT  14,  ZBAR 

l^    FORMAT  (2X,I4) 
43    FORMATC  ONDLX  EXIT  • 
RETURN  Z 
XS(L)=-XS(L) 
IF  (X(L)  .GT 
JAO=-X(L) 
KS=1 

GO  TO  17 
JAO=X( L) 
KS  =  -1 

X(L)=-X(L) 
CALL  RVPSNS 
JX  =  JX-l 

CALL  CGPSNS  (ZX,&1,£200) 
LTS  =  0 
RETURN  I 
200  RETURN  3 
END 


13,'  )  =  •  ,iin 


18) 


15 

130 


lb 


17 


SUBROUTINE  GTINEQ 

GENERATE  INEQUALITIES  FOR  AND-OR  NE 

THE  CARDS  MARKED  WITH  BLUE  ON  THE  T 

ADDITIONAL  INEQUALITIES  WHOSE  EFFEC 

THESE  CHANGES  WERE  PUT  IN  THE  PROGR 

MAKE  ALL  FOUR  VERSIONS  OF  GTINEQ  UN 

THESIS-  APPENDIX  B). 

IMPLICIT  INTEGER*2  (A-T,V-Z,$» 

INTEGER*^  NOGT,  NOVR,  NIPT,  NOPT 


T WORKS. 

OP  IMPLEMENT  CHANG 
TIVENESS  HAVE  NOT  B 
AM  ON  MAY  I,  197^  I 
IFORM  (SEE  HOHULIN 


COMMON     NOVK         ,    NOGT 
1      ,     U2  f    U3 

COMMON      ZBAR         ,    M 

1  •      KS  »    JX 

2  f     JTS         ♦    CLCF(20000 

3  t     ROW  (20000),    RWPT(  ^000 

4  f     PS   (  4000),    NS   (  <»000 

5  f     CLPT(  2000),    CLLT(  2000 

6  «     X    (  2000),    XS   (  2000 

7  f     SUB2(  2000),    SUP  (  2000 
COMMON     CLLS(  2000),    NCT2 

1  ,     NCT5         ,    NCT6 

2  ,      NUNDLX       ,    NFIXJ 
DIMENSION  T£MP( 100) ,VAL( I00),INPT(2 
EQUIVALENCE  ( TEMP ( 1 ) , X ( I ) ) , ( VAL ( I ) , 

102  )  )  ,(CMT3,X(203)  )  ,(  Sg,X(  20'»)  )  ,(INP 
DATA  Bl,B2,Gi,G2,OK  /'Bl', 

DATA   W,V,D,C,B  /'W*,* 

DATA   NR,ND,AO,GT,ST,IN,GT,OU,FI ,F0 


70) 

X(10 

T(l) 

•B2' 


NIPT 

U4 

N 

NCX 

COL  (200001 

RWLT(  4000) 

CLVT(  4000) 

OF   (  2000) 

S    (  2000) 

TP   (  2000) 

NCT3 

NGETYS 

RWLS(  4000) 


ES  IN  THE 
EEN  TESTED. 
N  ORDER  TO 
MASTER 


NOPT 

JAD 
LTS 

RWCF(20000) 

YS   (  4000) 

CLSPC  4000) 

XHAT(  2000) 

SUBK  2000) 

NCT4 
NNEWX 


I 


,  •  AD« 
/•CL« 


GT 
SI 


IN' 
/ 


,BLK 
,  'GT 


1)1  ,(CMT1,X(201)) ,(CMT2,X(2 
,XS(  I  )) 

,*Gl'  ,'G2' , 'OR'  / 

D',  'C,  '  •  / 

/ 
•      / 


DATA   CL  ,SI 
DATA  E   /'E'/ 
PRINT  1011 
1011  FORMAT  CI') 
C      GENERATE  OBJECTIVE  FUNCTION 
00  70  1=1, NOGT 
DO  71  J=1,N0VR 
OF ($V(  I,  J)  )=1 

71  OF(iW( I, J) )=1 
DO  72  JJ=1,N0GT 

IF  (JJ.EQ.I )  GO  TO  72 
QF($A(  I  ,JJ)|=1 

72  CONTINUE 
70  CONTINUE 

N0P1=NUPT+1 
U=NnGT-i-N0VR  +  2 
0NE  =  1 
CCCCCCCBASIC  INEQULITIFSCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
UO  100  1=1, NOGT 
CALL  CLEAR 
CMT1=AD 
CMT2=BLK 
CMT3=GT 
SQ=I 

TEMPID^O 
VAL(l)=-l 
TP(n=C 
TEMP(2)  =  $L(  I  ) 
VAL(2) =U 
TP(2)=C 

TEMP(3)=*P( I, ONE) 
VAL(3)=U 


TP( 3)=0 

TEMP(<^^  =  $W(I  ,ONE» 

VAL(4) =1 

TP(4)=V 

TeMP(5)=SV(  I,f3NE) 

VAL(5) =1 

TP(5)=W 

N=5 

DO  101  J=1,N0GT 

IF  (J.EO.I)  GO  TO  101 

N=N+1 

TEMP(IM)=$G(  J,I,ONE» 

VAL(N)=1 

TP(N)=D 

101  CONTINUE 
CALL  SPUNCH 
00  102  AA=l,N 

102  VAL( AA)=-VAL( AAI 
VAL(2)=-VAL(2) 
VAL(1)=U 

CALL  SPUNCH 

CALL  CLEAR 

CMT1=0R 

TEMP( 1 )=0 

VAH1)=2*U-1 

TP(1)=C 

TEMP(2)=$L( I ) 

VAL(2) =-U 

TP(2)=C 

TtMP(3)=$P(I,0NE) 

VAL(3)=-U 

TP( 3)=D 

TtMP{4)=$W( I tONEI 

VAL(4) =1 

TP(4)=W 

T£MP(5)=$V(I,0Nt) 

VAL(5)=1 

TP(5)=V 

N=5 

on  151  J=1,NUGT 

IF  (J.tO.I)  GO  TO  151 

N=N+l 

T£MP(N)=$0( J,I,CNE) 

VAL (N)=l 

TP(N)=0 

151  CONTINUE 
CALL  SPUNCH 
00  152  AA=1,N 

152  VAL{AA)=-VAL(AAJ 
VAL( 1)=U 
VAL(2)=-VAL{2) 
CALL  SPUNCH 

100  CUNTrNUE 
CMT1=AD 
CMT2=0R 
CVT3=0U 

00  200  I=l,i\OGT" 
DO  201  J=1,N0GT 
IF  (J. Eg. I )  GO  TO  201 
CALL  CLEAR 
SQ=I*10+J 


TEMP(1)=$A(I,J) 
VALC l)=-2 
TP(1)=C 

TEMP(2)=$P(I ,ONE) 
VAL(2l=l 
TP(2)=D 

TeMP(3)=$G(I ,J,ONE» 
VAL(3)=2 
TP(3)=D 

TEMP(4)=$B(I ,J,ONE» 
VAL(4)=1 
TP(4)=D 
N=4 

CALL  SPUNCH 
CALL  CLEAR 
TEMP(1)=0 
VAL( l)=l 
TP(1)=C 

TEMP(2)=$A(I , J) 
VAL(2I=3 
TP(2)=C 

TEMP(3)=$G(I ,J,ONE» 
VAL(3)=-^ 
TP(3)=D 

TEMP(*|=$P( I ,ONE> 
VAL(4)=-1 
TP(<»)=D 

TEMP(5)=i6( I ,J,ONE) 
VAL (5)=-3 
TPi5)=0 
N=5 

CALL  SPUNCH 
201  CONTINUE 
200  CONTINUE 
CMT3=GT 

DO  250  I=l,NOGT 
CALL  CLEAR 
SO=I 

TEMPI! )=$P( I ,ONE) 
VAL(1)=U 
TPa)  =  0 
N=l 

00  251  J=1,N0GT 
IF  ( I.EO.J)  GO  TO  251 
N=N+1 

TEMP(N)  =  $ti(  I  ,  J,ONE) 
VAL(N)=-1 
TP(NI=D 
251  CONTINUE 

CALL  SPUNCH 
250  CONTINUE 
CCCCCCtACH  GATE  HAS  AT  LEAST 
CMT3=IN 

DO  300  I=1,N0GT 
SQ=I 

CALL  CLEAR 
TEMPI! )=0 
VAL(l)=-2 
N=l 

DO  301  L  =  UNOVR 
N=N  +  l 


TWO  INPUTCCCCCCCCCCCCCCCCCCCCCCCCCCC 


TEMP(N)=$W(I ,L) 
VAL(N)=l 

N  =  N+l 
TEMP(N)=$V(I ,LI 

301  VAL(NI=l 

DL)  30^  J  =  1,N0GT 

IF  (J.EQ.I )  GO  TU  302 

N=N+1 

TFMP(N)=SA(J, I ) 

VAL(N)=1 

302  CONTINUE 
CALL  SPUNCH 

300  CONTINUE 
CCCC   ALL  GATES  EXCEPT  THE  OUTPUT  GATES  CONNECT  TO  AT  LEAST  ONE  GATE  CCC 

C'»^T3  =  0U 

DO  400  I=N0P1,NCGT 

CALL    CLEAR 

SQ=I 

TEMPC l»=0 

VAL(i)=-l 

N=l 

DO    401    J=1,N0GT 

IF     (J. £0,1)    GO    TO    401 

N=N+l 

TfcMP(N)=$A( I , J) 

VAL(N)=1 
401  CONTINUE 

CALL  SPUNCH 
400  CONTINUE 
CCCCCCTkAINGULAR  INEOUALITIESCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C       IJ  =  2  FOR  SINGLE-OUTPUT  NETWORKS;    IJ  =  1  FOR  MULTIPLE-OUTPUT 
C      NETWORKS.   PREVIOUSLY  ONLY  IJ=1  WAS  USED. 

IJ  =  2 

IF  (  NOPT  .GT.  1  )  IJ  =  I 

C^^Ti  =  AD 

CMT2=0R 

CMT3=ST 

su=o 

DO  600  I=IJ,N0GT 

00  600  J=IJ,NOGT 

IF  (J. EW. I )  GO  TO  600 

DO  601  K=1,N0GT 

IF  (K.EQ.I .UR.K.tO.J)  GO  TO  601 

CALL  CLEAR 

S0=SQ+1 

IF  (SO.GE.IOO)  SQ=l 

TEMPd  )=0 

VAL{1)=2 

TFMP12)=$A( I , J) 

VALl?)=-l 

TEMP(3»=$A(I ,K) 

VAL(3)=-1 

TEMP(4)=$A( J,K) 

VAL<4)=-1 

N=4 

call  spunch 
601  continue 
600  continue 

cmti=tw 

CMT2=SM 
CMT3=TY 


DO  900  I  =  N0P1,N0GT 

00  900  J=1,N0GT 

IF  (I.EQ.J)  GO  TO  900 

SQ=I*10*J 

CALL  CLEAR 

TEMP(l)=0 

VAL( 1)=0 

TEMP(2)  =  $LU  ) 

VAL{2)=1 

TEMP(3)=$L( J) 

VAL(3)=1 

TEMP(4)=$A(l, J) 

VAL(^I=-1 

N=4 

DO  901  K=l,NOGT 

IF  (K.EO.I.OR.K.EO. J)  GO  TO  901 

N=N-H 

TEMP(NI=$A( I ,K) 

VAL(N»=I 
901  CONTINUE 

CALL  SPUNCH 

VAL(I)=2 

VAL(2)=-VAL(2> 

VAL(3)=-VAL(3) 

CALL    SPUNCH 
900    CONTINUE 

CMTI=1 

CMT2=AD 

CMT3=1 

SQ^OR 

DO    1000    I=I,NOGT 

N=N+I 

TEMP(N)=$L( I ) 

1000  VAL(N)=1 
CALL    SPUNCH 
VALIU=NGGT-1 
DO     1001    AA=2,N 

1001  VAL(AA)=-VAL(AA| 
CALL    CLEAR 
TEMPI  1 1=0 
VALU)=-l 

N=l 

CALL    SPUNCH 

CALL    CLEAR 

INPT(1)=E 

INPT(6)=-99 

N8  =  6 

CALL     INTLP(NB) 

RETURN 

END 

SUBROUTINE    GTINEQ 
C  GENERATE    INEQUALITIES    FOR    NOR-ANO    NETWORKS. 

C  THE    CARDS    MARKED    WITH    BLUE    ON    THE    TOP    IMPLEMENT       CHANGES    IN    THE 

C  ADDITIONAL    INEQUALITIES    WHOSE    EFFECTIVENESS    HAVE    NOT    BEEN    TESTED, 

C  THESE    CHANGES    WERE    PUT    IN    THE    PROGRAM    ON    MAY    1»     1974    IN    ORDER    TO 

C  MAKE    ALL    FOUR    VERSIONS    0*=    GTINEQ    UNIFORM    (SEE    HOHULIN    MASTER 

C  THESIS-    APPENDIX    8). 

IMPLICIT     INTEGER*2    (A-T»V-Z,$) 

INTEGER*^    NOGT,     NUVR,     NIPT,     NOPT 


COMMON     NOVR 
I      ,     U2 
COMMON     ZBAR 

1  ,     KS 

2  ,     JTS 

3  ,  ROW  ( 

4  ,     PS   ( 

5  ,     CLPT( 

6  ,     X    ( 

7  »  SUB2( 
COMMON     CLLS( 

1  ,     NCT5 

2  ,  NUNOL 
DIMFNSION  TEMP( 
EQUIVALENCE  (TE 

102) ) , (CMT3,X(20 


NIPT         f    NOPT 

U4 

N 

NCX 

COL  (200001 

RWLT(  4000) 

CLVTI  4000) 

OF   (  2000) 
S    (  2000), 

TP   (  2000) 

NCT3 

NGETYS 

RWLSI  4000) 
100) ,VAL( 100) ,INPT{270) 

MP(1),X( 1) ) , (VAL(l) ,X(101) ).(CMTi,X(201 ) ) , (CMT2,X(2 
3)  )  ,(SQ,X(204)), (INPT<1),XS(1)) 


20000) 
4000) 
2000) 

2000) , 
2000) 
2000) 


NOGT 
U3 
M 

JX 

CLCF(20000) 
RWPT(  4000) 
NS  (  4000) 
CLLT(  2000) 
XS  (  2000), 
SUP  (  2000) 
NCT2 
NCT6 
NFIXJ 


JAD 
LTS 

Ri<CF(20000) 
YS  (  4000) 
CLSP(  4000) 
XHAT(  2000) 
SUBK    2000) 

NCT4 

NNEk^X 


OATA 
OATA 
OATA 
DATA 
DATA 


E/'E*/ 
CP 

bl,B2,Gl,G 
W,V,D,C,B 
NR,NO, AD, 


I 


OATA 
PRINT 
1011  FORMAT 


CL  ,SI 
1011 
I  •!•  ) 


/•CP'  / 

2  /'Bl* ,«B2',-oi- , 

/• W« ,«V», 'D* ,'C»,«  • 
GT,ST, IN,GT,0U,FI,F0,BLK 

,  • A0« ,»GT»  ,«ST«  ,»IN« ,  'GT* , •0U», 'FI* ,»F0» , • 
/•CL»,'SI«         / 


FUNCTION 


73 

71 


72 

70 


GENERATE    OBJECTIVE 
DO    70     1=1, NOGT 
on    71    J=1,N0VR 
IF     (NCT5    .NE.    CP)    GO    TO     M 
OF     (     $V(I,J))=1 
CONTINUE 
GF($W(  I,J) )=l 
00    72    JJ=1,N0GT 
IF     IJJ.EQ.I )     GO    TO    72 
OF($A(  I, JJ))  =  1 
CONTINUE 
CONTINUE 
NnPl=NOPT+l 
U=.NlUGT-»-N0VR+2 
0N£=1 
Tv<0  =  2 
THRFE=3 
CCCCCCCBASIC     INEOULITIESCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO    100     1=1, NOGT 
DO     10    AA=1,100 
TEMP(AA)=0 
VAL (AA)=0 
10    TP{AA)=BLK 
CMT1=NR 
CMT2=BLK 
CMT3=GT 
S0=I 

TEMP(1)=0 
VAL(  l)=U~l 
TP(1)=C 
TEMP(2)=tL(I ) 
VAL(2)=-U 
TP(2)=C 

TEMP(3)=$P(I ,ONE) 
VAL(3)=U 
TP(3)=0 


I 


TEMP{4)=$WII ,ONE) 

VAL(4)=1 

TP(4)=W 

IF  (NCT5  .NE.  CP  )  GO  TO  105 
TEMP15»=$V{I ,ONE» 
VAL( 5)=1 
TP(5)=V 

N=5 
105  CONTINUE 

DO  lOl  K=1,N0GT 

IF  (K.EQ.I)  GO  TO  lOI 

N=N+I 

TEMP(N)=$B(K, I,ONEI 

VAL(N)=1 

TP(N)=D 

101  CONTINUE 
CALL  SPUNCH 

VAL(  n=2*u 

DO  102  AA=3,N 

102  VAL(AA)=-VAL{AA) 
CALL     SPUNCH 
CALL    CLEAR 
CMT1=AD 

TEMP(  1)  =  0 

VAL (l)=-l 

TP(1)=C 

TEMP(2)  =  SL(  I  ) 

VAL(2)=U 

TP(2)=C 

TEMPO  )  =  $P(  I, ONE) 

VAL(3) =U 

TP(3J=D 

TEMP(4)  =  $,^(  I  ,IJNEI 

VAL(^)=1 

TPI^|=V 

N=4 

IF     (NCT5    .NE.    CP    )    GO    TO    150 

TEMP( 5)=tV( I ,ONE) 

VAL(5)=l 

TP(5)=W 

N=5 

150  CONTINUE 

DO  151  K=1,N0GT 

IF  (K  .EQ.  I )  GO  TO  151 

N=N-»-l 

TEMP(N)=$G(K,I ,ONE) 

VAL(N)=1 

TP(N)=D 

151  CONTINUE 
CALL  SPUNCH 
VAL(1)=U 

DO  152  AA=3,N 

152  VAL(AA)=-VAL CAA) 
CALL  SPUNCH 

100  CONTINUE 
CMT1=NR 
CMT2=AD 

DO  200  I=1,N0GT 
DO  201  J=1,N0GT 
IF  (J.EO.I )  GO  TO  201 


CALL  CLEAR 

S0=«I*10*J 

TEMP( 1  )  =  $A( I 

,JI 

VAL(l) =-2 

TP(  l)  =  C 

TEMP{2)=$P(I 

,ONEI 

VAL(2)=1 

TP(2)=D 

TEMP(3)=JG(I 

,J,ONE) 

VAL(3)=2 

TP(3)=D 

TEMP(^)=$fl{ I 

,J,ONE) 

VAL(^)=1 

TP(A)=D 

N=4 

CALL  SPUNCH 

CALL  CLEAR 

TEMPd  )=0 

VAL(l»=l 

TP(1)=C 

TeMP(2 )=$A{I 

,J) 

VAL(2)=3 

TP(2)=C 

TeMp(3)=tG( I 

,J,ONE) 

VAL(3)=-4 

TP(3)=D 

TEMP{4)=$P(I 

,0N£) 

VAL(4)=-1 

TP{4)=D 

TEMP(5)=$B(I 

,J,ONE) 

VAL(5)=-3 

TP(5)=D 

N=5 

CALL  SPUNCH 

201  CONTINUE 

200  CONTINUE 

CMT3=GT 

OG  250  I=I,NOGT 

CALL  CLEAR 

sg=i 

TEMPd  )=$P(  I 

,aNE> 

VAL{1)=U 

TP(ll=D 

N=l 

00  251  J=l,NOGT 

IF  (I.EQ.J) 

GO  TO  251 

N=N+I 

TEMP(N)=$8( I 

, J.ONE) 

VALlN)=-l 

. 

TP(N)=D 

251  CONTINUE 

CALL  SPUNCH 

250  CONTINUE 

CCCCCCEACH  GATE  HAS  AT  LEAST 

CMT1=NR 

CM-^2=A0 

CMT3=IN 

00  300  I=1,N0GT 

sg=i 

CALL  CLEAR 

TtMP( 1 )=0 

ONE  INPUTCCCCCCCCCCCCCLCCCCCCCCCCCCC 


VAL(n=-2 

TEMP(2»  =  $LCn 

VALC2l=l 

N=2 

DO    301    J  =  UNOVR 

N=N-H 

TEMP(N)=$W(I , J) 

VAL(N)=1 

IF     (NCT5    .NE.    CP     I    GO    TO    301 

N=N+1 

TEMP(N)=$V(I ,J) 

VAL(N»=1 

301  CONTINUE 

00  302  K=l,NUGT 

IF  (K  .EO.  I)  GO  TO  302 

N=N-H 

TEMP(N)=SA(K, I ) 

302  VAL(N)=1 
CALL  SPUNCH 

300  CONTINUE 

CCCC   ALL  GATES  EXCEPT  THE  OUTPUT  GATES  CONNECT  TO  AT  LEAST  ONE  GATE  CCC 

CMT3=0U 

DO  350  I=N0P1,N0GT 

CALL  CLEAR 

TEMPin=0 

VAL(  l)=-l 

N=l 

SQ=I 

DO  351  K=1,N0GT 

IF  (K  .EO.  I )  GO  TO  351 

N=N+1 

TEMP(N)=$A( I,K| 
351  VAL{N)=l 

CALL  SPUNCH 
350  CONTINUE 
CCCCCCTRAINGULAR  I NEQUAL I T  lESCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C      IJ=2  FOR  SiNGLE-OUTPur  NETWORKS;   IJ=1  FOR  MULTIPLE-OUTPUT 
C       NETWORKS.   PREVIOUSLY  ONLY  IJ=2  WAS  USED. 

IJ  =  2 

IF  (  NUPT  .GT.  1  )  IJ  =  1 

CMT3=ST 

so=o 

DO  400  I=IJ,NOGT 

00  401  J=IJ,NOGT 

IF  (J  .FQ.  I )  GO  TO  401 

DO  402  K=l,NOGT 

IF  (K  .GQ.  J  .OR.  K  .EO.  I )  GO  TO  402 

CALL  CLEAR 

SQ=S0+1 

IF  (SQ  .GE.  IJO)  S0=1 

TEMP( 1)=0 

VAL(  n=2 

TEMP(2  I^SAU  ,  Jl 

VAL(2)=-1 

TEMP(3)=*A( I,Kl 

VAL(3) =-1 

TEMP{4)=$A(J,K| 

VAL(4)=-l 

N=4 

DO  403  P=l,NOGT 

IF  (P  .Ey.  I  .OR.  P  .EO.  J  .OR.  P  .to.  K)  GO  TO  403 


TEMP(N)  =  $A(  J.PJ 

VAL(N)=1 
^03  CONTINUF 

CALL  SPUNCH 
402  CONTINUE 
401  CONTINUE 
400  CONTINUE 

DO    500    J=IJ,NOGT 

DO    501    K=1,N0GT 

IF     (K     .to.    J )     GU    TO    501 

DO    502    L  =  UNnVR 

SQ=S0+1 

IE     (SO    .GE.     100)     SU=l 

CALL    CLEAR 

TEMP(  1  )  =  0 

VALI l)=2 

Tt MP(2 )=$W(K,L» 

VAL(2l =-1 

TEMP(3)=$W( J,L) 

VAL(3) =-1 

T£MP(4)=$AIJ,K) 

VAL(4)=-l 

N=4 

DO    503    P=l,NOGT 

IF     (P    .Ey.    K     .OK,    f 

N=N+1 

TEMP(N )=$A(J,P) 

VAL(N)=1 
503    CONTINUE 

CALL     SPUNCH 

IF  (  NCT5  .NE.  CP  )  GO  TO  502 

TtMP(2)  =  $V(J»L) 

TE^^P(3)  =  tV(K,L) 

CALL  SPUNCH 
502  CONTINUE 
501  CONTINUE 
500  CONTINUF 

00  600  J=IJ,NOGT 

DO  601  K=1,N0GT 

IF  (K  .EQ.  J)  GO  TO  601 

00  602  L=UNOVR 

S0=.':Q+1 

IF  (SQ  .GE.  100)  SQ=1 

CALL  CLEAR 

TEMPd  )=0 

VAL{1)=2 

TEMP(2)=$W( J,L) 

VAL{2) =-1 

TFMP(3»=$W(K,L) 

VAL(3»=-1 

TEMP{4)=$A( J,K) 

VAL{4) =-1 

TEMP(5  )  =  $L( J) 

VAL(5)=1 

N=5 

CALL  SPUNCH 

Su=SU+l 

IF  (SQ  .Gt.  100)  SQ=1 

TEMP(N )=$L(K) 

CALL  SPUNCH 


EQ.  J  )  GO  TO  503 


IF     (     NCT5    .NE.     CP    )     GO    TO    602 

TEMP(2)     =    $V(J,L) 

TEMPO!    =    $V(K,L) 

CALL    SPUNCH 

TEMP(N)     =    $LCJ) 

CALL  SPUNCH 
602  CONTINUE 

DO  650  I=IJ,NOGT 

IF  (  I  .EQ.  J  .OR.  I  .EO.  K)  GO  TO  650 

SU=S0+1 

IF  (SQ  .GE.  100)  SO=l 

TEMP(1)=0 

VAL(U=2 

TEMP(2»=$A| I,K» 

VAL(2J=-1 

TEMPO  J  =  SA(  I, J) 

VAL (i)=-l 

TEMP(4)=$A( J,K) 

VAL(4)=-I 

TEMP(5)=$L( J) 

VALI 5) =1 

N=5 

CALL  SPUNCH 

SQ=SO+l 

IF  (SU  .GE.  100)  SQ=l 

TEMP(N)=$L(K) 
CALL  SPUNCH 
650  CONTINUE 
601  CONTINUE 
600  CONTINUE 
CCCC   GATES  WHICH  ARE  CONNECTED  TO  ALL  OF  THE  OUTPUT  GATES  ARE  NOT 
CCCC   CONNECTED  TO  ANY  OTHER  GATES  CCC 
CMT3=CL 

DO  700  I=N0P1,NCGT 
SU=I 

CALL  CLEAR 
TEMPil )=0 
VALI 1) =NOPT*U 
00  695  J=2,N0PI 
JMl=J-l 

TEMP( J)=$A{ I , JMl) 
VAL{J)=-U 
N=J 
695    CONTINUE 

DO  701  K=NOPl,NnGT 
IF  (K  .eg.  I )  GO  TO  701 
N=N*-1 

TEMP(N)=$A( 1 ,K) 
VAL IN) =-1 
701  CONTINUE 

CALL  SPUNCH 
700  CONTINUE 

CALL  CLEAR 

INPT( 1 )=E 

INPT(6)=-99 

N8=6 

CALL  INTLP(NB) 

RETURiNj 

END 
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CCCCCCCBASIC 
DO  100 
DO  10 
TFMP{  A 
VAL{ AA 
10  TP(AA) 
CMT1=N 
CMT2=B 
CMT3=G 

sa=i 
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1,B2,G1,G2,0R  /• 
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NIPT 
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0  00) 


73 
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72 

70 


CP 
CL  ,SI 

{  'l') 
1011 

TE  OBJECTIVE 

1=1, NOGT 

J=1,NGVk 

T5  .NE.  CP)  GO 

V(  I, J)  )  =  l 

UE 

I  ,  J)  )  =  1 

JJ=1,NQGT 

.EQ.I)  GO 

I  ,  JJ) )=l 

UE 

UE 

OPT+1 

♦NOVR+2 


/•CP* 
/•CL«  ,«<;i' 


FUNCT luN 


PT(270) 
(1)  ,X{10 
(  INPT(l) 
BIS  •B2« 

W« ,'V«,«I 
I ,FO,BLK 
•  IiN*  ,  'GT 
/ 
/ 


1  )  )  ,(CMT1,X( 201 ) ), ( CMT2,X(2 
,XS(  D) 

, 'Gl* ,«G2« , 'OR'  / 

l)«,»C«,«  •  / 

/ 
•      / 


TO  73 


TO  7  2 


2 
3 

I 
I 
AA 
A) 
)  = 
=  H 
R 

LK 
T 


NEQULITIESCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

=1 ,NOGT 

=1,100 

=  0 

0 

LK 


TEMPCl J=0 
VALC1I=U-1 
TPU)  =  C 
TEMP(2)=$L(I) 

VAL(2)=-U 
TP(2»=C 

TEMP(3I=$P(I,0NE) 
VAL(3»=U 
TP(3)=D 

TEMP(<t)  =  $W(  I, ONE) 
VAL(^)=1 
TP(4)=W 
N=4 

IF  (NCT5  .NE.  CP  )  GO  TO  15 
TEMP(5)=$V( I,ONEI 
VAL(5) =1 
TP(5I=V 
N=5 
15  CONTINUE 

00  20  K=l,NaGT 
IF  (K  .EQ.  I )  GO  TO  20 
N=N+1 

TEMP(NJ=$B(K,I ,ONE) 
VAL(N)=1 
TP(N)=D 
20  CONTINUE 

CALL    SPUNCH 
VAL(l) =2*U 
DO     30    AA=3,N 
30    VAL ( AA)=-VALIAA) 
CALL    SPUNCH 
CALL    CLEAR 
CMTl=ND 
TEMPd  )=0 
VAL( 1)=U-I 
TP(n  =  C 
TEMP(2»=$L( I) 
VAL(2)=U 
TP(2)=C 

TEMP(3)=$P( I ♦ONE) 
VAL(3)=-U 
TP(3)=D 

TEMP(4)=$WI I ,ONE) 
VAL(4)=1 
TPI4)=V 
N=4 

IF  (NCT5  .NE.  CP  )  GO  TO  35 
TEMP(5)=»VII ,ONE) 
VAL(5)=I 
TP(5)=W 
N=5 
35  CONTINUE 

DO  40  K=l,  NOGT 
IF  (  I  .EQ.K)  GO  TO  '^O 
N=N  +  1 

TEMP(N)=$G(K, I, ONE) 
VAL(N)=1 
TP{N)=D 
40  CONTINUE 

CALL  SPUNCH 
VAL( I)=0 


DO  50  AA=3,N 

50  VAL(AA)=-VAL(AAJ 

CALL  SPUNCH 

100  CONTINUE 

CMT1=NR 

CMT2=ND 

00  200  I=1,N0GT 

00  201  J=1,NUGT 

IF  (J.EU.I) 

GO  TO  201 

CALL  CLEAR 

SQ=I*10+J 

TEMP(1)=$A(I 

,J) 

VAL( l»=-2 

TP(1)=C 

TbMP(2)=$P(i 

,ONtl 

VAL(2)=1 

TP(2)=D 

TEMP(3)=$G( I 

,J,ONE) 

VAL( 3) =2 

TP(3)=D 

TFMP(^»=$B(I 

, J,ONE) 

VAL(4)=I 

Tp(4)=D 

N=4 

CALL  SPUNCH 

CALL  CLEAR 

TEMPd  )  =  0 

VALID  =1 

TP  {l)  =  C 

TEMP(2 )=$A( I 

,J) 

VAL(2)=3 

TP{2)=C 

TEMP13)  =  $G{  I 

tJ.GNEI 

VAL( 3J=-^ 

TP(3)=D 

TEMP(''t)=$P(l 

,ONE) 

VAL(4)=-1 

TP(^)=D 

Tf-MH(5)=$B(  I 

,J,ONE) 

VAL(5) =-3 

TP(5)=D 

N=5 

- 

CALL  SPUNCH 

201  CONTINUE 

200  CONTINUE 

CMT3=GT 

DO  250  I=l,NOGT 

CALL  CLEAR 

SQ=I 

TtMPd  )=$P(  I 

,ONE) 

VAL(n=U 

TP(U  =  D 

N=l 

DO  251  J=1,N0GT 

IF  lI.EQ.J) 

GO  TO  251 

N=N+1 

TEN'P(N)=tl3(I 

,J,aNE) 

VAL(NI=-1 

TP(N)=D 

251  CONTINUE 

CALL  SPUNCH 

250  CONTINUE 
CCCCCCEACH  GATE  HAS  AT  LEAST  ONE  I NPUTCCCCCCCCCCCCCCCCCCCCCCCCCCC 

CMTl=NR 

CMT2=ND 

CMT3=IN 

DO  350  J=1,N0GT 

SO  =  J 

CALL  CLEAR 

TEMPd  )  =  0 

VAL(l)=-l 

N=l 

00  310  K=l,NOVR 

N=N*l 

TEMP(NI=$W{ J,KI 

VAL(N)=1 

IF  (NCT5  .NE.  CP  )  GO  TO  310 

N=N+1 

TEMP(NI=SV(J,K) 

VAL(N) =1 
310  CONTINUE 

DO  320  I=1,N0GT 

IF  (I  .EQ.  J)  GU  TO  320 

N=N  +  1 

T6MP(N)=$A( I, J) 
320  VAL(N)=1 

CALL  SPUNCH 
350  CONTINUE 

CMT3=0U 
CCCC   ALL  GATES  EXCEPT  THE  OUTPUT  GATES  CONNECT  TO  AT  LEAST  ONE  GATE  CCC 

DO  400  I=NOPl,NOGr 

CALL  CLEAR 

TEMP(  1  )  =  0 

VAL(l) =-1 

N=l 

SO=I 

00  360  J=1,N0GT 

IF  ( I  .EO.  J)  GO  TO  360 

N  =  N  +  l 

TEMP(N)  =  tA(  I  ,  J) 
360  VAL(NJ=l 

CALL  SPUNCH 
400  CONTINUE 
CCCCCCT^AINGULAk  I Nt QUALI T I ESCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C      IJ=2  FOR  SINGLE-OUTPUT  NETWORKS;   IJ=1  FOR  MULTIPLE-OUTPUT 
C      NETWORKS.   PREVIOUSLY  ONLY  IJ=2  WAS  USED. 

IJ  =  2 

IF  (  NOPT  .GT.  I  )  IJ  =  1 

CMT3=ST 

S0  =  0 

DO  450  I=lJ,NOGT 

DO  440  J=IJ,NOGT 

IF  (I  .EQ.  J)  GO  TO  440 

DO  430  K=1,N0GT 

IF  (K  .EQ.I   .OR.  K  .EQ.  J)  GO  TO  430 

CALL  CLEAR 

SU=SU+1 

IF  (SO  .GE.  100)  S0=1 

VAL( l)=2 

TEMP(l)=0 

DO  405  AA=2,4 
405  VAL(AA»=-1 


TEMP12)=$A( I , J) 
TEMP(3)=$A( J,K| 

TEMP(^)=$A( I ,K) 

DO  410  P=1,N0GT 

IF  I    P  .EQ.  I  .UR.  P  .EQ.  J  .OK.P  .EQ.  K)  GU  TO  410 

N=N+1 

TEMP(N )=SA( J,P) 

VALlN)=l 
410  CONTINUE 

CALL  SPUNCH 
430  CONTINUE 
440  CONTINUE 
450  CONTINUE 

DU    500    J=IJ,NUGT 

DO    490    K=1,N0GT 

IF       (J     .EQ.    K)    GO    TO    490 

no    480    L=1,N0VR 

SQ=SQ*1 

IF     (SQ    .GE.     100»    Sg=l 

CALL    CLEAR 

VALm=2 

TEMPd  1  =  0 

DO    455     AA=2,4 
455    VAL{AA)=-1 

TEMP(2)=$W( J,L) 

TEMP(3i=$W(K,L) 

TEMP(4)=$A( J,K) 

N=4 

OU    470    P=1,N0GT 

IF     (P.EQ.J     .OR.    P.EQ.K)    GU    TO    470 

N=N+1 

VAL(N»=1 

TEMP(N)=$A( J,P) 
470  CONTINUE 

CALL  SPUNCH 

IF     (    NCT5    .NE.    CP     )    GO    to    480 

TEMP(2)     =    $V(J,L) 

TEMPO)     =    $V(K,L) 

CALL  SPUNCH 
480  CONTINUE 
490  CONTINUE 
500  CONTINUE 

00    600    J=IJ,NOGT 

DO    590    K=1,N0GT 

IF     {J     .FO.    K)    GO    TO    590 

no    510    L=l,NOVR 

SO=SQ+l 

IF     (SO    .GE.     100)    S0=1 

CALL    CLEAR 

VAL ( 1)=3 

TEMPd  J=0 

DO    515    AA^2,4 
515    VAL(AA)=-1 

TEMP(2 )=$W(J,L) 

TtMP(i)=$W( K,L) 

TEMP{4  )  =  $A( J,K) 

VAL(5) =-1 

TbMP(5)=$L( J) 

VAL(6) =1 

TEMP{6)=$L(K) 


N=6 

CALL    SPUNCH 
SQ=S0+1 

IF     (SQ    .GE.     100)    SU=l 
VAL(NI =-VAL(N) 
VAL{N-U=-VAL(N-1) 
CALL    SPUNCH 

IF     (    NCT5    .NE.    CP    )     GO    TO    510 
TEMP(2I     =    $V(J,LJ 
TEMPO)    =    $V(K,L) 
CALL    SPUNCH 
VAL(N)     =    -VAL(N) 
VAL(N-l)    =    -VAL(N-l) 
CALL    SPUNCH 
510    CONTINUE 

00    580     I=IJ,NOGT 

IF    (I     .EQ.J    .OR.     I     .EO.    K)    GO    TO    580 
SQ=SQ+1 

IF     (SQ    .GE.     100)    SQ=1 
CALL    CLEAR 
VAL(l)^3 
TEMP{1)=0 
DO    585    AA=2,^ 
585    VAL(AA)=-I 

TEHP{2)=$A(I , J) 
T£MP(3)=SAI J,K) 
TEMP(4)=$A( I,K) 
TEMP(5)=$L( J) 
VAL( 5) =-1 
TEMP(6)=$L(K) 
VAL(6)=1 
N=6 

CALL    SPUNCH 
SQ=SJ*1 

IF  ( SQ  .GE.  100)SQ=l 
VAL (N)=-VAL(N) 
VAL(N-l) =-VAL(N-l J 
CALL  SPUNCH 
580  CONTINUE 
590  CONTINUE 
600  CONTINUE 
CCCC   GATES  WHICH  AKE  CONNECTED  TO  ALL  OF  THE  OUTPUT  GATES  ARE  NOT 
CCCC   CONNECTED  TO  ANY  OTHER  GATES  CCCC 
CMT3=CL 

DO  650  I  =  'NJ0P1,N0GT 
SO=I 

CALL  CLEAR 
VALd)  =NOPT*U 
TEMPI  1 )=0 
DO  620  K=2,N0P1 
KM1=K-1 

TEMP(K)=$A( I ,KMl) 
VAL(K)=-U 
N=K 
620    CONTINUE 

DO  630  J^NOPl.NOGT 
IF  ( I  .EQ.  J)  GO  TO  630 
N=N+l 

TEMP(N)=$A( I, J) 
VAL(N)=-1 
630  CONTINUE 


CALL    SPUNCH 
650    CONTINUE 

CCCCCC^F^A^GATE    HAS    ONLY    ONE     INPUT.    THEN    IT    CAN    BE    SET    TO    A    NUK    GATECCC 

DO    700    I=1,N0GT 

SO=I 

CALL    CLEAR 

VAL(lJ=-2 

TEMP(l)=0 

VAL(2)=1 

TEMP(2)  =  $L(n 

N=2 

DO    660    L=1,N0VR 
N=N+I 
VAL(N) =1 
TEMP(N»  =  $W(I  ,L) 
660    CONTINUE 

IF     (     NCT5    .NE.    CP     )     GO    TO    667 
00    665    L=I tNOVR 
N=N  +  I 
VAL{N)=1 
TEMP(N)=$V(I ,L) 
66  5  CONTINUE 

667    CONTINUE 

DO    670    J=1,N0GT 
IF     (  I     .EO.    J)    GO    TO    670 
N=N*-1 
VAL{N) =1 
TtMP(N)=$A( J, I ) 
t70    CONTINUE 

CALL     SPUNCH 
700    CONTINUE 
710    CONTINUE 

CALL    CLEAR 

INPT(l)=E 

INPT(6)=-99 

N8  =  6 

CALL     INTLPINB) 

KE  TURN 

END 
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